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PREFACE 


This documentation is .prepared for the Goddard Modeling and Simulation 
Program (and others with some familiarity with general circulation modeling) 
to provide a detailed description of a grotind hydrology parameterization that 
can be used in the Atmospheric General Circulation Model (GCM) . The prepa- 
ration of this documentation is part of an ongoing program involving 
hydrologists with atmospheric modelers to develop realistic Interaction 
between the ground hydrology and the atmospheric components in the GCM. 

In 1973 the GISS-GCM was documented by Tsang and Kam (TK) . In this version 
the ground wetness is kept constant according to the particular monthly 
climatic surface relative humidity values for the globe. The ground wetness 
parameterizations were subsequently revised by GISS modelers to account for 
the effects of a variable ground wetness on the atmosphere. However, we 
have found a serious lack of information about how the parameterizations 
were developed. Our documentation is based primarily on the ground wet- 
ness and related sections of the computer coding (Subroutine C0MP35). 

This documentation is not intended to be self-contained. The reader 
should refer to TK for further details. 

We have developed an off-line testing program for diagnostic and 
sensitivity studies of the groxmd wetness and ground temperature parame- 
terization. As a result, some revisions and refinements of the Ground 
Hydrology Model (GHM) were incorporated in an off-line version, labeled 
as UCNI. 
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PART I 


MODEL DESCRIPTION 




I . Introduction. 


The ground hydrology portion of the Goddard GCM models the moisture 
transport processes related to the earth’s surface that were considered to 
be important to the GCM. The evaporative, sensible and ground heat fluxes 
that are calculated in- the ground hydrology are the principal- terrestrial 
output parameters for the GCM. The principal input parameters to the GHM 
from the GCM are precipitation, either in the form of rain or snow, and net 
solar radiation. In the hydrologic cycle the net change in ground cell 
moisture levels and surface water storage are generally small for long-term 
averages; then precipitation, evaporation and runoff are the major components 
of the moisture balance in the GHM. 

The present GHM parameterizations are conveniently divided into the 
following parts: 

a) ground wetness and soil parameters, 

b) precipitation, 

c) evapotranspiratlon, 

d) surface' storage of snow and ice, and 

e) runout (moisture which is unavailable for ' evaporation and for cell 
moisture changes) . 

In this documentation all of the above moisture transfer processes are 
considered. In addition, the effects of ground wetness on the ground 
temperature, the surface specific humidity and the surface temperature will 
be treated separately. The discussion focuses on understanding the computa- 
tional aspects of the GHM and how they affect the flux estimates. This 
documentation does not attempt to develop an extensive derivation or 
statement of assumptions for the parameterizations at this time. Although 



an attempt has been made to fully understand the current GHM parameteriza- 
tions we must point out that the presently available supporting docinnentation 
for this report was not completely adequate; therefore we made inferences 
from the GEIM computer code. 

Contrary to the traditional hydrology model which is based on the 
approach of drainage basin system, the development of the GHM parameteriza- 
tions are designed to capture, in a gross manner, the vertical fluxes for a ■ 
terrestrial cell 4 degrees latitude by 5 degrees longitude. This type of 
approach or modeling for hydrologic processes like evapotranspiration, 
runoff, infiltration, etc. is difficult because of the large size of the 
cell (which requires arbitrary lumping of hydrologic parameters) and the 
arbitrary divisions of the latitude and longitude slices that generally 
cut through drainage basins. The fundamental moisture coupling concept is 
based on the storage .capacity of the soil moisture near the surface to 
satisfy and therefore regulate the actual evaporative demand. The actual 
evapotranspiration depends on the characteristics of tine ground surface, 
soil moisture level and the potential atmospheric evaporative demand. 

The computer programs and the flow chart for the 'GHM are described in 
Sec. III. An off-line ground hydrology model (OLM) , devised for diagnostic 
and sensitivity studies, is discussed in Sec. IV. .Some results from numerical 
experiments are also presented in Sec. IV. The computer program description 
and flow chart for the OLM can be found in. Sec. V. Listings of .both -programs 
for the GHM and OLM are in Appendix B with symbols and notations described 
in Appendix A. 
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II. Formulation 


2- 1 ground Hydrology 


The present grovinS hydrology parameter! zations are developed for 
a ground layer of uniform thickness throughout the globe for each cell 
following essentially the approach outlined in Arakawa C1972). The funda- 
mental ground moisture concept is based on the availability of moisture 
in soil to satisfy the evaporative demand of the atmosphere either from 
bare soil surface or from vegetation canopies. • The equation of moisture con' 
servation is given 'by 
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where w' is the ground wetness in [fraction of saturation]’, 
is the density of water in [g cm"^], 

S is the porosity of soil, 

h is the deptli of the ground layer in [cm] , 

P is precipitation in [g cm”^ s~^]. 

Eg is evapotranspiration from the ground surface in [g cm“^ s“^], 

Rq is "runout" of water from the cell in [g cm~^ s“^], 

(This is not the conventional hydrologic definition of runoff (Sec 2.1) Runout)) 
and is the surface storage of snow or ice in [cm] . . • 

There are two additional conditions governing the ground wetness : 

(i) 0 < w' < 1 (2-1- lb) 




- (ii) w' = 1 , Sgi > 0 


(2-lrlc) 



Certain critical remarks concerning the fundamental concepc of the 
parameterization are necessary at this point. The present formulation ' 
of the GHM incorporates the folloirfing basic hydrologic assumptions: 

1. surface water storage other than snow and ice is negligible; 

2. the presence or absence of a ground water table is ambiguous; 

3. ground water storage or storage in the aerated zone below 
the bulk layer is neglected; 

4. horizontal moisture fluxes between adjacent ground cells - 
surface runoff, flow in the unsaturated zone, groundwater 
flow, deep vertical percolation, etc. -- can be eitKer lumped or 
calculated by a term in (2-1-la) , and 

5. the hydrology of individual ground cells are independent of 
each other, i.e. mass continuity of the horizontal moisture 
fluxes is not preserved between adjacent' cells . 

Another major concern of the parameter! rations for the GCM is how 
the actual evapotranspiration, E^, in (2-1-la) can be estimated, because 
Eg* also affects the ground temperature and the sensible heat flux ‘ 

(see Sec. 2.2 and 2.3) in addition to the moisture flux balance. 

Ground Wetness and Soil Parameters 

The ground wetness in the present parameterization is not clearly 
defined, causing a certain amount of confusion. Inferences had to be made 
from the specific values assigned to the soil parameters. The depth of 
the ground layer is not explicitly specified; instead, a maximum amount 
of water, w^ or WGfl, available for actual evapotranspiration is pre- 
scribed for two types of soil layer: WGM = 66.4 g for a "composite soil" 
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In addition, an impervious surface 


and WGM = 4.0 g for a "sandy soil." 
is identified by a "ROCK" factor. In reference to UCLA's model [Arakawa, 
1972], an amount of 10 g available for actual evapotranspiration 
was used. Therefore, it appears that WGM in this model is closer to the 
saturation moisture content than the UCLA's maximum available amount. If 
we interpret the WGM as the saturation moisture content, then the ground 
wetness, w' (WET), in this model is defined by volumetric moisture content 
as a fraction of saturation. We have adopted this definition throughout 
this documentation. ,Thus, w,„ = P^Sh in (2-1-la). For the "composite soil," 
Wjj = 66.4 g which is equivalent to specifying S = 0.52 and' h = 1.28 m. 
However, the maximum available moisture of 10 g in UCLA's model has been used 
as the range over which the actual evapotranspiration is scaled from the 
potential evaporation by a coefficient g. g is defined as the ratio of 
actual evapotranspiration to potential evaporation. For the "composite 
soil," the lower limit of this range is set in the neighborhood of the 
permanent wilting percentage, PIVP (USM) at a groqnd wetness, w’ , of 0,1, 
vifhile the upper limit has a value of 0.25 which is near field capacity, FC. 
The available amount of moisture for actual evapotranspiration (0 ^ g ^ 1) 
is then ASM = 0.15 x 66.4 = 10 g (=FC-PWP) [Salter and Williams, 1965]. The 
parameterization for evapotranspiration will be discussed later. 

For the "sandy soil", ASM = 0.065944 g and USM = 0.024278 g, which are 
consistent with the field capacity and permanent wilting percentage for a 
sandy loam [Salter and William, 1965]. The value of w^^^ = 4.0 g is not 
consistent with the parameterization outlined in the preceding paragraph.. 
Hoivever, if 4.0 g. is intended as the maximum available amount of moisture 
for a desert region, then the layer thickness should be adjusted- according 
.to the porosity of the "sandy soil." 


- 5 ' 



Another important aspect o£ ground wetness in the model roncems the 
transformation of the initial values of "ground wetness" from the climate data 
of Schutz and Gates (.19 72) . The ground wetness is first derived from the 
surface relative humidity data of Schutz and Gates according to 


GW 
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-'15 


85 


0 < GW < 1 


( 2 - 1 - 2 ) 


where r^ is the surface relative humidity and GW is the initial "ground wetness" 
input matrix in the'GCM. We cannot find the basis for this transformation. 

This "ground wetness" is coded as WAT, which is used only for computing the 
surface relative humidity and the criterion temperature (TK, p. 132). 

For the ground hydrology, another transformation is carried out according to 

WET = 0.5*C0,3*GW + 0.4) (2-1-3) 


which reduces the initial "ground wetness" from the range of (0,1) to 
(0.2, 0.35). Schutz and Gates' surface relative humidity may be considered 
as a long-term climatic average value which is proportional to the corres- 
ponding quantity near or on' the ground surface. It is well known that in 
moist soil the relative humidity in the soil remains high even if the soil 
moisture content reaches field capacity [Child, 1969]. This implies that the • 
atmospheric relative humidity near the ground becomes less than unity only 
when the ground wetness falls below the field capacity. Thus, it is reason- 
able to scale the range of the surface relative hxmiidity data of Schutz and 
Gates for the ground wetness between the field capacity and the permanent 
wilting percentage. Since the ground wetness in this model is considered 
as depth averaged value, the ground wetness on the ground surface could 
become less than the permanent wilting percentage during an extended drying 
cycle. Therefore, the range of the ground wetness scale should be .somewhat 
larger than the interval between the field capacity and the permanent wilting 
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percentage to realistically accomplish the transformation. The upper limit 
at WET = 0.35 is not unreasonable, but the lower limit at IVET = 0.20 should 
be reduced to a value at or below the permanent wilting percentage. 

Precipitation 

The precipitation is the primary moisture input in this 6HM in the 
form of rain or snow. The contribution of rain to the ground wetness is 
by means of the surface infiltration. There is no surface retention of .rain 
water allowed' in, the model. The surface runoff' and seepage are lumped 
together as the '’rimout." Snoiv is first accumulated on the ground surface 
as a form of surface storage. IVhether it .will contribute to ground wet- 
ness depends on the ground tenperature above or below the freezing tempera- 
ture. Therefore, the melting of snow will affect both the ground wetness 
and the ground temperature. 


Evapotranspiration 


The actual evapotranspiration in the GHM is scaled from the potential 
evaporation by means of a scaling factor, g. The potential evaporation, E 
(PREVAP) in [mm/hr] , is given by 


P 


^p ~ '^aS 


W (q* - 
s'-^g 




C-2-1-4) 


where p is the density of air in [gcm"^], 

Cp is the drag coefficient, 

Wg is the surface wind speed in (ms“^.), 

is the, surface specific humidity at a reference height above the 
■ ground surface, and 

q* is the ground saturation specific humidity. 
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The ground saturation specific humidity is a function of the. average 
ground temperature in' the diurnal layer. The actual evapotransporation. 
Eg (AET) , is scaled down form (2-1-4) according to 

E = BE 0 < B < 1. . (2-1-5) 

3 P 


For ocean, snow and ice, B = 1. For land, B is a function of ground 
wetness, soil characteristics, solar angle and the potential evaporation 
as the following: 

e- Vant" 

where B ^ ^ = 1. - exp { [-0. l*VPM*exp (16. 75*VPM)/PREVAP] *PLANT> 

plant (2-1-7) 


and 


^soil “ exp[-0.01*VSM*exp (16. 75*VSM)/PREVAP] , 

VSM = 2.*WET*(ASM + USM) . 


VPM = 0.i*VSM/USM + 0.2*(— — ) 


( 2 - 1 - 8 ) 

(2-1-9) 

( 2 - 1 - 10 ) 


It has been suggested that the parameterizations of (2-1-7) to (2-1-10) 
are' based on Van Bavel (1966) field study. However we are not able to pin 
down these relationships quantitatively. Therefore, a study of the range of 
variation of the above parameters was carried out as follows. The PLANT 
factor given as a function of solar angle, Z, 

PLANT = min (1.0, (2-1-11) 

is shown in Figure 1. It is interesting to note that for the 32° N Latitude, 
the ratio of daily maximum solar radiation at the ground between Dec. 21 and 
June 21 is 0.4 [Koberg, 1962] while the ratio of average PLANT factor derived 
from Figure 1 is 0.67.' Similarly, the ratios for April 21 and August 21 
were found to be 0.85 and 0.95, respectively. The comparison indicates that 
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the PLANT factor based on the solar angle is overestimated based on Koberg. 

For specified values of Ep and PLANT, the scaling factor 3 given in (2-1-6) 
is plotted in Figures 2(a) and (b) for the composite and sandy soil. These 
figures indicate that 3p2ant dominates the actual evapotranspiration scaling 
factor. The intent' of the specified values of field capacity and permanent 
wilting percentage in relation to 3 is not evident from the coded, parameteriza- 
tions. In an attempt to understand the effect of field capacity and permanent 
wilting percentage values we compared 3 with the parameterization of agricultural 
crops and a soil [Monteith, 1975] . This is presented in Figure 3 which shows 
that the current parameterizations represent a revision of the UCLA values and 
do not follow conventional agricultural type variations. For Ep = 0.1 mm/hr 
or approximately 1.2 mm/day if transpiration is shut off during the night, 

3 in the current model comes closest to the linear variation of the UCLA 
scaling factor between FC and PWP explained in Arakawa [1972]. 

In summary we see that the scaling factor for the calculation of actual 
evapotranspiration is a function of solar angle, Ep, soil type and moisture 
content and that in a gross sense the family of curves presented in the examples 
of Figure 3 is consistent with the previous UCLA formulation. However, there 
exists an inconsistency for the scaling factor in relation to the field capacity 
and permanent wilting percentage. If a parameterization embodies pertinent 
soil characteristics such as field capacity and permanent wilting percentage, 
the scaling factor should approximate the simple rule: a) once ground wet- 

ness has exceeded field capacity the actual evapotranspiration rate should 
approach the potential rate and b) at or below field capacity the scaling 
factor, 6, should fall below 1 and approach zero as the permanent wilting per- 
centage is reached. The parameterization in this model does not follow this rule, 
tending to scale down the potential rate excessively. If realistic evaporation 
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fluxes from the model are to be expected the potential rate would have to 
be significantly higher to compensate for the lower values of the scaling fac- 
tor g. On the other hand, the value of g is limited by this potential rate 
in the exponent of (2-1-7) and (2-1-8). Therefore, it appears that the 
actual evaporation rate is in general underestimated in this model. This 
may, significantly affect the calculation of moisture fluxes as an input to the 
GCM. See Section IV for furt.her details. 

Surface Storage of Snow and Ice 

The surface storage in the GHM is in the form of snow, ice or mixed ice 
and water. ' There is no surface storage of rain water on the ground surface 
during or after a rainfall. Thus the absence of surface storage inherently 
affects the computed values of ground wetness. The storage affects not only 
the ground wetness but also the ground temperature through transfer of 
latent heat. The thickness of snow and ice and the fraction of ice in 
mixed ice and water are updated every time step. ■ The surface storage in 
the ground wetness equation (2-1-1) and the latent heat transfer term in 
the ground temperature equation (see Sec. 2.2) are adjusted simultaneously. 
'Again the "ROCK" factor characterizing the impervious surface is appropriately 
incorporated in the adjustment of the transfer processes. 
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Runout 


The term, ’’runout," in the present GHM has a special meaning which 
is different from the "runoff" conventionally used in hydrology. Runoff 
is generally considered to be the continuous gravity jnovement of water 
on the ground surface and in channels with a characteristic magnitude 
which should be consistent with the intensity of precipitation, infiltra- 
tion and topography of the ground surface [Ward, 1975] The runout in 
this model is divided into two regimes, depending on the relative magnitude 
of the rainfall rate and the maximum surface infiltration rate. . Runout 
is defined as, in a gross sense, the sum of ail outgoing moisture fluxes other 
than evapotranspiration from a ground cell. In view of the assumption of 
hydrological iy isolated ground ceils, the runout moisture neither enters 
the surrounding ceils nor is being stored in a lower ground layer. There- 
fore, the amount of runout estimated in the model represents the portion 
of rain or melting snow which cannot contribute to the ground wetness and 
evapotranspiration. Runout consists of (a) surface runoff, (b) percolation 
and (c) any excess of moisture above the capacity of the soil layer (w* > 1). 

Before calculating the runout, an infiltration rate is first calculated 
as one half of 

BBFACT ~ 2. *FILTF*EXP(-1.75*(ADINF-0.5]3 in [cm/hh] (2-1-12) 


where 


FILTF = [1. 77. *EXP {-0.036* (TG - TICE)} + 0.44562]“^ (2-1-13) 


and 


ADINF = 1.25*(2 .*raT--0.2)-0.5 


(2-1-14) 



where TG and TICE are the ground and freezing temperature. FlilF is a 
temperature dependent parameter as shown in Figure 4 and ADINF is a function 
of ground wetness . Both ADINF and BBFACT as a function of ground wetness 
are shown in Figure 5. However^ the physical implications of ADINF and 
FILTF are not known at this time. IVhen.the rainfall rate exceeds twice 
the maximum surface infiltration rate (MSIR) , the excess contributes to 
the runout, or 

RUNOFF = RAIN - ,5*BBFACT In [cm/hhO ' -'(2-1-153 

When the rainfall rate is less, than twice the BBFACT', the runout is 
calculated by 

RUNOFF =.5* (RAIN) ^/BBFACT in [cm/hh] ■ (2-1-16) 

These formulas are shown graphically in Figures 6' and 7. Again, the hydrologic 
basis for (2-1-15) and (2-1-16) is not known to us. 

Given the ground temperature and soil moisture values one can first 
estimate the BBFACT from Figure 5. Tben, for a given rain input one can 
obtain the runout by using Figure 6. Soil characteristics or vegetal cover 
effects are not considered in obtaining runout, values. What are the charac- 
teristics values of runout? The following’ numerical example will help answer 
this question. Since typical global soil moisture contents are generally 
less than field capacity, WET = .25 should be a reasonable value for a 
field capacity or the upper limit of mean soil moisture values. Taking a 
typical ground temperature, say 20'’C, one finds, from Figure 5 a value for the 
BBFACT of 4.75. In Figure 6 we see the following for BBFACT =4.75: 
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RAIN 

RUNOFF 

[cm/hh] 

[cm/hh] 

1.0. 

,1 

2.0 

.4 

3 0 

.8 

4.0 

1.5 


A review of the computer output from the GOl indicates that the rainfall 
rates rarely exceed 1 [cra/hh] . The above table implies that runout is 
generally a small part of rain or that the major portion of rain con- 
tributes to cell moisture changes. However, a recent study indicates 
that, on a global scale, the ratio of runoff (which is a part of runout 
in GHM) to rainfall over the continent is 0.36 [Baumgartner, 1975]. 

Finally, runout versus soil moisture for a given temperature is plotted 
for various rain values in Figure 7. In the model, the runout is 
adjusted for the percentage of ROCK on the surface. The effect of a 
nonzero value of ROCK is to increase runout which reduces the surface 
infiltration rate. 

In order to balance (2-1-la) the runout is computed according to 
(2-1-12) through (2-1-16) . Two corrections of runout from the cell are 
made by (1) limiting the ground wetness to be less than unity as (2-1-lb) 
and (2) percolation by reducing the ground wetness in an amount of 5% if 
the ground wetness exceeds half of its maximum value. 

In summary, the runout parameterizations are a function of rain, cell 
moisture content and temperature. This represents a change over the previous 
UCLA model which used a linear scaling factor of the cell moisture to obtain 
a value of runoff. From the numerical example presented above it has been 
shown that runout is a small part of rainfall rate and hence plays a minor 
role in the calculation of cell moisture changes in this model. However, 
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Baumgartner has shown that runoff, which is less than runout, is not negligible 
hence the GCM calculations of cell moisture are questionable. 


2.2 Ground Temperature 

The ground .temperature in the G(M has been documented in TK. As 
the ground temperature is intrinsically coupled with the ground wetness, we 
intend to discuss only the terms in the ground temperature equation involving 
the ground wetness . The ground temperature is given by 
3T 

■ . C2-2-1) 


where T is the ground ten^ierature in the diurnal ground layer [Arakawa and 
Mintz, 1974] in [“K], 

C is the total heat capacity in [cal cm“^ deg”^] , 

Rj^ is the net radiation flux at the ground surface positive downward 
in [ly day“i], 

Fg- is the upward flux of sensible heat in [ly day“^]-, 

LEg is the upward flux of latent heat in [ly day"^], 

% is the heating or cooling due to conduction through sea ice from 
the ocean below in [ly day~^] , and 

is the transfer of latent heat due to melting of ice or snow or due 
. to freezing of water in [ly day"^] . 

The ground temperature equation is discretized with respect to time by an 
implicit scheme and, is updated in- the time step before updating the ground 
wetness. The finite difference equation is given in TK without the term, 

S^, as the following: 

,n 

“ . [ 2 - 2 - 2 ] 


T = T " . 


•Hjlt [IV,- 


LG +0 + S. ] 
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8R-, 8F 3E 30 

r + [Jl + -—L+ I -i- _ 

LaT a'T “ a'T aT -I 


■■aT 9T 
g . g 


3T 3T 
g g 


- 22 - 



where the superscripts n+1 and n indicate the time step and nj At is the 
time interval used in the subroutine. 

9F 9E 

s s 

In Equation (2-2-2), the terras, LE^, ‘3 ^" j 

. g g 

involve ground wetness. The' latent heat flux, LE^, or the actual evapo- 
transpiration, E^j.has been discussed in Sec. 2.1. The sensible heat flux, 
F_, will- be discussed in Sec. 2.3, because both the ground temperature and 
the surface temperature are affected by ground wetness. The total heat 
capacity is defined as 

C = C-^ ) in [cal cm“2 deg”^] (2-2-3) 

where c is the specific heat capacity in [cal cm'^deg ^], 

X is the termal conductivity in [cal cm“^ s“^ deg"^] ^ and 
• oj is the diurnal frequency in [s"^). 


values of C 

are given as the following: 


for ocean 

•C = » . 

(2-2-4ai) 

for snow, 

C = 2.3 . 

(2-2-4b) 

for ice. 

C = 5.1 . 

(2-2-4C) 

for frost 

9 



1/2 

C = [(0.331 + 0.075w*)(2; + 2.5w') * 0.001 * — ] (2-2-4d) 


for mixed water and ice, 

C = [(0.276 + (0.11 + 0.15w')(w*- 0.5w.))(l + w’ -0.25w^) 

* .002 * 24*5600] • C2^2-4e) 

27T 

where = the portion o£ ice in w’ = (WET)* 
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0 is replaced 

(2-2-4f) 

where VSM = 2 . 0*WET* (ASM+USM) , C2-1-9) 

ASM = available soil moisture in [% of saturation lOO] and 
USM = unavailable soil moisture in [% .of saturation -i- 100]. 

The proper use of C2-2-4f) requires an appropriate definition of WET. This 
has been discussed in Sec. 2.1. Therefore, WET in (2-2-4f) is the ground 
wetness defined in this model despite the fact that the maximum available 
amount of moisture used in this model is different from UCLA's value. 

The transfer of latent heat, Sj^, due to melting of ice or snow or due' 
to freezing of water is‘ updated at each time step according to the ground 
temperature below or above the freezing temperature. 

2.3 Surface Conditions In the Atmosphere 

The formulation of the atmosphere- ground hydrology interaction depends 
on parameterization. of both the planetary boundary layer (PEL) and ground 

hydrology. In TK, the boundary layer is taken .to be the bottom half of the 

( 

ninth layer in the GCM and the location of the "surface," the bottom of the 
ninth layer, with respect to the ground surface is not specified. We shall 
follow the same approach in this documentation. The atmospheric conditions 
that are coupled with the ground hydrology are related to the- source terms 
in the governing equations of the GCM (see TK, pp. 7-13); namely, 

(1) the surface wind and stability conditions 

# 

(2) the heat flux 

(3) the moisture flux, and 

(4) the solar radiation. 


For soil, the original UCLA formula, (2-2-4e) , with Wj|. = 
in the present model by 

C = [(.3 + VSM) * (1+ 10. *VSM) *0.001 * 
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Most of the discussions in the folloiving paragraphs can be referred 
to TK (Chapter XI, pp. 109-191). 

Suffice Wind and Stability Conditions 
Since the lower surface of the GCM or PBL is not explicitly specified, 
we assume that the '^surface" is located at the level of anemometer height, 
say, S ni from the ground surface. The surface wind components are extrap- 
olated from the values at the n-th and n-1 th layers with respect to their 
a coordinates by 


u - u 1 u -u 
n n-1 _ s n 

0 -a ~ a -o 
n n-1 s n 


(2-3-1) 


th 


where the subscripts n and s refer the quantities to the center of the n 
layer and to the "surface*' and u is the x-component of wind velocity. The 
y-component of velocity, v^, is computed in a similar manner. The magnitude 
of the surface wind velocity is then 


W = (U ^ + V 2)1/2 [jji g-1] 


(2-3-2) 


Special provisions are made for computing wind speeds at the north and south 

poles . The surface wind affects not only the surface skin friction but also 

the moisture and heat flvixes on the ground surface. These fluxes depend 

also on the drag coefficient, Cj^. For a neutrally stable atmosphere when 

the surface temperature, T , and ground temperature, T , are equal, 

s ^ 

Cy f Cp = min[(1.0 + 0.07 Wg2)/io00., 0.0025], for ocean (2-3-2a) 
o 

= 0.002 + 0.006(3^/5000), for not ocean (2-3-2b) 


where Zg is the ground surface elevation in [m] . 




In general. 


Cj. = /I1-7.0CT - T )/ W ^] , T < T 

^ Dq' ^ '■ g - s-^ ' S ' g s 

* Cf, • [1,+ /CT'“-'T^) /W , T > T . 

D '• ^ g s s ' g — • s 


(2- 3- 3a) 
(2- 3- 3b) 


These coefficients will be used in the moisture and heat flux computations. 
Heat Fltut: and Surface Temperature 

The- sensible heat flux on the ground surface is computed according to 


Pa fly day"M 


g 


(2-3-4) 


which is coupled to the heat flux in the PBL by the turbulent diffusion 
equation 


‘r r 1 / 

n ,P ® 3Z 


(2-3-5) 


PBL 


where Pg is the air density in [g cin"^]^ 

Cp is the specific heat at constant temperature in [cal g~^ deg"^], and 
K iS' the vertical eddy diffusivity in [m^ s“^]. 

The horizontal convergence in the boundary layer and geopotential difference 
between its top and bottom are neglected, then the surface temperature, T^, 
can be computed by equating the above two equations as described in detail 


in TK. However, we would like to point out one computation directly involv- 

9T I 

ing the ground -wetness for the.- temperature gradient, 


PBL 


, in (2-3-5) 


The temperature gradient is given by 


9T 


^^TB “ *^s^ " ^^TB ‘ ^s^crit 


PBL 


THJ 


(2-3-6) 


PBL 


('^TB - '^s^crit = . ^s^'V^s^nua. ^^-3-7) 
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where (^z)pgj^ is the thickness of the PBL and the subscripts, TB, s, Cd.a.) 

and Cm. a.) refer to the top of PBL, the "surface" or the bottom of PBL, dry- 

adiabatic and moist adiabatic, respectively. The relative humidity on the , 

"surface", r , is. derived from 
s 

2 w * X 

r = , and w' = WET (2-3-8) 

s w' + r„ 

where w' is the ground .wetness (see TK, p. 132). It appears that -the 
ground wetness is treated like the relative humidity on, the ground surface in 
(2-3-8) . The relation between the ground wetness and the relative humidity 
in a soil layer has been discussed in Sec. 2.1. In order to use (2-3-8), 
the following transformations are made: 

-WAT = 2. *(1. 66667*2. *WET - 0.66667) ' ’ (2-3-9a) 

for the "composite soil", and 

WAT = 0.2*(2.*WET - 0.26909) . (2-3-9b) 

for the "sandy , soil." WAT is used as w' in (2-3-8) instead of WET. In a 
way, it is an inverse transformation of (2-1-3) , which changes the ground 
wetness back to the "initial ground wetness" (Schutz and Gates* surface 
relative humidity) . 

On the other hand, the surface specific humidity is computed by another 
method using (2-3-10) as discussed in the following section. 

Moisture Flux ■ 

, Ihe moisture flux from the ground is equal to the actual evapotranspira- 
tiqn given by (2-1-4) and (2-1-5). In these equations, the surface specific 
humidity, q^, is derived by equating the flux in the PBL to the fiux from the 
surface as the followings 
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(2-3-10) 


q 


s 


. EVACO*q + EDV*q^ 
“ EVACO + EDV 


where q* and q^ are the saturation specific humidity at the ground surface 
and the specific humidity at the 9^ level in, the GCM. EVACO = an 

evapotranspiration coeffient in (2-1-4) . and • (2-1-5) . EDV = K/CAz)^^^^ , 
a surface layer eddy coefficient. Physically, (2-3-10) i^lies a constant 
moisture flux through the PBL similar to the assumption used in deriving 
the surface; temperature. However, in case of supersaturation, the surface 
temperature and specific humidity are adjusted to account for the heating 
effect of condensation and the releasing of moisture (see TK, pp. 135-136) . 

Solar Radiation 

The radiation model in the GCM has not been documented at this time; 
therefore, it must be referred to the forthcoming docimentation. In run- 
ning the GHM on-line with the GCM or off-line alone, the radiation terras 
in the ground temperature equation (2-2-1) are treated either as’ the version 
given in TK or as input quantities from the GCM output. 


- 28 - 



Ill . Description of the Computer Program for the GHM 
Subroutine COMP 35 

The ground hydrology model resides in subroutine C0MP35 in the 
GCM computer program. The function of C0MP35 is to evaluate the hydrology 
related source terms in the governing equations' and their effects on the 
atmospheric motion. It should be noted that COMP 35 is different from COMPS in 
TK in that the ground hydrology parameter! zations were later added. There 
are 16 steps in C0MP35 as in C0MP3 as illustrated in the flow chart. Steps 
i, 2, 3, 6 and 12 interact directly with the GUM. The major ground 
hydrology parameter! zations are incorporated in Steps 7 and 14. Step 7 
deals with the coupled effects of the ground hydrology and the atmosphere 
on the "surface" temperature and humidity as discussed in Sec. 2.3 and the 
.parameterizations of evapotranspiration for the latent heat flux. In Step 
14, the ground temperature is first updated according to the scheme outlined 
in Sec. 2.2 and then ground wetness change is calculated on the basis of 
C2-l-ij, 
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FLOW CHARTS 


C0MP35: 


COMP 35 


HYDOOOlO 

HYDOmO 


Step T 

Initialization of Radiation 
and Heating Parameters 


COMP 3 


Step 2" 

Surface Wind Magnitude 


COMPS 


Step. 3 


Surface Conditions 


Step 3A 

HYD01120. 

Surface Topography 

HYD01620 

Step 3B 

HYD01630 

Ground Temperature and 
Wetness 

HYD01670 • 

Step 3C 

HYD02090 

Surface Drag Coefficient 

. HYD02160 









Steo 4 


Set Up Vertical Arrays for 
Sotae Variables 


Step 4A 

HYD01680 

Pressures 

HYD01770 

Step 4B 

HYD01780 

Temperature & Convective Adjustment 

HYD01930 

j 

Step 4C 

HYD01940 

Moisture Variables 

HYD02000 • 


Sten 5 

HYD02170 

Vertical Diffusion of Heat & Moisture 

HYD02420 


Step 6 

HYD02430 

.Saturation Specific Humidity 


and Specific Humidity 

HYD02480 


Step 7 

HYD02490 

Determination of Surface 
Temperature and Specific 
Humidity 

HYD03960 










step 8 


Contributions of Source Terms 
at Level Nlay 

• 

Step ' 8A 

HYD03970 

Sensible Heat & Evaporation 

HYD04160 

Step 8B 

• HYD04170 

Surface Skin Friction 

HYD04450 




Steo 9 

I * 

* 

Moist Convection 

' 

Step 9A 


Preparation for Moist Convection 

COMP 3 

Steo 9B . 

* 

Vertical Advection of Moisture 

C0MP3 ' • 

Step 9C 

■ HYD02010 

Constrain- Specific Humidity 

. HYD02080 

Step 9D 


Special Moist Convection, 
for 9 Level Model 

•C0MP3 

Step 9E 


Types of Convection 

C0MP3 

Step 9F 


Distribution of CVTP and 
CVQP 

COMPS 


step 10 



HYD07330 
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C0MP3S, 
STEP 7: 


Step 7 


Determination of Surface Temperature 
and Specific Humidity 


HYD02490 
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COMP 35, 
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Flux Terms for Diagnotic 
Purposes 


HYD05860 

HYD05940 



2650 



















Snow Melt 


HYD05280 

HYD06310 



HYD06320 


Ground Wetness 

HYD0638Q 

for Snow Covered Land 

HYD06450 


Ground Wetness for 
Snow Free Land 


HYD06490 

HYD06680 


Ground Temnerature 
Update 


HYD06690 

HYD06820 


Ground Wetness 


HYD06870 

HYD06910 
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Ground Wetness Packing 

HYD06920 

HYD06960 



End of Step 14 

HYD06970 

HYD07330 
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PART- II. 


MODEL RESULTS & ANALYSIS 



IV,' Off Line Model (0L7t) 


Introduction 

The OLM provides a simple and efficient way to test various hydrologic 
parameterizations without the need for extensive simulation time. The hydro- 
logic parameter! zations are assembled in a parallel manner in the OLM and 
given the following identification: 

1. UCNI - a first revision by the University of Connecticut, 

2. UCLA - the UCLA Model (Arakawa, 1972; Arakawa & Mintz, 1974), 

3. GISSJK ~ Jack Komfield's work previously referred to as GHM, and 
discussed in this text. 

UCNI was developed because we felt that a hydro logically consistent model 
using only simple and well-defined parameterizations was needed for evaluating 
the existing GHM. UCNI represents our first step in developing a physically 
realistic multi-layered model. Specifically, the following ideas were incor- 
porated into UCNI: 

a) soil moisture is defined in fraction of saturation; 

b) cell moisture changes are calculated from (2-1- la) with respect to 
the amount of water at saturation; 

c) saturation moisture content is a function of given layer depth and 
porosity of soil type; 

d) the initial soil moisture range is linearly scaled between field 
capacity (FC) if the surface relative humidity (r^) equals 100% to permanent 
wilting point if rg equals 0; 

e) available soil moisture is the difference between FC and PWP; 

f) evaporation scaling factor, B, varies from unity at FC to 0 near PWP, 
with plant factor (2-i-ll) kept constant at unity regardless of solar angle; and 
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gD vegetation and soil type may change effective soil layer thickness . 
The OLM reads an output or data tape from a previously run GCM experiment 
with a certain known ground hydrology model. The ground hydrology parameteri- 
zation generally tised in GCM experiments to develop the data tape has been 
described in TK. Recalling these parameter! zations we see the following: 

1. $ = Min(l,2w0, 

2. . 0 < w' < 1, 

3. WGM = 10 gj and 

4. . w' = initial w' for all time = WET. 

The data set the OLM reads every half hour from the GCM data tape is: 

1. TAU, the model time in [hr]; 

2. PC, the surface pressure in [mb]; 

3. TSC, the surface temperature in [°K]; 

4. SlISC, surface specific humidity [g/g] for first time step the 
net radiation in [ly/day] for 2 through n time stepsj 

5. GTC, the ground temperature in [°K]; 

6. GWC, the ground wetness in [% saturation ^ 100]; 

7. UC, the X component of velocity in [m/sec] for levels 7, 8, 9; 

8. VC, the y component of velocity in [m/sec] for levels 7, 8, 9; 

.9. TC, the temperature of air in [°K] for levels 7, 8, 9 ; and 

10. SHC, the specific humidity in [g] of water vapor/[g] of air for 

levels 7, 8, 9 . 

The ground hydrology model being evaluated calculates the following diag- 
nos tic p arame te rs : 
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1. GW, the ground wetness in [% saturation v 100]; 

2. TS, the surface ten^jerature in [°K] ; 

3. GT, the ground teH 5 >erature in [°K]; and 

4. AET, (formerly E6DT) the evapotranspiration flux in [cm/hhj . 

With the GCM data tape and subroutines SDET (Solar angle calculation) , 

SURFWD (Surface wind calculation) , and RITE (Diagnostic parameters) the OLM 
simulates an on-line run and calculates the diagnostic variables. At the 
end of the experiment comparisons between GW, GWC, TS, TSC, GT, GTC and AET, 
AETC (from UCLA) were made and will be discussed in a later section. 

The OLM-UCLAl and 0LM-UCLA2 were incorporated in OLM to calculate the 
AET flux (actual evapotranspiration, missing on data tape) for the case 
of variable ground wetness (GW?5K, OLM - UCLAl) and constant ground wetness 
(GW=K, OLM - UCLA2) . The recalculation of AETC was also necessary because 
of the different eddy diffusivity parameterizations used in developing the 
GCM data tape versus OLM. The TK eddy diffusivity paramterization was used 
In off-line experiments because it required a substantially smaller data set. 

In figure 4-1 a schematic is presented of the parallel hydrology param- 
eterizations in the OLM. The accuracy of the OLM in simulating an on-line 
calculation can be evaluated by comparing UCLA(GW = K) with OLM-UCLAl or OLM- 
UCLA2. Close agreement between UCLA and OLM-UCLA2 was always found. In 
particular. Table 4-1 outlines how diagnostic variable comparisons are made 
for the different hydrologic parameterizations. 
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GCM with UCLA Formulation and NCAR Eddy Diffusivity (EDNS) 


2 day run 

I 


GCM Data Tape 



Reads and utilizes 
TS, GT, GW for 


Reads and utilizes 
all data sets . 


first two time 
steps only (LOGS: 
false) 

4 

UCNI 

OLM-UCLAl (GW i K) 
OLM- UCLA-2 (GW = K) 
GISSJK 


for all time steps 
(LOGS: true) 


UCLA (GW=K) = OLM - 
UCLA2 (GW = K) 

OLM-UCLAl (GW / K) 


FIGURE 4-1 Schematic of Parallel Model in OLM 



TABLE 4-1 


Model Diagnostic Parameter Read and Calculation Scheme 



R: Model reads data tape variable every time step 

C: Model calculation of this variable will be used for diagnostic 

puipos es 
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It should be noted that the surface and ground temperature are calculated 
by the following different schemes: 

TSC^ ->• ,, ->-TS^'^^...etc. 

GTC^ ^ ■* GT^‘*’^...etc. 

Thus, TS differences reflect variations during a time step while GT differences 
represent cumulative effects since second data tape read. 

Differences in Hydrologic Parameterizations 

In order to understand the significance of the coii^arisons discussed 
later we .will highlight their differences in this section. First, the UCNI 
differs, from the UCLA, parameterizations in the following areas: 

a) initialization of GW, 

b) - approach taken in scaling 0 as a function of GW, 

c) soil water content at saturation, and 

d) definitions of GW and how it changes . 

Secondly, the GISSJK parameterization differs from UCNI in the same areas, a) 
through- d) plus the following: 

e) plant factor (effect of solar radiation on 0), 

f) volumetric heat capacity, and 

g) runout . 

Among these differences, the principal factors that lead to significant 
variations in the' estimation of the evaporative flux were 0 and different 
definitions of changes in GW. Table 4-2 stiramarizes the differences between 
these hydrologic parameterizations (UCNI, UCLA, GISSJK). 
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TABLE 4-2 Comparisons of UCNI , UCLA 


UCNI 


1. Layer 

CS 

100 cm 

Depth 

■SS 

50 • cm 

2. Porosity 

CS 

.5 


SS 

.5 

3. Field Capacity 

CS 

■ .25 ' 

[% Saturation/ 100] 

SS 

.09 

4. PWP 

CS 

.10 

[% Saturation/ 100] 

SS 

.025 

5. Soil Water Content 

CS 

50 

at Saturation 
[g/cm^] 

SS 

25 


6. Available Water 
[g/cm^] 

7. GW Initialization 

8. Plant C^l) 

9. (Plant = 1) 

10 . C (GW=0) 

[cal (!eg“^] (GW=1) 

(GW=13 


CS 

7,5 

SS 

1.6 

CS 

SRH (FC-PWP) + PWP 


100 

SS 

t! 


[% Saturation/ 100] 

cs/ss 

-- 


WET- PWP 

cs/ss 

FC-PWP 

CS 

Cl 

CS 

C2 

SS 

C2 


GISSJK Parameterization 




TABLE 4-2 Comparisons of UCNI, UCLA and 6ISSJK Parameterization (cont.'d) ■ 




UCNI 

UCLA 

GISSJK 

11. Runoff [g/cm^] 
Examples: 

Rain( 

.1 + .9(GW)^-5) 

RainC.l + ,9(GW)l-5) 

See Text 

Rain = .1 [cm/hh] 

GW=0 

.01 

.01 

0 


GW=.5 

.042 

.042 

.05 


GW=1.0 

.1 

.1 

.1 

Rain = 1.0 [cm/hh] 

GW=0 

.1 

.1 

0 ■ 


GW=.5 

..42 

,42 

.45 


GW=1.0 

1.0 

1.0 

.9 ' 

12. EDV 

[m sec“^] 


TK (1973) 

NCAR • 

TK (1973) 


Notations: CS: composite soil; SS: sandy soil 

SRH: surface relative humidity 

EDV: surface layer eddy coefficient 



Experiments 

To conduct the OLM experiments a 2 day GCM run (DSNAMES: I719/VRM0, 

VRH0, VEP0, VRR0, 11/76) was executed from which the input GCM data tape 
variables, discussed in the introduction, were extracted for the general US 
surface and neighboring ocean cells for levels 7, 8, 9 of the GCM. The ground 
surface and ocean cells are shown in Fig. 4-2. . This compacted or condensed 
GCM tape information was then read by the-OLM for the following cases: 





' US 

Globe 

Area 


1. 

Dry 

(J>I) 

(4,4) 

(J,l) 

(33,15) 

Dry Southwest 

Cell Types. • 

2. 

Average 

(3,6) 

(32,17) 

Southern Great Plains 


3. 

Wet 

(7,12) 

(36,23) 

Canada/Wew England 


V 


in combination with the following ground hydrology parameterizatlons: 

UCNI 

OLM-UCLAl CGW/K constant) 

0LM-UCLA2 (GW=K constant) 

UCLA (GW=K constant) 

GISSJK 

These experiments took about 1 minute on IBM 360-50 at the University of 
Connecticut Research Computer Center while the Amdahl at GSFC took about 
half this time to execute. 

Results 

The OLM provided a simple comparison of the variolas ground hydrology 
parameterizatlons that were reviewed in this section. In particular, the 
following parameters are discussed and summarized: 


- 51 - 



Soil moisture 


Surface, and ground temperature 

Potential § actual evapotranspiration and their ratio, g. 

Soil Moisture 

The following comparisons were made about soil moisture levels and 

> 

variations for various cell cases (see Fig, 4-3 to 4-5) : 

1. GISSJK moisture levels were significantly higher for all 
three moisture cases; 

2. UCNI and UCLA moisture levels were close in all three 
moisture cases; 

3. OLM-UCLAl (GWj^K) moisture levels changed significantly 
for typical ET and rain fluxes; and 

4. In general the change of moisture levels was not significant 
during two day test period. (Therefore longer OLM runs of 
duration greater than 2 days should be made.) 

Surface and Ground Temperature 

For temperature comparisons the following assessments were made (see 

Fig, 4-6 to 4-11) : 

1. The ground and surface temperature calculations for the UCLA 
' model and UCNI were consistent for all cases with better . 

surface temperature agreement during daylight hours-; 

\ 

2. Temperature predictions were generally in phase for most 
models and cell types with the exception of GISSJK; 

3. The ground and surface temperature were generally higher 
for GISSJK, sometimes as much as 10®C higher for- the 
different cell types; and 
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A. The, ground temperature for GISSJK demonstrated a phase . 

shift of about 3 hours with a gradual upward drift in magni- 
tude from those of other models .■ 

Potential .(ETj and Actual (AET) Evapotranspiration and Their Ratio, 3 . 

The OLM provided the following results related to evaporation rates 

I 

for the three parameterizations (UCNI, UCLA, GISSJK) (see Fig. 4-12 to 

4-14 and Table ’4-3) : 

I; ■ Potential evaporation rates were generally highest for 
GISSJK for all cell types; 

2. The potential evaporation scaling factor, 3, was lowest 
for GISSJK and highest for UCLA for all cell types; 

3. The actual evaporation rate was lowest for GISSJK for 
all cell types; and 

4. The UCNI generally predicted values of PE and 3 that 
appear reasonable when compared with GISSJK and UCLA 
models . 


Summary 

'fhe major findings discussed in previous sections individually will 
now be summarized collectively; 

Different .definitions of soil moisture, although significant 
between models-, did not cause major differences in the calculation of AET. 
This suggests a certain consistency within each model. 

■ Surface and ground temperature results are similar for all 
models for the range of moisture levels tested with the exception of GISSJK. 

Ground temperature difference of +I0'’C between GISSJK and other 
models is believed to be due to the underestimation of AET by a factor of 2 
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or more. 


Surface temperature variation of 10 between time steps had 
appeared to be due, to separate (stable, unstable) parameterizations for the- 
surface layer eddy coefficient* The variations had a deleterious effect 
on* AET when they occurred. When the OLM used the eddy dif fusivity in TK, 
the oscillations were greatly reduced. 

Peak AET seems to occur when surface and ground temperatures 
peak' at mid-day. 

3 as a function of PET&GW in GISSJK led to a ”flat’\- fairly 
constant AET for. most of the daylight hours. 3 as a function of GW in 
UCNI & UCLA produced significant variations of AET during daylight hours. 

In OLM-UCLAl with, variable ground wetness noticeable changes in 3' occurred 
during significant evaporation or precipitation periods. 

‘AET estimates for three different soil moisture levels, wet, 
average. and dry case, were remarkably similar in magnitude for the 2 days. 
This is due to the compensating factors of a large PET when 3 is small for 
the dry case or vice versa for the wet case. 

Small negative evaporation rates generally occur during evening 

hours. 


UCNI appears to be. a reasonable revision ^of currently available 
ground hydrology models. 

In conclusion the OLM has provided an efficient way to compare three . 
different Hydrology parameterizations. The effects of the- differences between 
these parameterizations on soil moisture and temperature, although not inter- 
active with the GCM, can be readily evaluated. The long*-term calculations 
(greater than '2 days) although not attempted at this time should be conducted 
later as soon as major revisions are completed and a data tape is available. 
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FIGURE 4-4 SOIL MOISTURE CONTENT FOR' AVERAGE CELL 



FIGURE 4-5- SOIL MOISTURE 


NOON 



FOR WET CELL 


fraction of available 



OLM-UCLA 1 (GW 4 K, RUN # 619) and OLM-UCLA 2 (GW = K, RUN #674) 


FIG. 4-6 SURFACE AIR TEMPERATURES (TS) FOR DRY CELL 








UCLA (GCM DATA TAPE)-, UCNI (#218), GISSJK (#468) 


FIG. 4-8 SURFACE AIR TEMPERATUp: (TS) FOR WET CASE 
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UCNI(#97) 
QLM - UCLAlCf/ i 
0LM-UCLA2 (#6 

GISSJK(#782) 
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UCLA (GCM DATA TAPE) , UCNI (#218) , GISSJK (#468) 


PIG. 4-11 GROUND TEMPERATURE (GT) FOR WET CASE 









FIG. 4-13 EVAPOIRANSPIRATION FOR AVERAGE CELL 








FIG. 4-14 EVAPOTRANSPIRATION FOR WET CELL 
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.033 

HH 

HHI 

.-83 

. ;83 

HHI 

;027 


GISSJK 


HH 

.08 

IHI 

HH 

.08- 

,08 

miM 

-008 


2l8 

UCNI 



.025 

m 


IIRSI 


HHHl 

.018 

' 

668 

OLM-UCLAl 




HH 

^Hl 

HH 


.01 

.027 

HHi 

0IM-UCLA2 


s 

.02 

HH 

|^H[ 

1 

1 

.01 

> 

.02 

WET 

356 . 

UCLA 


HH 

RSI 

I^H 

|R| 

1 

m 

RHH 

.025 

l»68 

GISSJK 


HH 

■1 

HH 

IHR 

.6 

HH 

.006 

.012 


NOTE: % sat/100 = (55 avail/100 ) x. (FC-PWP) + PWP 


■EABEE U - 3 

DAYTIMB AVERAGES OF PE, GW, R AMIJAET VERSUS DRY, AVERAGE AMD WET CELL COHDITIONS 















































































V. Description of the OIM Computer Program 


OLM Program 

The OLM program contains a main section wkere parallel hydrologic 
parameterizations are ■ evaluated. The first two subroutines, SDET and 
SURFWD provide needed information to the main hydrologic section that 
would normally be available in a GCM run. The third subroutine, RITE 
contains the bulk of Fortran statements needed for parameter evaluation. 

There is a total of 9 steps in the OLM as illustrated in the flow 
chart. Through step 2, Parameter Specification, the OLM receives the 
necessary data to simulate an on-line GOI run. The major ground hydrology 
parameterizations are incorporated in steps 3A, 7 and 14. Steps 7 § 14 
are similar to those in C0MP35 and discussed in Section III, while step 
3A was added to parameterize different soil characteristics. 




step IB 

Program Constants 


Step 1C 

Do Loop Parameters for Specific 
Cell Locations 


OLM00350 

0LM00700 


0LM00710 

OLM00920 

0LM00930 


0LM01220 



Parameter Specification 


Step 2 A 


0LM01240 


Variables That Are a Function of Time Step OLM01330 


Step 2B 


Step 2C 


Step 2D 


0LM01340 


Data Specification for Initial Runs 0LM01430 


0LM01440 


Initialization of Cell Variables 0LM01610 


0LM01620 


Set Sums for All Averages Equal to 0 OLM01700 


Step 2E 

Initial Data Tape Read 


0LM01710 

0LM01860 


Step 2F 

Begin Cell Loop of Model 


0LM01870 


0LM02390 


Step 2G 
Do .Loop Read 


0LM02400 


0LM02460 


Step 2H 

Solar Declination Calculation 


0LM02470 
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0LM02560 


Step 2J 

Zenith Angle Calculation 


0LM02790 


Step 3 


Surface Specifications 


Step 3A 

0LM02810 

Surface Condition and Elevation 

0LMO357O 


Separate off-line deck 
•for GISSJK .only with 
the physics described 
in C0MP35, (see COMP- 
35 flowchart^) 


GISSJK ? 


GW(J,I) = SRH + *(FC X Cl - PWP x C2)/100. ■ 

h PWP*C2 

Step 3B 

OL-M03580 

Ground Temperature 

OLM03640 

Step 3C 

0LM03650 

Surface Drag '.Coefficient 

0LM3810 


0LM03120 


HYD01240 - 1261 
HYD03000 - 374 


3 - 












Step 4 


Vertical Arrays for Some Variables 


Step 4A 

0LM0383O 

Pressure • 

0LM0396O 

Step 4B 

• OLM03970 

Temperature and Convective Adjustment 

0LM04140 

Step 4C 

0LM04150 

Moisture Variables 

0LM0423O 


Step 5 

0LM04240 

Vertical Diffusion of Heat 

0LM04460 

and Moisture (Air-Earth Interaction) 



Step 6 

0LM04470 

Relative Humidity 

0LW4570 


Separate offline deck 
for GISSJK only with 
the- physic’s described 
in C0MP35. See COMP- 
35 flowchart- 
CHYD02550-260Q) - 


GISSJK ? 


SEE 

C0MP35 
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i 


Step 7 



Determination of Surface Temperature. 

; 



step 7A 

0LM04580 


. Initialization 

0LM04660 


TIE (NLAYPl) = TS(J,I) 

0LM04620 


Step 7B 

0LM04670 


Relative Humidity Scale 

OLM04700 


Step 7C 

0LM04710 


Bulk Aerodynamic Coeff 

0LM04760 


Step 7D 

0LM04770 


Surface Diffusion Coefficient 

0LM04810 


Step 7E 

0LM04820 


Eddy Diffusivity (see;TK;and NCAR) 

0LM05080 


T&K: 

• 

EDNS = AMINKlOO., ED x EXP (.32 *DTS/WMAGNxx2) 

EORIG = EDNS 


ZLN = ZLNCO X TIE (NLAYPl) x SP/PS 
. New EDNS formulation 

N.CAR w/no countergradient- (Used in GCM Data Tape Creation) 

NCAR: . 


EDNS = 10. + 100. *(1-EXP (1200.*((TH(NLAY) - THG)/ZLN))) 

Can't use above calculation for.TH(9) >THG. Best that can 
be done offline with data from levels 7,-8 & 9. 

¥ 


▼ 
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0LM04940 


EDNS = EORIG 


EVE = 1 

0LM05040 

Step 7F 

COLM05090) 

Evaporation Scaling Factor 

0LM05360 



0LM05110 


EVE = AMIN 1 (1,(WET-PWP *C2)/{FC * Cl - PWP *C1)) 




UCLA? 


0LM05130 


EVE = AMIN 1 (1, 2 *6WC(J,D) 


Step 7G 

OLM05370 

Potential Evaporation 

0LM05410 


PREVAP < 0 and EVE < .999 


Redo Step 7F 
with EVE = 1 













0LM05420 


Step 7H 

Potential Evaporation Using LOGS Input 0LM05550 

Note: 

Potential Evaporation Calculation Differences 


Offline (Step 76) 

Wet from UCNI or UCLA 
Eve from UCNI or UCLA 
Step 7F adjust SHLE (10) 
where SHLE = f(EDV, ■....) 

SHG = QSAT (T6,PS) 
for LOGS False: 

TG from previous calculation 
for LOGS true: 

TG from 6CM tape 


UCLA ((ilCM Data Run - Step 7H 

Wet from Data Tape 

Eve from UCLA 

Step 7F adjust SHLE (10) 

where SHLE = f(EDV., . . ... ) 

SHG = QSAT(T6,.PS).- 


TG always from GCM data tape' 


step 8 

0LMO556O 

Sensible Heat and Evaporation Fluxes 

0LM05740 


Step 14 


Ground Update of Moisture and Temperature 


Step 14A 

, 0LM05760 

WI Calculation 

OLM06050 

Step 14B 

0LM06060 

Ground Temperature Update 

0LM06360 

Step 14C 

0LM06370 

Actual Evaporation 

OLM06400 

Step 14D 

0LM06410 

1 

Ground Moisture Update 

0LM06780 

Step 14E 

0LM06790 

Soil Moisture and Freezing 

OLM07220 

1 Snow Covered Land 

0LM06810 

2 Snow Free Land 

0LM06890 

3 Ice 

0LM07150 





MAIN END 


OLM07230 
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Appendix A 


Symbols arid Notations 



FORTRAN 

SYMBOL 

MEANING 

UNITS 

PROGRAM 

LOCATION 

ADINF 

antecedent soil moisture indicator 


HYD 06560 

ADDW 

water which may enter the soil 
to be evaporated 

cm/hh 

HYD 06610 
06620 

ASM 

available soil moisture 

.15 - con^josite soil 

.066 - sandy soil (desert) 

% sat 

HYD 03040 
03090 

BBFACT 

effective infiltration rate 

cm/AT=30 min 

HYD 06570 

BETA2 

scaling factor for bare soil evapotrans- 
piration rate (3) 

I 


HYD 03310 
03660 
03730 

BETA2P 

scaling factor for plant influence on 
transpiration 


HYD 03640 

BTA 

= BETA2 + BETA2P (0 < BTA < 1) 


HYD 03670 
03680 
03720 

BTA2 

argument of BETA2 coefficient 


HYD 03540 

BT2P 

argument of BETA2P coefficient for 
transpiration 


HYD 03560 
03620 

CD 

surface drag coefficient 


HYD 02140 
02150 
02160 

CZH 

(^ ) , bulk heat capacity 

10 

where X = thermal conductivity 
c = specific heat 
(0 = diumal angular frequency 

cal/ cm^- deg 

cal/ cm-s e c-de g 

cal/cm^-deg 

l/sec 

HYD 05040 
05050 
05060 
05090 
05120 
05150 
05160 
05270 


A-1 



FORTRAN 
SYMBOL _ 

MEANING 

UNITS 

PROGRAM 

LOCATION 

DENOM 

storage variable used in ground 
temp . eq . 

ly/day-cm^-day 

HYD05540 

DRAW 

surface wind drag parameter 

m/sec 

HYD02690 

02700 

04000 

04010 

DRNG 

water which is- unavailable for 
ET due to field capacity being 
exceeded loss of water 
through percolation when WETR > . 5' 


HYD06090 

06670 

DRRN 

see DRAW 

' 

HYD03490 

03500 

DSHS' 

Aq„ = q* - q saturation specific 
0 . g g 

humidity at the ground minus specific 
humidity at the ground 

■ 

HYD03990 

• DTRBET 

9/3T of (evapotranspiration) 


HYD03700 

03710 

05420 

05430 

05440 

05460 

DTSI 

difference between ground temperature 

"C 

HYD03480 

DTS 

T_ - =• AT = difference between 

G S ' g 

ground and surface air temperature 

“c 

HYD02680 

DU 

temporary storage variable 


HYD02810 

02860 

02890 

' DV 

temporary storage variable 


HYD02830 

02870 

02'910 

ECAR 

= E6DT/C1-R0CK) 

(Note: EVE in E6DT is equal to 

(l.-ROCK)-EVE so ROCK affect cancels out 
in evaporation calculation) 


■ HYD06630 
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FORTRAN 

SYMBOL 

MEANING 

UNITS 

PROGRAM 

LOCATION 

E6DT 

evaporation rate from soil at grid pt 

gCfl.2?) • /iiT=3o ■ 

cm2 

HYD05950 

ETFLX 

evaporation flux for diagnostic 
purposes 


HYD05900 

05910 

05920 

EVACO 

=EVE-DRAW 


HYD03780 

03790 

EVAL 

= 600 , latent heat coeff land 

680 " •’ - snow, frost 

= 680 + 80WI, FOR MIXWI CONDITION 

cal 

g-CHa®) 

HYD05030 

05080 

00110 

00140 

00180 

00260 

EVE 

soil surface wetness coeff 


HYD02980 
• 03740 
03770 
03830 

EXPTI 

argument of evapotranspi ration 
coeff 


HYD03530 

03550 

EXDTR 

storage variable 


HYD03690 

EXPTP 

does not permit BT2P to become too 
large 


HYD03630 

EXPTS 

keeps BTA2 from becoming too large 


HYD03650 

FILTE 

infiltration -rate derived from viscosity 
relation -fCt) only not- soil type or 
antecedent precipitation 

cm/AT=30 min 

HYD03220 

03270 

FREEZ 

amt of water that freezes or melts (%) 


HYD06050 

06770 

06810 

FLDRK 



HYD06080 

FXTRN 

not used 


HYD03020 

03450 
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FORTRAN 

SYMBOL 


MEANING 


UNITS 


PROGRAM 

LOCATION 


GAMCN 

L q* 

s 

« 1 + G K T^ 

r ^ CAZ) E ! 

1 + 54181; 

C t 2 
P s 

HYD03360 

GAMCN 

see GAMC 

HYD03360 

GAMSN 

same as GAMS 

HYD03410 

GDH ■ 

logical variable, if true CD=.5CD 

INPUT 

GNFLX 

variable for diag. purposes 

HYD05940 

GSW 

how often (xn hrs) print out will occur 


GT 

ground temperature °K 

HYD05570 

06920 

GW 

ground wetness sat.) 

=FL0ATCISNR) +.5 WETR 

HYD01210 

01230 


01260 
01270 
• 01280 
01290 
01300 
01330 
01350 
01420 
01430 
01440 
01450 
01460 
06960 

IGVET =13 

=21 -.dry land, GW = 0 

=23 • . ' INPUT 

=24 - snow thickness - 
. function of albedo 

ISNR integer storage HYD06200 

06950 
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FORTRAN 

SYMBOL 

MEANING 

UNITS 

PROGRAM 

LOCATION' 

ITC 

print variable - allows a portion 
of J § I ceils to be printed when 
GSW occurs 



lALB 

storage for I 


HYD01090 

PLANT 

a variable which decreases transpiration 
from 'plants when the solar zenith is 
greater than 75® 

• 

HYD03570 

03080 

03600 

POTALP 

precipitation 

cra/time step 

HYD07080 

PREVAP 

potential evapotranspiration 

mm-gCH^O) 

HYD03210 



cm^-hr 

03510 

QSAT 

saturation specific humidity function 

, g/g 

HYD00890 

rAdtrm 

net radiation term 

iy/day 

HYD05560 

05880 

RAIN 

' 

cm/hh 

HYD06030 

06190 

06540 

RHS • 

surface relative humidity 


HYD02530 
02550 
' 02580 

02650 
02660 

RNIT 

•=RADTRMC6.9 X 10“'^) for diagnostic 


HYD05890 

ROCK 

areal- fraction devoid of soil, 
in 5 )ervious soil 

% 

HYD03030 

ROKCJ^I) 

same. as ROCK 


INPUT 

RSURF 

decimal value for albedo 


HYD01520 

01550 

01560 

01570 

RSURFF 

variable that determines soil type 
from albedo array 


HYD03060 

3070 
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FORTRAM 

SYMBOL 

MEANING 

UNITS 

PROGRAM 

LOCATION 

RUNOFF 

water that does not infiltrate soil 
col. 

cm water 
AT=30 min 

HYD06060 

06440 

06580 

06590 

06660 

RUNOW . 

water that runs off impervious 
area and is unavailable for 
evapotranspiratibn 


HYD06100 

06650 

RICH 

Richardson No . 


HYD02930 

SAT I 

saturation specific humidity 


HYD03290 

03390 

SHG 

saturation specific htunidity at 
the ground 

. g/g 

HYDO3I90< 

SHLE 

sat. specific humidity f(DRAW) 

*' fCBETA2) 

" f(EVACO) 

g/g 

HYD03750 

03880 

03950 

SHSATS 

sat. specific humidity at the 
surface air temperature 

g/g ■ 

HYD03800 

SLBEDO 

albedO' • 

7 = ocean 
14 = land 
35 = desert 
39' = 

40 .snow 
50 = snow 
70 = snow 

.5 cm 

1.0 cm 

2.0 cm 

HYD01500 

bl490 

01510 

SMELT 

snowmelt 

cm/hh 

HYD06070 

06280 

06290 

SNFAL. 

snow fall 

cm/hh 

HYD06040 

06180 

SNFLX . 

sensible heat flux 


HYD05860 

05870 

05930 

SNOTEM- 

storage location' 


HYD01540 
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FORTRAN 

SYMBOL 

MEANING 

UNITS 

PROGRAM 

LOCATION 

SNR 

thickness of snow or change 

cm or 
cm/hh 

HYD06000 
06210 
06220 
06240 
06310 
• 06940 
06510 

SNRl 

snow thickness 

cm or 
cm/hh 

HYD06230 

TERMl 

short wave radiation 

ly day”^ 

HYD05470 

TERM 2 

sensible heat 

-2 * -1 
cal cm day 

HYD05480 

05490 

TERMS- 

evaporation flux 

If 

05500 

05510 

TBRM4 . 

sea ice heat flux from ocean 

n 

05520 

TERMS 

long wave radiation 

tf 

05530 

TGM 

=WAT, used for printout of WAT 

• 

HYD02640 

•07040 

TON 

used for printout of WAT 


HYD07030 

07050 

TGR 

ground temperature 

OR 

HYD06020 

06160 

06300 

06720. 

06820 

06900 

THG ’ 

potential temp of air 



TOTALP 

model .predicted precipitation 

m/hh 

HYD06170 

TPMN 

print variable 



UGTC 

storage variable 


HYD03240 



FORTRAN 

SYMBOL 

MEANING 

UNITS 

PROGRAM 

LOCATION 

USM 

J 

unavailable soil moisture (%) 

,025 sand soil (desert) 

.1 composite soil 

Ref. Salter § Williams '64, Table 1, p. 

•312 

HYD03050 

03100 

VALBFG 

model control variable 



VIGT 

viscosity correction factor 


HYD03250 

03260 

VGtC 

centigrade temp, of ground 

deg ■ 

HYD03230 

VPM 

composite soil moisture 


HYD03180 

VPl 

% of available soil moisture 
(temp ,stor) 


HYD031S0 

03160 

03170 

VP 2. 

% of unavailable soil moisture ( ) 


I 

VSM 

total soil moisture 

%. sat 

HYD03140 

WAT 

soil moisture available for 
evaporation scaled up 

% sat 

HYD02590 

02600 

02610 

02620 

02630 

WET 

ground wetness 

% sat 

HYD01580 

01670 

04900 

04910 

05340 

WETIT 

water available 
same as WAT 

% sat 

HYD04880 

04890 

.WETR 

amount of water 

% sat /. 

HYD05990 

06430 

06450 

06640 

06680 

06700 

06710 

06910 

WGM 

max water holding capacity or 
sul moist cont at saturation 

g/foo;t soil 

i 

HYD03000 

03110 

03130 


= 2 sandy soil 

= 2- 16,6 coB^JOsite soil g/foot- soil 



FORTRAN 

SYMBOL 

WI 

WIR 

WMAG 

IVMAGN 

WTER 

WTSTO 

XRS 


MEANING UNITS PROGRAM 

LOCATION 


% o£ ice in mixed ice and water 
1 

%WET 

HYD05190 

05210 

05220 

05230 

05240 

05250 

05350 

05730 

05740 

tf 

It 

HYD06010 

06780 

06790 

06800 

Wg wind speed at ground 

m/sec 

HYD02100 

02120 

02130 

M 

It 

HYD02710 

02720 

% of saturation for ground in 
desert regions 

% sat 

HYD01240 

01250 

stores WET 

% sat 

HYD04870 

used for point out of RHS 


HYD03200 

07070 
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Appendix B 


Code Listing - GHM § Off-Line Models (OLM) 



GHM: 

FILf: I635TFST SSSCSN A '«*« GCDDAPD SPACE FLIGHT CENTER *** 


bd 

I 

H 



SUQPCUT INE CCMP3S(J» JCFNT-, JUP.KCN^KCENT ,«LP) 

C 

C . THIS SUERCUTINE WAS ADDED TO PREVENT CCNP3 FRCV EXCEEDING THE 
C COMPILER. LIMIT FCR OPTIMIZATION (CPT=?). IT CCMFRISES THE INNER 

C LCCP CF CCMP3» 

C 

C***#: THIS IS CCMMCN ELCCK CMS! CREATED MAY £ 74 
LCG1CAL*4 FLAGS, SFINUP.GDH 
REAL+4 KAPAtLAT 

DIMENSION’ LAT(4e),DXU<46),OXP{4e),CYLM€),OYP(46),SINL{4 6>, 

1 COSL (4F),0XYP(46),F{45),C(3CC),CI<3rC),r.LMNV(7?) 

C IMONSICN PH IS I 4 £, 7?) ,MAP (4 6,2) ,F2(4e,7S ) 

CIMEMStCN 0(72,S,11,4> 

DIMENSION CT (7? ,9,1 1,4) 

DIMENSION U{72,c,n),v(72,9,ll>,TI??,c,n),SH(72,C,n) 

D IMENSI ON UT ( 72,9,1 1 ) ,yT{ 72,9, 1 1 ) , TT ( 72 ,9 1) ,SFT ( 72 ,9, 1 U 

DIMENSION P( 72 , I I ) ,PT( 72, 11 ) 

DIMENSION PU(7r,c,.->),pv(72,9,2),'=VLP(72.S.2),PUCL0(72,0,2) 

DIME NS I ON C0NV(72,9),SU<72,9),PIT(72) ,SC1<9) ,CCNSP(9) 

DIMFNSICN NMSAVEI46) , ALPS AVI 46) 

DIMENSION PH I ( 72 ,P, 2), PHI OLD! 72, 9, 2 ),SPA(72,9,?),SFACLD(72,9,2) 
DIMENSION PU2( 72.9,2) ,PU2CLO('T2,9,2) 

DIMENSION PU3(72,9) 

DIMi^NSION US (72 ,2),VSI72,2) ,05(72,2,2) ,\JINCSG(72 ,2 ) 

DIMENSION UF IL (72 ,9 ) , VF IL (72,9 ] , TF IL ( 72 ,9 ) ,SHF I L ( 72 ,9) , PF IL (72 ) 
DIMFNSICN GRLP( 46,72 ,4 ) 

D [WENS I ON SHS( 46 ,72) , GT (4 6, 72 ) ,GW ( ,7 2 ) ,TCFCC (46 ,72 > , T S (46 , 72 ) 

COMMON, 

A JSP»JNP,IM,NLAY ,PTROP ,ISTART,JSPP1«J'NFM1,FIM,.NLAYM1 .NLAYPl , 

4 J1 , JM ,KM,T AUT , IROT ,M.ROT,TALTHT ,MRCTHT , 

* NR, JAYSI 12 ), 1NCS( 1 1 ), JSE, JNO,DLAT ,DL C N . J TE ST , HE ST ♦ 

* DT,TAU, ITAL ,X1NT ♦ 1 DAY.JCAY.TCFDAY , JTATE , JMONTH (2 ),JYEAR.NSTEP, 

* N CY CL E, NC0NP7, NHOG AN ,TAUP, TALI, TALE, TALC, 

* P1,GRAV,RGAS,KAFA,PSL,FD,FML-,NFLV',R3F,MRCH,RSriST,£IND,CCS0, 

A NC 3T, M R rC 3T,LC3T, PR CJAR, TOT AES, SO. ISPACE, SOM ULT .CUMMYI (5) , 

* 1 VER , CDh ,RHMAX ,FCU,SP INLP.MACh IN , 

A I ALTER , NNI M E , F T E , PT E2 , Nl IMPT S , 5ETMM,MINJ,MAXJ,IRANC,IRTYPF,ISS, 
A SS , ID lAG.PEP ICD ,ALPHA(4) .BETA, GAMMA ( 4 ) , T I ME 0 ( 4 ) , NT FB , R I NFLU , 

A MODEL T ,LR AC ( 16 ) ,NMCT, NRADC )6).KTFA,JTAFC,COX.FLAC-S(10), 

* MINI, MAXI , CLH, CLMM.Y2{20) .TINT, 

* XLABEL(20) , S I G ( 20 ) . DS I G ( 20 ) , ? I G E ( 2 1) , C £ I GO { 1 9 ) , 

A JlPS(ll).JMPS(ll),JiUS(n),JKLS(ll)»JIF , JMP , J I L , JML , J 1 PV, JMPV, 

A INC . J IPY, J 1PM , JMPY , JMPPl , JMFP IX , J U > .JlLMl.JltMlX, JMUX, F INC , 

A INC2, INCH, 

* LAT , DXU ,DXP, CYU .DYP ,S INL.COSL ,DXYF,F,CLMMY,NTAE.FHJS,MAF,F2, 

A P, 0 ,PU, FV, FUCLD'.PVUP, C0NV,P1T,FHI,PHICLC,SPA,SFACLC,PU2, 

* PU2CLC , FU3 ,P1TSF,SCI,NMSAVE, ALFSAV 
EOUIVALENCE (J£F,C(1)) 

EOUTVALENCE ( P ( 3 , 1 ) , C U 1 ) ) 

EOUIVALENCE ( P ( 1 , 1 ) , PT ( 1 , 1 ) ) 

EQUIVALENCE ( 0 T ( 1 ♦ 1 , 1 , I ) , C ( 1 , I , 1 , 1 ) ) 

EOUI VAL ENCE (0( 1 , 1 , 1 , 1 ) ,U ( I , 1 , 1 ) ,L T ( 1 , 1 , 1) ) 

EQUIVALENCE t Q ( 1 , J , 1 , 2 ) , V (1 , 1 , 1 ) , V T ( 1 , 1 , I ) ) 

ECU I VALENCE (0(1,1, l,2),T{l,l,n,TT(J,l,l>)- 
FOUI VALTNCE (C ( ) , I , 1 ♦ 4 ) ,SH( 1 , 1 , j ) , SHT ( 1 , 1 . 1 ) ) 


HYDCOOIO 
HYDCQC20 
HYDC0030 
HYD00040 
HYD00050 
HY000060 
HYD00070 
HYDOCOBO 
HYD0C09O 
HYDCC too 
HYDCOn 0 
HYDOC t20 
HYDOO 130 
HYD00140 
HYDOO 150 
HY000160 
HYDOO 170 
HYD'^0 IBO 
HYDCC190 
HYD0020C 
HYD0021 0 
HYDCC220 
HYO0C230 
HYC0C240 
HYUCC200 
h!YD00?60 
HYDC^pyO 
HYD002B0 
HY000290 
HYD003C0 
HY000310 
HYD0C32O 
HYU0033C 
HYD00340 
HY000350 
HYDOO 360 
HYDOO 370 
HYOC03B0 
HYD00390 
HYD00400 
HYD0041 0 
HYD00420 
HYDO 04-30 
HYD00440 
HYD00450 
HY000460 
HYD0C470, 
HYD00480 
HYDOCARO 
HYOOObOO 
HYOrOS! 0 
hYDOOSc’O 
HYOoos'm 

HYDO 0 54 0 
YDOr 


ON-LINE LISTING CGHM) 



no no 



w 

I 

tsJ 


FQUIV/iLENCe {SCI It t >,CONV(l, n ) 

E;QU IV ALENCE (SC HU, CGNSP£ J ) J 

EQUIVAL ENCS <QS ( 1 ♦ 1 » 1 ) »US,< i,l),P*U()ti,U>tCVS(l.l>,C£<l,l, 2 >, 
I FU ( 1 t 3 . 1 ) > 

S cur VALE N''cfe U'INCSQC i;i > ,PU( 1 . 5,1 n 

ECU IV AL ENCS £UF IL < I , 1 ) , U( X » 1 , 1 ) > , < VF IL ( 1 , U , \ < 1 , I , 1 ) ) , 

1 (TFtL(l»l),T(l,l,l))t( SHF 1 L( I t U X SF ( I . I . 1 ) ) t ( PF I L (1 ) , P U , 1 ) ) 
ECUtV ALENCE < 6 FUP ( 1 , 1 , 1 > ,EHS( I , 1 ) ) , ( GF CP < I , 1 . 2 ) ,G 1 < 1 , 1 ) ) , 

1 ( GPUP( 1, 1 , 3 ) . G\\( 1 , 1 n . (6RUP( 1 t 1 , A ) . TCFCG ( 1 , n ) 

EOUIV AL ENCe (GPUP( l,l,l>tPU{l, 7 ^ji),(I$(Ul),U{l, 3 , 6 )J 
C*##* ARRAYS FOP SPLIT GRID 

l.CG ICAL J IPL ( AC ) , J lUL (AO . JFPL ( 4 gUjFLL(Ae) 
iNTEGER ■I^CP (A-' )» INCU(4e> 


HYD00560 

HY000570 

HyOOOSSO 

HY000590 

HYD00600 

HY000610 

HYOC0620 

HYD00630 

HYD006A0 

HYDOOeSO 

HYD00660 

HYD00670 

HYD006ftO 


CCMMON/SPLIT/ J 1 PL , J 1 UL- , J I^PL. , J|v!LL i t ^ CP , I ^CU HVD0C6f)0 

C . . HYD00700 

I.CGICAL LAND «OC€ AN, I CE» SNCW.MtXtt t »FFC£T .FREC HYD00710 

COMMON/ FACCON/ AS ( 1 5 ) . RE ( I C ) , PL ( 1 5 > . F LE ( 1 C ') ,PL K ( 1 E > , PL KC ( 1 6 ) , T L ( 1 5 ) hYD00 720 

* ,TLE( 1 6 } ,TG ,TH( 1 E) ,SHL ( 1 5 ) , SI- L^ ( I 6 ) .SE C .CLOUD ( 1 E ) ,C C S/ , S A , SG , C XL , HYOC C 730 

* OCEAN, ICE.SNCV ,SHSAT (15>,GAM{1?. )tFMlE),£S5(l£),SE£F(15),Hh(15), HY00C7A0 

* HHF( 16 ) ,HHS ( 15 ) ,CVT( 15) ,CVQ( IE) ,CXDE{ IE ) ,PREC( IE ) ,EXL< 15) , HYOOC75C 


* CGE{ 15 5 , HYD0C760 

Ci'TiViVON/ C.VFoCM/TftCFI , CTC3 . 50AY ,2 LnCC,TICE»FM,CTI,CTIC,>-1CE, HYDCC770 

* GAMFAC ,C0SF ,CCFFS,RSCIN,ROT,Plleb,ENCViN,£NChS,CCtFi ,C0EF2 , HYOOOVaO 

* WMAGSP ,WMA<jNP,Ml .yE.NSsNCL.LCL J ,LCL2»£1EG«CP HYDf^C790 

CCMMCN/GRCHYC/ ICVET , ITC t GSW.CCN 1 tCCN2 ,CC N3.CCN4 ,CCN5 iCGN6 ,CCN7 HYD^OFOO 


CCMMON /ALPCCM/ R SURF » SLB EDO ( A C ♦ 7 ? ) , V AL E F G , J AL6 , I A L E 


HYOOO SIO 


CCMMON /SOOTAE/ S COT ( 46 , 7^ , 8 > . MN M X T S 
DIMSNSICN RON( 46,72) 

CCMMCN/FOKK/FOk 
1NTF.GER*2 SLEECC 
LOGICAL V ALB FG 
LOGICAL 44 GWSWCI- 
C44-** LrCICAL FFAJ 


QSAT(X,Y)=,622*( 1 C , * * ( 9 * 4 C51 -2 3 C: 3 . > X ) )XY 


HYDOC 820 

HYDC0830 

HYD00840 

HYDC0650 

HYoooaao 

HYDocero 

HYDCC880 

HYD00890 


SWSWCF = .FALSE, 

GWSWCH = AMCOi ( TAU-TAUI+GSV»)/G! 5 V., ) . ) .LI .tCCOl 
IF( J.GE ,34 )GVtS\S CP = . FALSE. 

IF(J.LE .le)G\^SWCH=, FALSE, 

IF (GWSV»a-,ANOt( J-l + ITC)/lTC*nC»EC. J4 nC -1 ) V.P ITF (3 ,2555) 
IMINC=IM 
OC ang 1 = 1 , IM 

SET THESE GLANTITIES TC ZEPr-feQR USE IN FPEPAFING COMct taPE 
WITH SPLIT GRID 
T ( I t 1 , M E > = C . 

T{I, 2 ,M 2 ) = 0 , 

SET THESE OLANTITIES TC ZCRC FCR LSE IN PREPARING C0MP3 TAPE 
surface FFICTICN TERMS 
P ( I ,M2 ) = U( I , 1 ,M3 > 

' Pn,M3) = U(It2,N3) 

U( I , I ,M3 ) - C, 

80S U(I,2,M3) = C, 

JALB=J 

00 ei 0' i=jNc , iM , snc 

1AL0=1 

HACC5 = CC3r'*CCS(R0T + ( I-l )^DLON) 


HYdCCSOO 
HYDCCOl 0 
HYDC0920 
HYD00930 
HY000940 
HYDOCSSO 
HYDOC060 
HYD00970 
HYD00P80 
HYOCC9O0 
HYDOIOOO 
HYDO 1 Oi 0 
HYDC 1020 
HYDO 1030 
HYDO 1040 
HYDO 1050 
HYDO 1 060 
HYDC 1070 
HYDO 1080 
HYDO 1090 
HYDO 1 100 



F ILF* 


ie35T£ST 5SSCSN 


A 


GCDDAPC SPACE FLIGHT CENTER ' 


w 

I 




AO ) 


C0S2“SI HL ( J ) ^S I^C-^CCSL( Jj^^HACOS 
SURFACE CEND IT UN 
OCEAN = T C PEG ( J . I ) *GT * I ^ 

ICF=:TOPCC ( J » 1 ) .1 T *-9*PE*S 
LAND= .NET ^ ( I CF. ,CP .OCEAN ) 

C + FfrAJ= .NCT .DC FAN* AND# ( IGVET'.EQ# 1? ) 

ir (GW (4e ♦ VD “1 C . ) 2SG4 ; 2 55 4 I 2E5 3 

P:b5 2 SNOW = L Ar.O. ANC.CW( J» I ) .GT. 1 #0C 
TF{ IGVET .EQ. 24 ) CC TO ?55? 

<=;N0W •= LAM:#ANC^( CW( Jf 1 ) .CT.l »C.CP.SLeFCC(JfX) .CE 
IF( SNOW ) GW( J t n - AMAXK10.5 $ GW(J.n) 

GC TO 2552 

2554 IF(GW( J.I ).LT,(-0.n')GW( Jtl ) - 1 .» 0 

WTER=GW ( J • I )^5 .CQ-f. 2 690 91 795 6 
WTER=.54WTER 

IF(LAND } GW( J» r )-,5i^(GW( J. n^.3+^4)‘ 

TF( SLBECC( J . I ) .EG. 2 5 ) G W ( J . I) - A H I M ( 1 . , WT E F 

Gw’( J » n “ CW i J , I )>/2 . 

rF(LAND»AND. IGVFT.EO.Pl)GW(Jfl) = C.C 
ir(CCEAN) GW(J»n = 0.50 
IF( IGVET .FQ. 24 ) GC TG 2552 1 
rF(OCEAM GO TC 2552 
IF (LAT ( J ) .GE .SNCUN )GW( J t I)==GW ( J , I l-MM ( ( L A T C J ) -S NC WN ) / ( L A T ( J NF ) YDO 1330 
* SNOWN) ^3^E4) HYD01340 

IF( LAT( J ) .LE .SNOWS j GW ( J t I) GW(Jtl) I N T ( { L A 1 ( J ) - S NO W S ) / { L A T ( 1 ) HYDO 1 350 


) 


HYDO 1 no 
HYDO 1 120 
HYDO 1 1 30 
HYDO 1 140 
HYDO I 150 
HYDOl 160 
HYDO I 170 
HYDO 1 130 
HYDOl 190 
hYOOl200 
HYDO 1210 
HYDO 1220 
HYDO 1230 
HYDO 1240 
HYDC 1250 
HYDO 1260 
HYDO 1270 
HYDO 1230 
HYDOl 290 
HYDO 1300 
HYDO 1310 
HY001320 


2S5?l 

^ sir 

C « 4 iV* 


-SNQViS) *3*64 
GC TO 

IF( SLHEDC< J • I ) 
SNOVJ THICKNESS 
PESIDES F4VINC 


) 

.LP,39) GO TO 2552 

IN CMS. = 0.5 » X .0 »2.,C FCF 4 L EEC C=4 0 t 5 C , , 70 PERCENT* 

ANC FERNNM ICE* 


GivI J , I >. = 2* 
iriSLBECCI J* 
IF(SLBEDC< Jt 
IF{ ICE) CW ( J» 


2555 


2556 


25566 




THICKNESS CF 30OS 47 ECilFFCLE 
'AITF MAXM ALDEDO CF 8C PERCENT* 
c c r E C 3 +0 • 5 

n'.LTlsO) GWCJ.I) = C.5C0EC34 0.5 
1).LT,70> GW(J.I) = 1.0F02 + C*5 
I ) = 50.0F03+C. 5 
IF{J,Le.9*ANC.L/NC) GW(J, I)=3 0.CEC2 + C*f 
SNOW = LAND. 4NC,CVk{ J* I ) *GT. 1 *0 1 
TF(VALBFG) GC 7C 2556 
SLBEDOI J. I ) = 14 

IFCDCEAN)SLBenC< J.l ) = 7 
IF{ SNCW oGR. 1 CC ) SLaUDCIJ.I) = 7C 
PSURF = 0 .01 =tSLEECO( J . I ) 

IF{ IGVET *NE. 24 ,CR.GW( J* I ) *LT. 1 *C 1) GC TC 55566 
SNCTEM = SORT ( (GW (J » I )-0.S )/lCCC . ) 

RSURF = 0 .50 *ShCTEM 

IFISNOTEM.LT 0 1 . )RSURF = R SURF + 0 * C 1 * SL B E CC ( J , I ) =j« U . -S NCTE N )‘ 
RSURF = AMN 1 (RSURF .0 .80) 

WET = AFI N) < I, t ( A^OD(GW( J .,I ). 1 . >*2 ) ) 

LAND=LAND »ANC , .NOT. SNOW 

MIXWI = GT(J.I> .LE. 1.1 . AND. (LANC.CR.SNCVn) 

FPOST=L ANC. AND «(GT(U»I).LE.TICE).ANCo.NCT.)/I>WI 
IF( .NOT .OCEAN ) Z =FH I S ( J . I } /GR A V 
GROUND TEf/PEFATURF AND WETNESS 
IF(OCEAN) TG=TCPCG(J.I) 
ir( ,NCT .CCEAN) TG=GT(J,I) 


HYDO 1 360 
HYDO 1370 
HYDO 1330 
HYDO 1390 
HYDO 1400 
HYD0I41 0 
HYDO 1420 
HYDO 1 430 
HYDO 1440 
HYDO 1450 
HYDOl 460 
HYDO 1470 
HYDO 1430 
HYDO 1490 
HYDO 1500 
HYD01510 
HYDO 1520 
HYDO 1530 
HYDO 1 540 
HYD01550 
HYDO 1560 
HYDO 1570 
HYDO 1 530 
HYDO 1590 
HYDO 1600 
HYDO 1610 
HYD01620 
HYDO 1630 
HYDO 1640 
HYDO 1650 
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^ 2S 

G) 
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S3 
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F ILF- 


IfUSTEST l-SSCSN 


A . *** OCDOAFn SFACfc FL16HT CFNTEF 


M 

1 

-P* 


IF(M IXW I ). TG = T ]CF 
IF{ *N'CT-»L ANC ) V.ET=IoO 
C 
C 

PFESSUW E5 

SF = P( r.KGENT ) 

CCLM9 = PN'/£P 

DC -30 L=1-»M.AY 

FL(L )=SJf'(L)>‘SF + PTNOr' 

PLC. (L )=EICE(l. )*5P-«PTPrP 
^0 PLK(L )=EXPIJYK(FL<L) ) 

PlE{NLAY4!) = SP-t'n PUP 

CAtfA# TEf^PERATUHFS AND AUIADATIC ADJUS1WENT 
DP 35 L= 1 ,KL AY 
3f« TL(U) = T ( I ♦!. f KCEM » ' 

DC AO N=1 ,3 
DC 4'^ L=2iNLAY 
LK/1=L~1 

IF <TL (L-1 )/PLK (L- 1) .G£. TL (D/PLK (L ) ) CC IC 4C 

TXETA = (TL {LN< 1 ) «DS IG( LWl )+TL<l- ) x-DSIG^L 2 )/ (PLK {L^' 1 )4CSIG( LM ) + 

« , PLK(L I’f^rSICIL ) ) . 

TL(L. - 1 ) =TPET A*PLK<L-1 ) 

Tl-{ L) =T^ETA»^FLK(L ) 

40 TH(L- n*:TL{L-l )/PLK(L-l ) 

TH(NLAY )=TL ( NL A Y ) /PL K I NLA Y ) 

C SAVE CHANGE I TEMPERATURE CAUSED PY ACIAEATIC ADJUSTMENT HEATING 
DO 4 £ L = ! «NL AY 

45 U(I»L,M) = TL ID -T( !♦ L ♦KCENT ) 

MCISTURF VARIAEIES 
CC 50 L= 1 tNLAY 
SHL ( L ) = CHI f , L tKCENT) 

50 CLOUDIL ) = 0 , 

CLOUD lNLAY-1 ) = f 
CLOUD (NLAV+? )=r« 

CLOUD (NLAY + 3 )=r . 

C**** ELIMINATE NEGATIVE SPECIFIC HUMIDITY 
C.C 202 L = 2,NLAY 
LM]=L-I * 

IF(SHL( LM U «CF.C« ) GO TO 20? 

SHLIL 2=£HL IL ) + SHL (LMl )*CS IC(LM1 }/CS JG <L )' 

5HL ( LMI ) = r , 

202 CONTINUE 

IF(SHL(NLAY) ,LT,0.) SHL(NLAY) = C. 

drag COEFICIENT 

IF ( J«GTo JSP »ANC*J.LT* JNP > IvMAG = .SV £ CRT { WI NDSG U . JCEN T ) + 

1 WINOSQ { IM INC t JCENT )4WIMDEC( I ♦ JLP) 4V I NC S C t IM I NC , JL F ) > 

IF( J,.eQ.JSP) WMAG=WMACSP 
IF(J*EG.JKPJ WMAG=WMAGNP 

IF(OCEAN) C0=AMIN1{( l»4*0 7*WMAG)T'*CCl«.rC25) 

IF(,NCT»CCr.AN> CC = oOC2 + ,Ct:e*2/ECrC» 

IF( GDh ) CD=C .5’^GD 

‘veRTI CAL CIFF.LSICN CF HEAT AND MCISTLFE 
70 LMIN = 2 

C CHANGE THIS FROM »LMIN = NLAYMl* 

DC 75 L = 1 , NL AY 



HYDO 1660 
HYDO 1 670 
HYDO 16S0 
HYDO 1690 
HYDO 1 700 
HYDO 1710 
HYDO 1720 
HYDO 1730 
HYDO 174 0 
HYOC 1 750 
HYDC 176 0 
HYDO 1'770 
HYDC 17f»0 
HYDC 1 790 
HYDC 18C0 
HYDC 1 UlO 
HYDO 1«20 
HYDO 1830 
HYDO 1840 
HYDOl £50 
HYDC I860 
HYDC 1870 
HYDO 1880 
HYD0I890 
HYD0190O 
HYD01910 
HYDO 1920 
HYDO 1930 
HYDO 1940 
HYDO 1950 
HYDOl 960 
HYOO 1970 
HYDO 1980 
HYDO 1090 
HYD020CO 
HY002010 
HYD02020 
HYD02030 
HYD02040 
HYC02050 
HYDO 2060 
HYD02C70 
HYDO 2 060 
HYD0S09.0 
HY002100 
HYD02110 
HYD02120 
HYDO.? 130 
HYDO 2 140 
HYD02 150 
HYD02 160 
HY002170 
HYD02180 
HY002 190 
HYD02200 



F U F 


{fj^STCST 


S S C S ^ 


A 


*** GGCCAPC eFACG fLIGt-T CEMEF >'** 


W 

I 

Ol 


C IF TFC VEPT rCAl riFFUSICN FE/^TING TER^S TC EE SAVED CK 
C ARil NTT CCMFLT[:0 IN TFE CODE THEY I LL OFF/LLT TC 2EFC 
7 ? U(I,L»W?)'r: To 

DC SO L=L^iN.NLAY 

CTFTA = ( ThO:i^ 1 )“TH (L ) ) v(PLML^/l )■» FL K ( L ))=«■. £ • 
D7UP=SP G {L-N^H >^RGASX'TL ( L V 1 ) / (Ft (L Y 1 > AV > 
D7DN=SP> *CS1G (L)*FGAS*TL(L )/(PL (L)>fGRAV ) 


TFE CDP/F3 TAF-C HYD02210 

HYD0,?r2C 
■ HYDD2230 
HYD02240 
HYD02250 
HYD022&0 
• H YD 02 270 
•HY0022H0 


CDLE= o2^FC 

TEMr> = DTC3>' (DSIC(LM)+OSIG(L))/(C2tF40?rN)>»*2 

FLF=-2o *'EEL EYC7CTA*1EMh 
TL(LN'l )=TL(LM }4rLE/CSIG( LMl) 

TL(t )=TL(L J-FLE/CSIG(t ) 

U(TtLWi,N‘2) = t'( I .LMl , > + FLE/OSlG(LM J 

C SAVE VERTICAL OIFFUSICN HEATING TEP.VS 
U(r,L*f^2) = -FLE/DSIGCL) 

TH( Lf/ 1 ) =TL ( L M )/FLK (,L Y I ) 

TH{ L-)=TL ( L ) /FL K (L ) 

DSH=SHL (LP 1 ) “SHL n. ) 

FLF=-?.=!'ErLE=<'DSFA T£Mp 
SHL (LW 1 ) = SHL (LM ) ■♦ELE/DSI C( LMl ) 
eO S»-'L(L) = £HL('L )-CLE/DSiG<L) 

SATUPATICN SPECIFIC HUMIDITY 

DC 100 t= 1 , NLA Y 

SHSAT (L )=CSAT(TL(i > ,PL(L) ) 

GAMCLlsGAMFAC^SHSATfLl/TLfL)**^ ’ 

FH(L)=SHL(L)/SHSAT(L ) 

ICO CONTINUE 

C*v.** DETFfVI hATI ON CF SURFACE TEMPERATURE 
FS = FLE { NL AY+ 13 
PSK^E't.PEYKC PS 3 
TLEfNLAYP 1 ) = TS { J . I ) 

RHS = 0 » 

THG = T£( J, I 3/FSK 

r F( WFT+RH ( NL AY ) « NFo Co ) PH S = 2 o •« E T Y RH ( M. A Y > / < IaET + PH ( NL A Y ) ) • 

C RH SCALE 

I F( oNCT ,LAN0 3CC TC 101 
RHS = P o 

WAT=2 .* ( I ,fre f6eC'Ce*2 .*'WET-.eee e ceccE 3 

IF ( SL0ECC ( J f I ) <»ECo35 ) h A T= » 2 =* < 2 * =} ViE T- <> 2 6-3 C S X 7 S 5 E 3, 

IF( WAT, CT , 1 , )W AT= 3 . 

IF(WAT,LT«0‘)WA7 = r, 

\vAT = WAT*( 1 ,-ROK ( J , I ) 3 
TGM=VvAT 

IF( WAT+RH (NL AY), G7o0 ) RH S= Vi A 7 *RH ( NL A V ) ✓ ( W AT + PH ( KLAY ) ) 
RHS=2»*RHS 
101 CONTINUE 

OTS=TG-TLn(MLAYFl ) 

I F< DTSo CE ,0 ) CPAW=CCF{ WM AG 4SCP 7 ( D 7 5 ) ) 

IF (OTS.LToOo) CR A-A = CD*WMAG*( Wf'AG’f'WM AC ) / ( ( WN AG* VtR A G 3 -7 . *D T £ ) 
«MAGN=1 0 

IF(0T£ ,LT ,o . ) \^^'AG^=»2* (WMAG+ ,0 1) 

7LN=ZLNCC*TLE( NLAYP 1 ) ♦SPyPS 
IF ( TH( ^LAY) oGT oTHG 3 GO 7C 110 

ECNS = 10 V+ 1 CO .A ( 1 o-FXP ( 12CCo* ( (TH{ ^L AY )-THG )/2LN3 ) ) 


HYD0 2 2<)0 
HYD023CO 
MYDC2310 
HYD02320 
HYDC2330 
HYD02340 
HYD02350 
HYD02360 
HYD02370 
HYD02380 
HYD02390 
hYD02400 
HYD0241 0 
HYOC2420 
HYD02430 
HYD02440 
HYD02450 
HYD02460 
HYD02470 
HYDC2480 
HYD02490 
HYD02500 
HY002510 
HYD02520 
HYD02E70 
HY0 0 2f>4 0 
HYD025S0 
HYDC2560 
HYD02&70 
HYD02580 
HYD02590 
HYO026C0 
HYD0261 0 
HYD02620 
HY002630 
HYD02640 
HYD02650 
HYD02660 
HY0C2670 
HYDOP660 
HYD02G90 
HYDO27C0 
HYD0271 0 
HYD02720 
HYD02730 
HYD02740 
HYD02750 



no 
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lei^STEST SSSCSN 


A 


'K** GCDDAPC SPACE FLIGHT CEMEP *** 


W 

I 

o\ 


GC TO 1 1? - • HYD02760 

CC^^PUTAT1^N CF DIFFOSION CCEFFtClEhT EEHsEEN SURFACE AND LAYER HYD02770 
•NLAY» WITH NCAR FCP I^LLAT 1 ON GUT NC CCL H TE F C H AC I EN T HYDO27G0 

no IF (J,eC<»JSP) GC to 113_ . . HYD02?90 

IF (J»EC.JNPJ .Gf-TO 1 l'4‘ " " HYD02600 

DU = ,2E*(U{ ♦NLAY»KUP)~US( IMINC ,JLF )+U( UNLAY ,»<UF)-USU .JUP) + HYD02810 

1 U( IMrNC.NLAV,KCENT>-US< I M NC » JCENT ) 4 1 < UNL AY » KC EN T ) -U S ( 1 tJCENT) ) HYD02820 
DV = .25*<VC1NINC tNL AY ♦ KUP ) - VS I I V I NC » JL F ) + V ( U NL AY , KLP ) - VS ( I , JUP )+HYD0 2S30 
1 VC I VINC,KLAY,KCENT)-VSC I M NC ♦ JCENT ) 4V II * N LA Y ♦ K CE N T ) - V S ( I , JCENT ) ) HYD0 2840 
GC TO I 18 

113 OU = »54 (U( I M NC t hLAY«KUP }-USC IMINC » JLF > tU ( I .NL AY , 1<LF )-US ( U JUF) > 

OV = ,5*C VC IMNC.NLAY.KUF)-VSC IMNC,JLF)4VCI,NLAV,KLP)-VSCUJUF)) 

GO TO 1 18 

llA DU = ^5*( U( 1 MNC» NLAYtKCENT >-LS( IMNC .JCENT) +UC 1 ,NL AY.KCENT)- 
1 use T .JCENT ) ) 

DV = .5*(VC1MNC.NLAY.KCENT)-VSIIMNC,JCEM>4VCUNLAY.KCENT>- 
1 VSe I , JCENT) ) 

IIP PIC'-< = C-RAV« C(THCNLAY)-THG)/'ZLN)y'nC* C <DL/ZLN)«*2+ {DV/'ZLN)4*2t 
1 1 .E- 1 ? ) ) 

EDNS = IT ./{ lo4AC«*RICH)42* 

112 CONTINUE 

ECV=EDNS/ZLN 

fVE= AMIM C 1 • ,? 4*WET ) 

**#*4*4, N£;W BETA FGR GF • FVAFOPATICN ******** 

WGM=2 .* 16»C6F6f 

IFCSLPEDCC J. [ ) .LT.IO) GO TO 1C2 
FXTPN= 1 
RCCK = ROK ( J, I ) 

ASM= ,15 
USM = . 10 

RSURFF=( SLHECC ( J . 1 )+.01-.35) 

RSUPFF= ABSC RSUFFF 

IFC A3SC RsSURFF ) «GT .-,01 ) GC TO EAIS 
ASN=T ,065944 
US:.<='T »0 2A 278 
WGN' = 2. 

6415 CCNTINUE 

WGM=WGM4S. 

VSM=24WET*< ASM4USN) 

VPljrAMAXl <0 . .( 2.. *ViET*(ASM4,USM)~LSV )/ASN > 

,2.4VIET*CASM + LSM)/L£V) 

,VF2) 


10 ? 


VP?=AMI N 1 { U 
VF2= AMAX 1 <0 
VPM= , 14VP2+. 2>*VF 1 
SHG=OSATC TG, FS ) . 

XRS = 0 . 

PFEV AP=:C , 

FILTF = 1 

VCTC=GT < J . I ) -T ICE 

UGTC= 1, 76 25*BXPC- .036AVGTC ) + . A45 62 
VIGT=U/LCTC 

1 FCVGTC ,LT ,n ) VIGT= UOE-21 
F ILTP = F ILTF4VI GT 

SHLE(NLAYPn=C DFAW*SHG + EDV’»:£HL (NLAY ) )? CCRAW + EDV ) 
SAT I = QSAT CTLE ( NLA VP 1 ) ,PS) ■ 

CCNTINUE 


HYD02850 

HYDC2860. 

HYD02870 

HYD02880 

HYD02890 

HYD02900 

HYD02910 

HYD 02920 

HYD02930 

HYOO29A0 

HYU02950 

HYD02960 

HYD02970 

HYD02980 

HYD02‘390 

HYD03000 

HYD0301 0 

HYD03020 

HYD03030 

HYDC3040 

HYD03050' 

HYD03060 

HYD03070 

HY003080 

HYD03090 

HYD03100 

HYDO3U0 

HYD03 120 

HY003130 

HYD03140 

HYD03 150 

HYD03 160 

HYD03 170 

HYD03130 

HYDC3 190 

HYD03200 

HYOC3210 

HYD03220 

HYD03230 

HYD03240 

HYD03250 

HYD 03 260 

HYC03270 

HYD03280 

HYDC3290 

HYD03300- 



FILE leaSTEST SSSCSN 


A 


GOOCARC SPACE FLIGHT CENTER 


FETA2=1 _ 

IF( LAND .AND. *NCT .SNCVOGO TO 123A 
GC TC 1235 
1224 CONTINUE 

TEMP=SATI/TLE<NLAYP1 > ■■ 

GAMCN=7LN*RHS + 9 . giK 1 1 . 4CLH^TEHF7K AP A ) / ( I , 4GA |Vf AC * TEM F/ 

ITLE(NLAYPl)) )*KAPA/RGAS 

T NEH= { OP A\v*T G+ECV*( TH INLAY J*PSK-G A HCN) (CRAU+ECV ) 

SATl=C«AT(TNEh tFS ) 

IRISHLE (NLAYFl ),LE.SAT1 » CO TO ei2C 

GAMSN = G AMFAC >ISAT I/TNEH*«2 

TEMP= ( SHLE( NLAYF 1 ) -S AT I ) / U • +• G A (x SN ) 

TNEH=TNEH4TEVP*CLF • 

ShLEC NL AYFl JsSFLE (NL AYPl J-TEMP 
FXTRN=1 ./EET A2 
6120 CCNTINUE 

RV-OS=PS/'( FGAS*TNEh) 

DTS -.5* { CTS+TC-TNEH ) 

IF( OTS I .GF .0 «) CRRN=CCt (WNA6 + SCRT( Cl SI)) 

IFIDTSI <LT.n,)nRPN =:CD*WJ^AG* (W^A6*V»VAG)/{ (HMAG4 ViNAG )-7**DlSI) 
PREV AP=CRPN*RHCS* 10 6 . ♦ ( QS A T ( TG , F S ) -SHL E ( N L A Y F I ) ) >«2ecr. 
IF(PREVAF.LE .0 . ) CO TO 1225 
EXPT I=AN! IM ( 17C at ICt 7S*VSM) 

ETA2=»0 l*VSHAEXF(EXPTn/PREVAP 
EXPrr = Ar/INt( 17C..I6.7 5-«:VPM) 

?T2P=,10*VPH4EXF<EXPri)/PFEVAP 
PLANT = AMM( I , , C CSZ / « 25 E £ ) 

I F( COSZ .LT tO ) FLANT = Oa . 
tF( PSURFFaGT *.C 1 ) GO TO 4C0C 
FL ANT= 1 . 

4000 CONTINUE 

BT2P-=FL ANT-«!ET2P 
EXptP=-A.VIK1 { l?C . .QT2P.) 

PETA2P= 1 t-EXFI EXPTP) 

CXPTS = “AYTN1 (17C »,,3TA2> 

PTTA?-! .-EXF<eXFT£) 

BTA=0rTA24BETA2P'' ' 

IF{0TA , GT * 1 * ) ETAs I . 

EXTDR={ EXPTSAE>F{ EXPTS)+n)(PTPAEXF( EXPTF) ) 

CTRBET= ( I .-RCCK )4<EXT0R 
IF( ETA , GT . .9 90S )DTRBET=0 . 

PTA -1 I t-ROCK) >»ETA 
FETA? = BT a 
EVS=i:FT A2 

SHLEI NLAYFl )=SHL INLAY )-BETA2«( SHL I NL A V ) - £ HLE I NL A Y F 1 3 ) =('F XT R N 
12?5 CCNTINUE 
EVE= 1 « 

EVACO=EVE*DR AW 
CVACn=tANAXl I EVACC i 1 *E-4 0) 

SMSATS=CSAT ITLElNLAYPDtFS) , 

SHG^OSATITGtFS ) 

C NEW RETA CONTINLEC 

EVE=8ETA2 

TEMP=SHSATS/TLEINLAYP1) 

GAMC = ZLN*RHS’*9 »8 + (1 «- I 1 , + CLH* TEI^F/KAPA )/(l*4GAYFAC*TEKP/ 


HYD03310 
HYD03320 
HYDC3330 
HYD02240 
HYD03350 
HYD03360 
HYD03370 
HYD03380 
HYD03390 
HY003400 
HYD0.7410 
HYOC3420 
HYD03430 
HYD03440 
HYD034S0 
HYD03460 
HYD03470 
HYD03AQ0 
HYD03490 
HY0035CO 
HYD03510 
HYD03520 
HY003530 
HYD03540 
HYD02550 
HYD03560 
HYD03570 
HYD03580 
HYD03590 
HYD03600 
HYD0361 0 
HY003fc20 
HY003630 
HYD03640 
HYD03650 
HYD03660 
HYD03670 
HYD03680 
HYD03690 
HYD03700 
HYD03710 
HYD03720 
HYD03730 
HYD03740 
HYD03750 
HY003760 
HYD03770 
HYD03780 
HYD03790 
HYO03S00 
HYD03810 
HY003820 
HYD03830 
HYD03840 
HYD03850 



FILE I635TEST SSSCSN 


A 


*** GGODAfiC SPACE FLIGHT CEMEF **♦ 


* TLE( NLAYPI ) n>*K AP.A/PGAS 

TL.r( NLAYP 1 >= (ORA«*Td+E0V* ( Th ( NL A V ) * F S K -,C- A HC ) ) / ( C F A ^» + 1 D V ) 
SHLFIMLAYFI ) s:( EVAC0*SHG+E0V*SHL (NLA\ ) )/ (EV.ACC + ECV) 

SHSATS--CSAT ( TLE (NLAYF 1 ) ,P£ > 

IF{ SHLEINLAYFI J.LE.SHSAtS ) GO TO'lZC 
. GA!v*S=GAMFAC’!!SH£ATS/TL£<NLAYPn*^2 
TE,Yp= (S1-LE( NLAYFl )~SHSATS )/fl . 4CAMS ) 

TLE( NLAYP J )=TLE( NLAYP 1 )+TEMP*CLH 

IF(FHAJ )TLE( NLAYP n=TLE( NLAYPl ) - .7 £ * T E N F*CLh n # j} * 4 «>)c »i< * * * a ** 
SHLFI NL AYFl ) =Sl-LE (NLAYP 1 )-TEMp 
l?0 FMOS = PS/( FGAS + TLE INLAYPl > ) 

SFNSIRLE HFAT (LY/DAY) ANC EVAF0KA1ICN <Clv/C f'’*'*p/CEC ) 

OTS= ♦ 5* ( DTS+TG-7LE( NLAYPl ) ) 

nSI-5 = St-G-£HLE<NLAYPl ) . . 

IF(DTS*CE.0, ) CPA\s=C0*(WVAG4SCFT{01S) ) 

IF {OT5.L7»0») ^RAW = CO*WMAG»(WWAG*V^IVAG)/((V^^/AC+Vl^AG )-7.*DTS) 

FSURF=CCEFS4RHCS*CRAW/SP 

TL(NLAY ) = TL INLAY )4FSURF+DTS 

C * 4 4 Jc 

CCNl ■= T^-'(NL «Y )4PSK-GAMC- 
CCN? = DRAW*7G ' 

TF( IGVET«EQ. 17) TL(NLAY) = TL(NLAY) 4 F SLPF*DTS 4 C ♦ 5 
C SAVE CHANGE IN' TEMPERATURE IN LAYER »NLAY» CALSEC BY SURFACE 
C INTERACTION HEATING 

V(If9,M2) = FSURF+DTS 

IF( IGVET.E3. 17) V(I,9,M2) = V(I,S»M2) 4 FSLRF * 01 £ +C.5 

C S\VE EVAPORATION IN WM H2C 

- PIItM) = FSLRF1<EVE4DSHS'?'I00 .*D£IG (NL AY ) + SP/GRAV 
C*4A* IF(T HAJ )P ( i ♦ M )='» 25*P ( I , M 1 ) 4 4 4 4 44 * 4 *4 4 * 4 *** * « 

TH(NLAY )=TL(NLAY)/'PLK(NLAY) 

SHL( NLAY)=SHL(NLAY ) 4FSURF*EVE*DSHS 

C4=*AA I F( FMAJ )SHL ( NL AY ) =SHL ( NLA Y )-.7£*FSLRF4EVE*0£HS 4 4 44*44 4A.-1'4*4<’;' 

IF(SHLl NLAY) ,LT.2«> SHL(NLAY) = C* 

C SAVE EVAPCRATICN ANC SENSIBLE HEAT FLUX FCR DC E A N { L Y /C A Y ) 

TF ( ,NCT, OCEAN) GO TO 122 
EVAL =600. 

TEMP= 10 -.*CRAVi*FHn£'««SDAV 
V(I ,3,M3) = EVAL<EVe>i'TEMP*DSHS 
V{ I ,4, V7)=CP*TEMP*DTS 
C SAVE PRECIPITAELE WATER VAFGR( G/CM**2 ) 

122 PRECWV=0. 

CO 1 23 L=1 . NLAY 

123 FRECWV =PFF.C ft V 4 S HL ( L ) *DS I C ( L ) 

V(I, l,M2)=SP4opecftV*lO*/GRAV 
FSURFO= .2E>^FSUPF ‘ 

CC 1 AO 1 ♦ 2 
IF (J.EC.JSP) CC TO 135 

G (I , NLAY .KCENT ,K ) = 0 ( I . N L AY . K C E N T , H ) - F £L RFC *03 ( I , J C BN T . K ) 
">(t,M+K) = F(IiMl+,<) 4 FSI RFO*G£ ( I ,^CENT , K ) 
0(IMINC.NLAY»KCEN7,KJ = Q(IMINC.NLAY.KCEM»K) - FSLFFG* 

1 CS( IMI NC , JCENT . K ) 

P( IMINC-.M 14K) = P( IMINC.Ml+K ) F S L F F G *Q S ( I M I N C » JCEN 7 , K ) 

135 CONTINUE 

IF(J*E0.JNP) GC TC 140 

O ( I .NLAY .KUP ,K ) = 0 ( I . NL A Y ♦ KUF » K ) - F £ U F F C4 C £ ( I » JL F . K ) ' 


HYD0.3860 

HY003870 

HYD03880 

HYD03890 

HYDC3900 

HYD03910 

HYD03920 

HYD03930 

HY003940 

HYD 03950 

HYD03960 

HY003970 

HY003980 

HYOQ39«0 

HYDC4 000 

HYD04010 

HY004020 

HYD04C30 

HYD04040 

HYD04050 

HYD04060 

HYD04070 

HYD04060 

HYD04090 

HYD04 100 

HYD04 110 

HY004 120 

HYD04 130 

HYD04140 

HYD04150 

HYD04 160 

HYD04 170 

HY004 ISO 

HYD04 190 

HYD04200 

HYD0421 0 

HYD04220 

HYD04230 

HYD04240 

HYD04250' 

HY004260 

HYD04270 

*JYD04280 

HYOC4290 

HYD04300 

HYD0 4 31 0 

HYD04320 

HYD04330 

•HYD04340 

HYU043S0 

HYD04350 

HYD043.'0 

HY004380 

HYD 04 390 

HYD04400 



FILE I635TEST SSSCSN 


A 


«**.GCDOAfiC SPACE FLIGHT CENTER*** 


Cd 

) 


U(I.K»N3)'= L(I.K,M3)I + FSURFC * G S < I » JLF ♦ K ) 

Q{ I VI NC .NLAY ♦KUP tK.). = Q( I VINCi NLAY ♦>< UP * K ) -FSURF C *C £ ( 1 V I NC t JUF .K > 
U( IVINC.K ,V3) = U< IMINC.K.V3) 4 F 5L P FC *QS ( I V I NC » JUF * K ) 

140 CONTINUE 

SOOT{ J, I , 3> = T0.TALF + Sbdt't j;i »3) 

C:**** ELIMINATE NEGATIVE SPECIFIC HUVJDITV 
00 410 L=£.NLAY 

. LM = L-1 ■ ■ 

IFISHLC LV. 1) »GE .0* > GO TO 4 1C 
SHHL)=SHL(L >4SFL <LMl )*DSIG(LV1 )/D£ IGIL ) 

SHLILVl ) = C. 

410 CONTINUE 

IF(5HL( NLAY) «LT.O<. > SHL{NLay)=C. 

IFISHLE {NLAYFU .LT*0» ) SHLE I NL A YP 1 )= C . 

C***4 R/JOIATION SUEROUTfNES 
IFINFLW *EC»0 )GC TC 4?0 
IF ( *NCT .FLAGS ( 1 1 ) GO TO 420 
XDAY=FLC AT( JCAY462) 

XLAT = ( J- .5*{ JSP4JNP ) ) *DL AT* 1 ?a ./P 1 

LATO = XLAT 4 ,£ 

CALL SGLAFl I NL # Y «XOA V , XL A T ) 

OSAPEA = CXYF I J 1*RS0 I N 

IF ICDSZjGT.sOI) PROJAR = PR0JAR + C£AREA>1CGSZ 
TCSASS = £G 
nC 412 L = 1 .NLAY 
412 TDSAtiS = TDSAe£+AS(L) 

TGTAOS = TOT AG£+TCSABS*OSAREA 
CALL LINKhO ( NL A Y ,NFL lA , JO A Y .L A TC ) 

C SAVE SOLAR RADIATION INCIDENT ON TCP CF A T NCSPHE RE I L Y ✓ D A Y ) ' 

C*v** V ( I , 5 . M 3 ) =S0 X*CCSZ •<:♦+*** It *♦*:»*:»##* Cj_ q CARD REP L ACE D » * ♦ *'* >» 4 ** ♦ * * 
VII,5,M3)= AVAX I ( C* , SCX*CCSZ ) 

C SAVE SOLAR PAOIA.TICN ABSORBED 8Y ATMC SPHERE ( LY/RA Y ) 

V( I . 6 . V3 )=T CSAES-SG 

C SAVE. 30LAA RADIATION A6SURBED 8Y G R C L N C < L Y /C A Y > 

V(I ,7,V3)=SG 

C SAVE NET LONG~\aAVE PADIA.T.ICN AT TOR CF VCCFL A T V.O SP H ER E ( L Y /'D A Y > 

V{ I , e. N3 ) = RE{ 1 ) ■ •• 

C SAVE NET LONG-ViAVE RADIATION AT GPCUNC ( L Y/CA Y ) 

V( I , 9 , V3 > = RE( NL AYPI » 

ShS( J, I )=SG“RE(NLAYPl ) 

C * * SAVE WIND VAGNITUOE 

V( I .9 , n ) = W VAG 
C 

C**** PREDICTION OF CRCUND CONOITICN 
420 CONTINUE 
C ' NEW 

a-TSTO = WET 

W = TIT = V«ET*3.333'2333 3-«e6eeeEee€ 

WFT IT=2 ,+W5JT IT 

WET = A, VI Mdo.WETIT) ' 

IFIwET .LT .0. )WET = C. • 

C SET DEFAULTS FOR CCVP3 TAPF 

V{ I,5.V?) ■= C« ' “ " ■ 

V( I ♦ 6 . V2) = C . 

V(I,7,V2) = C« 


HYD04410 
HYD0A420 
HYD04430 
HYD04440 
HY004450 
HY004460 
HYD04470 
WY004480 
HYD04490 
HY004600 
HYD04510 
HY004520 
HYDC4530 
HYD04540 
HYD045S0 
HYD04560 
HYU04B70 
HYD04580 
HYOfi4S90 
HYD04600 
HYD0461 0 
HYD04620 
HYD04E30 
HYOC4640 
HYD04650 
HY004660 
HYD04G70 
HYO04eSO 
HYD04690 
HYD047C0 
HY00471 0 
HYD04720 
HYDC473b 
HYD04740 
HYD04750 
HYD04760 
HYD04770 
HY004780 
HYD04790 
HYDC4800 
HYD0481 0 
HYD04820 
HY004830 
HYD04840 
HYD04850 
HYD04860 
HYD04870 
HYD04880 
HYD04890 
HYD04900 
HYD04910 
HY004920 
HYD04930 
HYD04940 
HYD04950 



FILE T635TEST 5SSCSN 


A 


♦ GCDDARC SPACE ‘FLIGHT CENTER 


bd 

I 

H 

O 


V( I * E • ) - C ♦ 

T( I , 2.P2) - C . 

TF( OCEAN } GO TC 6^0 
IF(SNOW) GO TO € 1 C , 
r F( ICE) GC TC‘ €2C 
IFCFROST) GO TC 630 
. IF( M rXW I ) GO TC 6A0 

cVAL=60C* ’ , . 

CZH:=SCRT( (»3ec+-. 15*WET)*( 1 • + WE T ) * 2 ♦ E - 3 ^ S C A Y/ T VfiGP I) 

C7H=SQRT( ( .34VSy ) *( 1 • + I 0* VSM ) * 1 . OE- C 3 S C A V /T WCP I ) 

CZ»-'-CZH^( 1-RCCK5 4P0CK’^6« 7 
GC TO 660 
610 eVAL=:66C. 

CZH=2,3 
GC TO 660 
620 EVAL=68C* 

C?H=:n • 1 
GC TO 660 
630 EVAL^r68C. 

CZH=:5 0PT( (•331 + .C75^WET)#(2*42*6^\aE1)41 ♦ 3* SO A V /T W C P I ) 

C2H=CZHJf ( l-RCCK ) 46 •T^ROCK 
GC TO 660 
640 EVAL = 6BC-i 
W1:=GT(J^I) 

W I . . . , 

W T=2 > 4VV I 

WI = VM ^66 66^6 66 6- *66 6 666 6 66 

WI = 3 * I 

WI=:AM IN 1 ( 1 * * WI ) 

wT^AMAx i( c* ♦ y n 

IFlWET»GT*l*E-6) EVAL=6 00*4B0*5*yi/yET 

C2H=S0RT( (^2764(» U4-»15^WET)i{ WET - *6^y I ) ) ^< ( 1 * 4 W E T- W I 4 l . 26 W I) * 

^ E-T«SDA Y/T WCP I ) 

ir( ( GT ( J. I ) *LE * 1 * 1 ) ♦ AND* (LAND ) .AND* ( *NCT *SNOW ) ) CF2=CH2* ( l-RCCK ) 

1-^6«7AR0CK 
6HO TEM-^, 

NEW 

WET==WTSni 
W T = GT{ J * I ) 

t F< ICE* ANCoZ *LT • * I ) T EM = C T 1 0 /F I CE 
TGS0=TG4TC 

nPAD = 4* *STBO*TCSQ*TG 
DS0G“54 1P**shg/TC5Q 
TFMP= 10 * *CR AV4|;h0S^^SDAY 
NEW BETA TERN 

IF< EVE * CE ♦ 1 • CE- AC >0 TR EET= C TRBE T / E V E 
IF(EVE*LT * 1 *06-40 )E)TROET^C 
TF(EVF*GT#*909<; VCTREET^Oo 
OSQG-nSQQ^{ I i4CTREET) 

DTRFET=C 

TFRMl SG^NHOGAN 
TERM? = -CP4TENF*CTS 

IF( I GVET* FO* 13 ) TERM2 = TERM2 ^ 1*6 
TFRM3 “ -EVAL^EVE+TEMP'^OShS 


HYD04960 
HY004970 
HY0049WO 
HY004990 
HYD05000 
HYD05010 
HYDC5020 
HYD05030 
HYO0bO40 
HYD05050 
HYD05060 
HYD05070 
HYD05C30 
HYD05C90 
HYD06 ICO 
HY0061 1 0 
HYD06120 
HYD05 130 
HYD05 140 
HY005 150 
HY005 160 
HYD05170 
HYDCSiaO 
HYDC5 190 
HYD^E 200 
HY0D521 0 
HYD0h220 
HY006230 
HYD06240 
HY006250 
HYD05260 
HY005270 
HYD05280 
HYDC5290 
HYD06 30 0 
HYDC631 0 
HYD05320 
HYDC5330 
HYDG6340 
HYD06350 
HYQ05360 
HYD06370 
HYD06360 
HYDO 5 300 
HY005400 
HY00641 0 
HYD05420 
HYD05430 
HYD06440 
HYD06450 
HYD06460 
HYD06470 
HYD06430 
HYD05490 
HYDCSbO'^ 
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FILE leiSTFST SSSCSN- A 


+ GCODARC • SPACE FLIGHT CEMEF ♦ + * 


w 

t 


C**4* IF(FHAJ )rERM;:=TePKT+»25 + + 

TFRM4 = TEM* <T ICe-rGr 
TERH5 = -FF( NLAVP U*NHOGAN 

DENOM = SDAY*C2H/crC3 + DRAG - TEN 4 TFNF*EOV* ( CP/ ( DFAW4EDV ) + 

1 eval*eve*dsc 6/ (evAco+Hov ) ) " 

RADTRN=SHS( J t I ) ' ■ ' 

GT(Ji I )=TG4 ( RACTRN + TeRM2+TEFM3 + 7eRN4 )/CENCN 
IF t GT(J,n .C-T» 100*) GC TO 775 
WRITS (El E70 ) J, I ,CT( I ) 

WRITE (6,774 ) J , I , { TL ( L ) , L= I , NL A Y ) , T L F ( H L Ayp 1 ) , TG 
b70 FGPMAT( 10.x . J I GT(J,I) * * + A’) • 2 I 5 , E I 5 , 7/ / ) 

774 FPR.M at ( IX ,21 : . 1 IE 1 1 »3 ) 

77‘5 CCNTINUE 

C SAVE EVAPCRATICN. SENSIBLE HEAT FLLX , AHC CCNCfCTICN 
C THRCUGH SEA ICEILV/DAY) 

V( I , N3 )=-TEFN2 
V( I , 4, W3 )=-TEFN2 
V( 1 ,2, N3 )= TEFN4 

C SAVE GRCUNO TEMPERATURE ADJUSTMENTS 
V(r,5,N.?) = TERMl/DENCN 
V(I,6,H2) = TEFM2/CENCN 
V(I,7,N2) = TEFM3/bENOM 

VCI,8,N2) = TEFM4/DENGM 
T(I,3,N2) = TEFME/DENON 
fj90 CONTINUE 


CRCUND HETNEES EECTICN-GISS ♦*»**:^* + K. 


SNFLX=TERN?^« e<,G4 4A44 4 444 4 6-04- 
SNFLX=-SNFLX‘ 

RADTRV=SHS( J . I ) 

RNIT=RACTPM*C„S4444444444F-04 ' 

eTFLX=TERM3+e,SA4A444444E~04 : • 
ETTLX^-ETFLX 

IF ( OCEAN ) ETF LX = V { t , 3 , M3 ) + e« 94 4 < 4 A4 444 AE-C4 
IF(OCCA N )SNFLX=V { I , 4, M3) ♦ C , R4 4 44 4444 4E -C4 
GNFLX-RM T-S NFLX-ETFLX+TEPM4*6,S4'4 4 444 44E-C4' 
EeDT = DSH£*EVE*CTC3i'0RAW*RH0S« 1 C .C 


WETR='AET 
£NR=0 . 

W IR = 0 » 
TGR=TG 

2550 RAIN=0. 

SNFAL=0 . 
F RFH 7-50 , 


HYDOfsblO 
HYD05520 
HYDOSSiO 
HY005540 
HY005550 
HYD05560 
HYD05570 
HY005580 
HyD05590 
HYDOSfcOO 
HY,00‘J610 
HYD0562Q 
HYD05630 
HYD05640 
HY005650 
HY005G(>0 
HYDCbS70 
HYDD5630 
HYD0G690 
HYD0570D 
HYD0571 0 
HYD05720 
HYD05730 
HYDC574 0 
HYD05750 
HY005760 
HYD0b77C 
HYDO’5780 
HYDC57Q0 
HYOC5EOO 
HYD05FI0 
■•~HYD05e20 
HYD05830 
- HYO0Sa40 
HY005G50 
- HYD05860 
HYD05870 
HYD05E80 
HYD05890 
HYDCrs900 
HYDC591 0 
•HYDC5920 
HYD05930 
UYOC594 0 
HYD0595r> 
HYD05960 
HYDC5970 
HYD05980 
HYD05990 
HYDC6000 
HYDC6010 
HYDC6020 
HYD06030 
HY006040 
HYD06050 



rvnnri ronnnnnnonn n r>nn 


Flue I635TEST SSSCSN 


A 


*♦* GCDDAPD SFACe FLIGHT CEMEF ♦> * 


RUNOFF=0* 

SNELT=0 . 

FLDRK=0 . 

□FNG=0» • 

RLNOVV=0 . 

NEW 
N EW 

IFIOCEAMGO TO 26E0 
TGR = GTI J, I ) 

IFITCTALP.LT ,.0 )TQTALP = 0, 

IF(TLE( NLAYP 1 ) *LE »T I CE ) SNF AL = TC 7 At F * I ^C* 

TF{ TLE( NLAYP I ) .CT 4T ICE) RAIN=7C7ALF >U C C • 

I SNR = 7F IX ( GW IJ t n ) 

SN REFLOAT ( 1 SNR )/ 1 COO.,0 
IF(SNR<.LT<.«0C1) SNR=0* 

SNR1 = SNF+SN'F AL 
SNR-SNR 1-E6CT 
I F( SNRl ftLE^O CC TO 26EC 
TF( SKP .LE »0 * > CO TC 2650 

IFITGR «LF.4TI Ce )GC TO 2600 
S KELT=C2H’i-(TCR-7 I C£)/80« 

SHELT = AHIM ( EHELT ,SNR) 

TGR=TGR-3C«*FNFLT/C2H 
S^R = SNR-S^'ELT 
I F( ICE) GO TC 557f 


SOIL (VOISTURE ANC FREEZING 


SNOW COVEFEO LANC 


NEW .... 

600 CONTINUE 

WETR=:WFT+(RA 1N4£VELT)/WGM 
RUNOFF=AMAX 1 1 0 .C ♦ WETR- I • ) 

WETR=WETR-RUNCFF 

GO TO 2580 

NEW 

SNOW FREE LANE 

2650 SNR=0« 

in ICE) GC TO 2E7C! 

C NEW 

IF(RAIK,LT-.0.)RAIN = C*C 
TFI RA IN .UE 1 1 4E-eO )G0 TO £140 
A0INF=1 .25*( *WFI-.2)-.5 

FEF ACT= 2.*4F ILTF*rXP{ - U7E*( ADI NT- , e ) ) 
IFIPArN,Cr»BEFACT )RUNOFF=RA 1N-.E*EEFAC7 
1 r( RAIN iUT^ REF A CT )RUNCFF= , 5* R A I N* R A I N / E E F AC 7 
Fla'' c r K T ! M 1 1-: 


HYD0606D 
HYD06070 
HYD06080 
HYDC6090 
HY006 100 
HYDQ6110 
HYD06 120 
HY006130 
HYD06I40 
HYD06150 
HY006 160 
HYD06170 
HYD06180 
HYD06190 
HYOC6200 
HYD0621C 
HYD06220 
HY0C6P30 
HYD06240 
HYD06250 
HYO06260 
HYD06270 
HYD06280 
HYD06P90 
HYD06300 
HYDC63I 0 
HYDC63.20 
HYD0e330 
HYDC6340 
HYD06350 
HY006360 
HYD06370 
HYD06330 
HYD06390 
HY006400 
HYD06410 
HYD06420 
HYD06430 
HYD06440 
HYD06450 
HYD06460 
HY006470 
HYDOfcCiaO 
HYD06490 
HYO06S0O 
HYDC6510 
HYD06520 
HYDOeSOO 
HYO0GF40 
HY006550 
HYD06500 
HYD06870 
MYD06580 
HYD06F90 
HYD 06600 



FILIt I635TEST SSSCSN A *♦* GCDDAFC SP.ACE FLIGHT CENTEF 


w 

1 

M 

W 


ADDVi'=K A TN-RUNQFF 
ADDVi'=ADC’A>f-( 1 »-PCO<) - 

FCAR = E6Cr/( 1 •- PCCK) 

•aETR=V»ET + (ACCW“ECAR)/1«GM 
RLNOVV=FCCK*PA r^ 

RUN0FF= FUNOVV + (1-ROCK )*RUHCFF 
ORNG=0. C5 55-AM AX 1(0«0» WETR- *5) 

WETR=WETR-DRNG 
C NEW 

I F( WETR «GT<, 1 • ) ViETR= 1 * 

I F ( WE TR «L T . C -. > V« E T R =0 . 
r F(FRCS7 ) TGR = TGR + P0•*AOD^^/CZH 
IF( .NCT .M IXVv I ) WI=0. 

I F[ FRCST ) W I=ViETR 

IF( *NCT *(A'1XVII «CR.(FR0ST.AND-.TGR.GT.TICE).CR*(TG*G1 .TICE .AND. 

1 TGR «LT »T ICE ) > )GC TC HEEC 

FREEZ = C7f-=HTlCF-TCR)/'80,/V.GH 
'A rR-=Wt + FREEZ 
IF( W I R.CT .WE TR ) Vv I R = W£TR 
IFIW IP *LT»0. )W lP = r. 

FREFZ = W IR-W I 

TGP=TGR-fFO.*FREEZ*WC-M/CZH 

C 

C ICE 

C 

C ICE 

C 

C 

GC TO 2 eec 

2670 TGR=Ayi M (TGP.TICE) 

WCTR= ! « 

26EC TF( IGVET.GE. 10 ) GT(J.l) =. TGR 

IF(WIR«GT..O .ANC.WIR.LT.WETR) GT(J.l) = VIP 
SNR = AN INK £NF .2C.0 ) 

ISNR = IF IX( AV/IX 1 (Oo , SNR ) * loco. > 

IF( IGVET.GE.20 ) GW(J.I) = FLCAT(ISNR) 4 C.50+WETP 
C*****s/\VE GRCUND ViE TNE SS ♦ RUNOF F . SMEU T 
TI 1 , 1 .M2 ) = GW { J. 1) 

T(I,2.M2) = SMELT 
T(I .3. M2 ) = FUNCFF 
T ( I ,4 ,M2 ) = CT ( J. I ) 

T(I.5.M2) = TCTALF+ICO. 

TGN=GT( J . I )-T ICE 

IF{ RSURFF .UT..Cl)TGM=TGM + e 

IF{TGN.LT.~270. )TGN = GT( J, I ) 

TPMN=A , ’JGW{ J . I ) 

XRS=RHS 
PCT ALP=RA IN 

GWSWCH= AM('D( ( T AL-TAUI + GSW )/GSW . I . ) .L T . ,C CO 1 
I FIOCFAN )GWShCH=.FALSE> 

IF( J.LF . IF )GWSW CH=,FALSE. 

IF{ J ,GE .271GWSWCH = .FALSE. 

IF(I.GE«46) GW£WCH=*FALSE * • 

IF{GWSWCH.ANO.(l/ITC*ITC*EO.I.AND«<J-HnC) / ITC’* ITC .EO. J4 ITC-1 ) ) 
1 GO TC GOES 


HYD0661 0 
HYD06630 
HY006630 
HYD06e40 
HYDC6650 
HYD06660 
HY006670 
HY006680 
HYD066R0 
HYD0670C 
HYDC671 0 
HYDC6720 
HYD06730 
HYD06740 
HYDC6750 
HYD06760 
HY006770 
HYD06730 
HYD067<70 
HY00680C 
HYD0681 0 
HYD06820 
HYD0683C 
hYD06340 
HYDC68S0 
HYDoeaao 
HYD06870 
HYD06830 
HYD06890 
HYD06900 
HYD06910 
HYD06920 
HYOC6930 
HYDC6940 
HYD0695C 
HYD06960 
HYD06970 
HYDC6930 
HYD06990 
HYDQ7000 
HYD0701 0 
HYOC7020 
HYDC7030 
HYD07040 
HYD07050 
HYD07060 
HYD07070 
H YD 07 080 
HY0C7C9r 
HYDC7100 
HYD071 10 
HYD07120 
HYDC7130 
HYD07140 
HYD07 150 







nnnnn 


FILE I675TEST SSSHSN 


A 


*** GCDDARD SPACE FLIGFT CEN7EF *** 


I3d 

4S 


IF{ TAU* EG *TAU I «ANC . I .GE. 1 E. AND< I*LE.22.ANC.J*GE.31.AND,J.LE«36) 
t GC TC FOES 

Gc TC aces 

8055 CCMINUE 

CHANGED UNITS GF' PREV AP (MW/HOUF TC CM/HAIF hCUR EY D I VD BY 20 
PFEVAP=FPFVAF/2C . 

WRITE (6 « 2 551 ) J. I » GW { J ♦ I ) » GT ( J V I ) ♦ F Cl A L F -.F LNCFF , FF E V AP , UCT A 2 , 


) HYD07160 

HYD07170 
HYDOViaO 
HYDC7190 
10/1HYD07200 
HYDC7210 
. HYDC7220 


lECAR* F I LTF. BEF ACT tRHOS, DRRN*QSSH*DFAW ,CSHS»PLANT ,FCCK,R0FI6V\, DTSIHYD07230 
2551 FCRMATI 1X,2I 3 .BXtCcei 1 ,4, jx ) ♦/ , ICX (E 1 1 ,2X n HYDC7240 

2555 FCRMATI //^tax. • J I W E TNE S S /RHC S CFTENF/DRPN FC1ALF/CSSH FUN C F HY DC-7250 

1/DRAW PREVAP/DSHS BETA2/FLT ECAR/RCCK F IL TF/ROF I CW EEFACTHYDC7260 

2/CTSI !♦///) HYD07270 

P0C5 CONTINUE HYD07280 


* « I»* If ground V»EINESS SEGTICN-ENCS *>|if«*)f*** 


C*A*fc TOTAL HEATING AND MTISTURF ADJUSTMENT 


DC 520 L=],NLAY 

SH( [ ,L« KCENT ) = SHLIL) 

C SAVE CHANGE IN TEMPERATURE FRCM SOLAR PACIATICN HEATING 
TM.L.Ml) = AS(L)*COE(L> 

C SAVE CHANGE IN TEMPERATURE FRCM LCNG-WAVE FADIATICN HEATING 
SH(I,L,M1) = <FE( L+1 J-RE (L ) ) fCCe (L 1 
520 T(T,L*KCENT) = TL (L ) + ( AS ( L ) +RE < L 4 1 )-Rf { L )) *CCE { L ) 

TS( J« n =TLE( NL AYP ! ) 

SDPTIJ.I.n = ami M ( SDOT( j', I ♦ 1 ) -.TSIJ ♦ I ) ) 

SDOriJ,!,?! = AKAXU SDOT( J, I ,2) ,TS (J . 1 ) ) 

61'^ IMINC=I 
RETURN 
END 


HYOC7350 
HYD0736C 
HYD0737Q 
HYD073R0 
HYD07390 
HYD07A00 
HYD0741 0 
HYD07420 
HYDC7430 
HYD07440 
HYDC7450 
HYD07460 
HYD07470 



si-*a 


om; 


'15Tm5TW 

COMMON 


TOTF57TFinnci:¥7Tirrr‘zr 

PHIS (08,14) ,FD(0a,14) ,TOPOG (08, 14) ,SINL(46) ,COSl(46) 


$ 

T 

$ 


PC (8,14) ,TSC (8,14) ,SHSC(8,14) ,G3?C (8, 14) ,GWC(8, 14) , 

UC (8, T4 ' ;3) ' • ,VC ( ’ 87T473n ' T ^T B7T47 37T5HClF7T473TTP'rK'rr0)' 

SGC(8,14) , OS (08,14) ,VS(08,14) 




~w 

* 


PLE(IO), TLB (10), BH (10) ,TH (10) ,DSIG (10) ,SHL (10) , 
S18L"B'C10) ,Bel ( 1 0) ,5H'{8', 14 , J) rFITTOTT I 

GI(08,14) ,SHS (08, 14) ,T1 (10) 


SLBEOO (08,14) , Tg ' CWn^) ' 
IHB (10) ,SIGE (10) ,SXG (10) 


UMTO tr;T3-)”7BTjirTJ87i"^' 


OLM00010 
OLM00020 
OLM0OO3O 
OLM00040 
, OIMOO'OSO' 
OLH00060 
OIM00070 
OLM00080 
"CniHOOO'90 
OLM00100 
UI'HOOTTO'“ 
OLM00120 
■■'OXHOO'I 30 
OLH00140 
0TM‘00150‘ 
OLM00160 
OLM00170 
OLM00J80 
olm56i'90 
OLM00200 


~w- 

* 




#yHus,C'z 

* 

DIHENSI 
1 


2 

3 


EHS, DTS,CD,SDV,ZLN ,DRAM,PaSV&P,PLaNT,E VE,B6DT, 

H , L'TH'B'ET , TEKTriVi’-HIB M2 ,TBB M3 , TljyTlTr7T‘Him5",”'T5E NCmVB'A'WBTT,' 
SNFaL,BAIN,OVES,HSTB, WET,TAO,HMAG_,E^ 

ON PLTGTC (200) ,PLTGT (200) ,PLTGWC (20O)'7pX'IGW (200) , 
PLTSHL (200) .PLSHE (200) ,PLTX (200) ,PLTIS (200) 


PLTPB (200) ,PLTAF (200) ,PLTEVE (200) ,PL18GS (200) 
PIEGCM (200) .PLT9TG (200) .PABGO (2Q0) 


L03ICAL + 1 STAB/** V#PLOS/* + V»POUND/*#V» CROSS/* X' /» PSHIOD/’ 
LOGICAL O C SAN^SNOW.ICE.LAND^MlXWI.FaoST.GCH 


'/ 


-1.1 - hjlxwx , rauax , 

UCLA , UCNI,UCNII ,UCiJJJA, A3Sia,TINPT, VALBFG 
WETCSI , A7CSIL ,DRYCEL ,LOG8 , UCLAVW, V5P4, GDH 
a t>a 


OLM00210 

OLM00220 


C 

C 


LOGICAL 

LOGICAL 


BFAL44 KAPA 
INTEGER »2 SLBEOO 
QSAT(X,Y) = 


BEOO 

.622* (10.«*(9.4051-2353./X) )/Y 


REVISED MODEL STARTED 1/78 
E CMS HYDRO NEW AT GSFC 


OLM00230 

OLM00240 

OLM00256 

OLM00260 

OLM00270 

-Q.WlPO.2eo 

OLm6029O 

OLM00300 


C 

_C_ 

C 

C 


C 

C 


DATE 


CORRECTIONS 


STEP 1: DESIGN OF EXPBRI-MSNT 


GCM=-TBUB. 
GCM-. FALSE. 


A SET xMODSL LOGICAL PABAMET35S 


OLM00310 
OLH00320 
OLM00330 
OLM00340 
OLMOOSSO 
^L MOO 3 60 
OLm6037O 
OLM00380 


■OFF-ILINE MODEL (joLM)] 



B-16 


c 

UCLA=.FALSB, 


OLM00390 

OLM00400 


■ UCLA=.TEU^. 


bLH60410 


OCLAVW=. FALSE. 


OLM00420 


UCLAVW=.TRUS. 


OLHObitSb " 


DCNI=.TSUEi 


OLM00440 


UCNI=, FALSE.. 


OLH00450 


OCNJJA=. FALSE. 


OLM00460 


.UCNII = .THUE. 


dLM6047b 


aCNII-. FALSE. 


OLM00480 


ASSTM=. FALSF. 


OLMO 649 O 

c 

TINPTrTfiOE IP SURFACE MD SPEED INFO IS AVAILABLE 

OLM00500 

c 

ON INPUT DATA TAPE 

* 

OL MOOS 10 


TINPT=. FALSE. 


0LM00520 


LOG8=. FALSE. 


dLM00530" 


LOS 8=. TRUE. 


OLM00540 

c 

LOGS TRUE, OLM WILL 

READ AND USE TAPE VARIABLES, 

OL MOO 5 50' 

c 

GWC,GTC, EVERY TIME 

STEP 

OLM00560 


VALBFG=.TEUF. 


OLm00570 

c 



OLM00580 


HETCSL =. FALSE. 


“~OLM00 590" 


WETCSL=.TRUE. 


OLK00600 


AVCELL=.rRUE. 


'OL'HOOeiO 


AVCELL=. FALSE. 

. 

OLM00620 


DRYCEL=.TR“nE. 


'OLM00 630 


DPYCEL=. FALSE. 


OLM9*)640 


VRP4=.TSUE, ' 


b'LM006'5~0 


VRP4=. FALSE. 


OLM00660 


GDH=.TPUE. 


" “OLM0T670 

, 

GDH=. FALSE. 


OLM00680 




■'OLK006'9'0 

__c 



OLM00700 

g*™ 

; 3 PROGRAM' CONSTANTS 


OLMOO'7’1'0" 

_c 


• 

OLMO072O 


NLAY = & ' ~ 


OLK007'30 


NLAYP1 = NLAY + 1 


OLMO074O 


JMIN=1 


OLH0075'0 


JMAX=8 


OLM00760 




■ 5lm0‘077'0 



n n o 


I a 1. 

jarN0s=3O 
'JM A X US^TT 
IMINnS=12 
IMAXUS=25 

' nlTymi=nlaT-T 

JWP=JHAX-1 

I!lC=TfiTN 

JNC =1 

'lM=iaAX 

J1P=3MIN+1 

’JHPP1=JMP+1 


C DO loop' PAEAM'F;.TPPS""foE Sp'?cfFIc''ciLL" LOCATIOKS 


11 

12 


13 

J1_ 

J2 

J3 


= rMIN+1 
= I».AX_ 
=“ INC' 

= jaiN 


JXAX-1 

JSC 


= Q« 

DO 7 0""l='1 ,NLAY 
nSIG (L) =« 1 11 1111 1 • 


SIGE (L+1) =STGE (L) +DSXG (L) 

70 SIG ( X) =.5* { SIG S (L) ^-STG S(L +1) ) 


PGAS=287- 
KAPA=.28G 
GRAV = 9.8 
■ffCE = 273'. 16 

■ "f WO'PI” =''"2 . * 3 1 4 1 6 
PI=3, 141592654 
DPI=PI 
CTTD = 432. 



■ 

o 

. 


■S'-rs- 

^ & 
•*«. fn 

y 

— 




’oia00780 
OLH00790 
OLM00800 
OLMOOSfO 
OLH00820 
OLM00830 
OLM00840 
OLM00850 
OLB00860 
6 lm 60870 
OLH00880 
OLM00890 
OLB00900 
0LH009i0 
OLB00920 
OLMOOPSO 
OLB00940 
OLM00950 
OLM00960 
OLM00970 
OLB00980 
OLH00990 
OLM01 000 
OLM01010 
OLH01020 
OLB01030 
OLBOI 040 
OLB01050 
OLB01060 
OL«01070 
OLM_0_10J0 
OLH01090 
OLB01 100 
0LM01 110 
OLM01 120 
OLM01 130 
OLB01140 
OLM01 150 
OLB01 1 60 



o jO n a n n n ^ n a n ri ci \ • o o o ol n 


HIC^'= 300.' "• 

STBO = 1. 173-7 • • ' ■ 

'-Cp'*= 75 a " ■ ““ * 

CLH=600./.24 

'GAM FAC =’“!TrT'F^LH 

a yap ah 3 F FnNCII0N"W*TlM3"TT3P ' 

'■ WC3=‘”TB7)'0'r'“^ — — • ' 

NH0GAN=1 

rpTA'$STffr"’E’r^3^3T^:'*'3T •" 

IF(ASSTH) NH0GAN=1 

^ j[^t«-7-|}-^S~jr--3- -jjrjiju- ■■Tia^™r'imFTSX 

C03FS=G5AV*D TC 3. /DSI G [N^Y) 

B DATA SPECIFXCATION FOE IMITIAL RONS 


total? = 0 . 

"flU'NT''H'=3 

IGVST=1 5 
IGVST^T”~ 
IGVFT=23 
■PTEOT^si 0 . 
FD=10. 


C INITIALIZATION OF CELL VARIABLES 

EOK=1: OCEAN, 2:L0AMY SOIL, 3;SANDY SOIL 

" • "8i;si)‘{ 5 ; Trsy • ( PVT)' iu^~, 8) ,r=i / 1 4 ) " " ' 

SLBEPO FOE JULY I527ALBD BUT GC« TAPE DS3D14 
P.7AD (5, 146) ( {SLB3D0 (J,I) ,J=1,8) ,1=1,14) 

BUT FOB THRE E ‘TEST C ELLS IT READS PHIS ARBAY 

SEE TOPOG PRINT OUT • 

SEADfS. 145) ( fPHIS(J.I) ,J=1,8) .1=1.14 ) _ 

DO 148 L=8, MONTH ' 

D O 148 1= 1 ,14 _ 

'"“’"read fs, 14 7)" m 7 (TOPOG (j/l) ,J=1,8") 


OLN01 170 
OLH01 180 
OLM01190 
• OLM01200 
0LM01210 
OLH.01220 
0LH01230 
OLH01240 
OLH01250 
OLM01260 
OLM01270 
OLM01 280 
OLM01"290 
OLM01300 
■ '6LM01310 
OIM01320 
0LM01330 
OLM01340 
OLM01350 
OLM01360 
OLH01370 
OLM01380 
OLM01390 
OLM01406 
OLMOiaiO 
OLM01420 
OLM0143O 
OLM01440 
OLM01450 
OLM01460 
OLM01470 
FOR LANOLM01480 
01M01490 
OLM01500 
OLM01510 
. _ OLM01520. 

OLMOl 530 
OLM01540 
OLH0155O 



:n o' n o|n o‘n oi'o of |0 { in oi jo oin 1 { 1 nn 

.< i i r • { I • I : I. r 


1il8 GONTINTL, ^ ; 

1U5 ''format (8 (PI 0.5')') •" ■" 

14 6 FOP.HAT f 8 fI2) ) . 

’147 F0SMAT{I8,8.(F9.3)') 

G BKAD(20’) SLBSDO'.PHIS.TOPOG ^ ^ 

D SST SOWS FOE ^ LL AVSEAG BS f !1W 0. 

AVGT=0/" ' “ ' 

S U!3AV =6. . 

AVTS=0". ■■ ‘ ~ - — 

A/GW = 0.. ;;; 

AVGTC=0. 


OLW01560 
0IW01 570 
^LM01580 
OLfJOVsOO 
o;lmoi66o 
■ bLM016’‘l0 
OLM01620 
O1M01630 
OLM016U0 
OLfibl 650 
OLH01660 
brab"f67o 


AV TSC ^O. _ - OLM01680 

'AV3WC=0. - - “6L«bl690 

OLK01700 

~ ‘ " P INITIAL DATA TAPE HSAD"“ “ ” " OLM617i6 

; OLM01720 

READ (1)8,95) 'TAO,PC,TSC,SHSCrGTC“,GWC,OC, VC,TC,SHC,FLOXC " OLM01730 


9 5 FO RM AT ( 0 26 (255 A4) ) 


HOOP. = TAU ' • 

NHH -1Q 

NHH=96 

l'I_NC ^5 

'"l^’(ASSIM) TINC=3.0 

D b“'20o" ih'h"=”T,"n h'h 
. H0rjE= HCOR+TINC' 


f" BEGIN CELL LOOP" OF'mGDSL'' 
VARIAB LE LOCA TION 


OLN017U0 

■ OLM01750 

OLM01760 

0LM01776" 

OLM017^0 

6l«01790 
OLNOI 800 
OLH01810 
OLM,01820 
OLN01830 

OLM01840 

•" OLM01850 

OLK0 1860 

OLH01'870' 
OL H0 18 80 
' OLH01890. 

OLH01900 
OLM01910 
TL,T(L_AYER 90LM01920 
'OLK01930 
0LH01940 



w 

1 

M 

O 



IS^TL'2/PL'S OLM01950 

DLa01960 

^ sMc;jc 4c GT\P / TS^^'S HS ^ OL H 0 1 9 7 0 

OLM01980 

OLD Ww ‘ SQUIVALInT SOUFCE ‘ NAJ3B ‘ " ' OLM01990 

OLM02000 


'"NOTS i 
TS(10, 10) 
sTs XTJr"i ->) 
GT (1 0, 10) 


N0T2 1 
TS (8, m) 
~5hs (1J7T07 
GI(8,1«). 

G» (8TJ4) 


PC 

TSC (8,14)-' 
SHSC (BTI^T 


IAP3 

TAPP 


NOTF 2 

•suTTar'— 

TL33(10) 
■^h'TTO] — 
PLF (10) 
IH (1 0) 
DSIG (10) 


DC{8, 14, 3) 

TnT5;Tr,iT 

NOTE 2 
3KTB7''nf7Jy 
TLP(IO) 


_GTC (8,14) 
GWC(8',T4') 

PC (8 ,14,3) 

TcTF7Tr,3P' 

TC (8,14,3) 

3iin'8TTr;3T 


TAPP 

T APF 

fAPF 

TAPE 


TAPE 

TAPP 


TAPE 

STEP 


HH 

PL? (10) 
TH (1 0) 
DS.IG(IQ) 


STEP 

STEP 

STEP 

P80G 


■5~ 

4 A 
7 


OLM02020 

"p {’^,1:)- '=' PI '^■ps-?OLa02030' 
SUHFA^ TEHPSBATDEOLH02 040 
"¥E T ■ E A D rA T 10 N ( L Y / D O' L M 0 2 0 5 0 
GBOOI^ _TEMPE_RATUBSOLM02060 

■"G BOUND" WSTNES'S' DLW02070 

X VSL COMPONENT (M/OLM02080 
COMfONENT (?l/OLM02090 
LAYER TEMPERATUBE OLM02100 
" S’PSC HUM OF "L'AYSH OLM02110 
T (L) OF INTER I-1,OLM02120 
■ RELAfiVE HUMIDITY" OLM0"2130 
PE OF INTER. FACS, LOLM02 140 
POTENTIAL TEaPEEATOLM021 50 


3LBED0 (10,1 0) SLBEDO(P,14) 


NOTE 1: PL & PLE ARE SPKCIFTSD IN 
NOTE 2 TL IS SPECIFIED IN STEP 4fi 


IF (WITCEL) 
IFf WETCFL) 


„O0_pNF CFLL FOR NOW 10J22/T7 

> 1"C= 1 2 

J1C=7 


IF(AVCSLL) 
IF (AVCELL) 

"if(d'byc?lT 

JFpRYC_fL) 

j3C=l“ 

l3Cf1 

I2C=I1C 
J2C=J1C 
NCFLL=0' 


I1C=6 

J1C=3 

I1C=4 

J1C=4 


CONST. =1/NL yr= . 1J 1 1 1 1 OL M02 1 60 
SPECIFIED 'earth ALBEDO PBOMOLM02170 

OLK02180 

step” 4.A FROM PC OLM02190 

IF. GCM. INPUT NOT AV AILABOLM02200 

OLH02210 

OLM02220 
OLH02230 

OLM02240 

OLK02250 
OLM02260 
OLM02270 
OLM02280 
OLM02290 
OLM02300 
OLM02310 
OLM02320 
OLM02330 


^0- 


i 



o o o o, onnoo » n nn ^ *no 


DO 300 ’=I1C^I2C,T3c 
DO 300 J=J1C, J2C, J3C 


NO='LL=NCFLL+1 


DO LOOP BEAD 


rF(IHH.GT.I) 5SAD (08,95) TA 0 , PC ,TSC,SHSC , Gf C ,G1JC, UC, VC,TC, SflC 
*,FLUXC 


IF(NCFLL.EQ.I) 


1CALL 3DBFWD(NLAY,NLAyHl,JKP,ZNC ,I K, J1 P, JMPP1 , TI SPT, ASS,I«) 


H SOLAS DECLINATION CALCULATION 


IF(NCFLL.FQ. 1) 

1C ALL 3DET (PI, DPI ,COSD , SIND , ROT, DION , ASSIM , IHH, LAT ,JDAY, NCELL , 
■fTOFDAY) ■ . ... 


I OFFLINE INTEEFACE CONSIDERATIONS 


INITIALIZS GW(J ,1 ) WITH G i iC (J ,I) , _ V AL UF 3 

I?(IHH.FQ.1) GW (J,I) =GWC (J,I) 

■■ "iJTr~RB'SD"TWlC'B"'BSClUS’E'“Sc‘r'fSDTATION NOT AVAILABLE 

IF (GTC ( J, I)_.LT. 2.) GTC(J,I) =TTCE . 

lf’.'2) "‘^T (J’,'iySGTC (J,'I)“ - - --- - - 

TS (J , X) =TSC (J,’’) 

■3'?lx'oi37?~A'?m7Tirtnn~^G¥TT;i)'4r ~ 

IT='(LOG8) GT (J,I) =GTC (J,I) _ 

_J ZENITH ANGIE CALCULATION 

IUS=ININUS+T-1 

-JUS=JNTNDS+"J=n ‘ ' 

HACOS=COSD»COS(SOT+ (IOS-1) *DLON) 

COSYl'^^Sr! 

Z =AECOS(COSZ) 


QLa02340 

OLK02350 

OLH02360 

OLK02370 

OLM02380 

OLH02390 

OLM02400 

OLM02410 

OLM02420 

6Lil0243"6 

OL«02440 

OLM62450 

OL.N02460 

OLK02470 

OLM02480 

OLH02490 

OLH02500 

OLM02510 

OLM02520 

OLM02530 

OLH02540 

OLM02550 

OLM02560 

OLM02570 

OLM02580 

IHOIN02590 

OLM02600 

OLH02610 

OLM02620 

OLH02630 

OLM02640 

OLa02650 

OLM02660 

OLa02570 

OLH02680 

OLH02690 

OLH02700 

OLN02710 

OLM02720 



bd 

1 

to 


C 

C 

C 


Srira-=COSir*SITT (?OT+mS-l) *D101T)/SIN'(^) ■ - OLM02730 

CHOriP=TOFDaY- (37. -luS) +2tt./72. OL«027i»0 

— -:n?-xi:injirR'.'LT;"DT} "•cirTJF"='t:HTmFy2^'r - -- - • olh 02750 

IF(NC'E!LL.?Q. 1) WRITF(6,152) COSZ, SISA , CHO0H OLM02760 

t572-FOB7ia'Tnoor;-^so'LiTr"z^rrH awrFr"co'sz“'*=‘*","Trr2;Q,5x, ’izrjio-TH auGiis (+»olho2770 
1,--3 OF S) , SINA = • ,312.4, 07X, ’LOCAL CELL HOUH = ’,F6.2) OL«02780 

„ . . . .OLM02790 

STEP 3: SURFACE SP2CIFICATI0SS OL«02800 

OLM02810 

3A SURFACE COKDI.TldN OLM0282P. 

OLM02830' 


OCEAN=TOPOG f J .11 . GT. 1 . 

IC3=T0P0G (J,I) .LT.-9.9B5 

LAND-. NOT. fIC3.0R.0C3AN1 

FQ.1) GO TO 54 ' • - 

. 53 SNOH=LAND.ASD.GW (J.I1 . GT. 1.Q 0 

IF(IGVST.SQ.24) GO TO 52 

_ _ SNOW = LAND. AND. (GW fj . II .GT. 1. 0 ^ B- .J LMMi J * .401 
IF(SNOW) GW (J,I)=AMAXi (10. 5, GW (J,I)) 

GO TO 52 

54 SRH=85. *GW (J,T) +15.' " “ 

_ I F(P.QK(J,I) .20 . 1.) G O TO 55 _ 

IF(HOK(J,X) .EQ.2.") GO TO 44 

IF (H OK(J, 7) .2 Q.3 .) GO TO 42 ' __ 

44 PC-725 ‘ ‘ “ " ' 

. PHP = . 10 

C2=1 . 0 

'?OFE=.S 

SO ILLB = 1 QO. 

WCAS=POEE*SOILLE ' 

GO TO 49 

4 2 FC=,09 ' "" - — - - - 

PMP=._02^_ _ 

cY=i."o 

; . _ c2=i .0 _ _ _ 

I ' ' P0H3=.5 

_ S OIILR^OSO. _ _ 

WCaS="PORE*So‘lLLR “ “ ’ * T' • 


OLtl02_8_4Q. 

6lm62850 

DL.f302860, 

OLH0287'0 

07.M02880 

OLM02890 

OLM02900 

OLN0291 0 

OLM02920 

OLMOiSSiO 

OLM02940 

6lM02950 

OLM02960 

OLM02970 

OLH02980 

OL«02990 

OLM03000 

OL«03010 

OLf303020 

OLH03030 

OLflO3O40 

OLn03050 

OLM03050 

OLH03070 

OLM03080 

OLN03090 

OLM03100 

OLM03110 



w 

1 

hO 

Ij) 


C 

C 


l i9 GW (J,I) =SKH» (FC*C 1-P WP*C2 )/130- +PWP*C2 
"5S GW {J,D =GW (J, I) /2'' 

if(oct:an) gw(j,i) =.5 
IF {IGVBT.EQ.24) SO 3?0"21 
IF (OCEAN) GO TO 52 

P.FF T&K CMP12980 


PT180 = P7/180. 

SNOWN=(6Q.-15.*COS(.9863» (JPAY- 24 .6 6B) » P1 180)^) *PI 1 80 

SNOWS=SNOMN-2.*PI/3. 

JNP=46 


C 

•L 

C 


C 

r 


IF (LAT (JnS) .GE.SNCWN) GW (J,I) =GW (J,I) +IST( (LAT (JUS) -SNOWN) / (LA 
*) -SNOWN) *3.54) 

"IF (LAT (JUS)". LE. SNOWS) GW (J,I) =GW (J,I) +iNf ( (LXT "(JO Sf-SNOMS) / (L 
♦-SNOWS) *3.54) 

GO TO 52 ■ " ■ 

21 IF(SLBBDO (J,I) .LB.39) GO TO 52 
GW (J,I) =2.00503 +7T 
IF (SLBEDO (J,I) .LT.70.) GW (J ,1) =1 . OE03 +.5 
IF(SLBSDO (J,I)'.LT.50) ■ GW(J7J) =.3I"0'3+."5 ■ 

IF (ICS) GW (J,I) = _^.E03 + .5 
■■ lF("JUS.‘LE.87AND.LANDy“W('4',r) =30VOF03 + .5 
SNOW = LAND-ASD.GW (J,I) .GT. 1.01 
"5 '2 ‘"iTTvTlbfg) go ■ 

SL33D0 (J ,I) =14 

IF (OCEAN) SLBSDO (J,"i) =7: " ‘ 

SNOW = LAND.AND. (GW (J,I) .GT. 1. J.OF.SLBSDO (J,I) .GE.40) 

IF ("SNOW, of: ICE)" SLBEDO fJ,!) ‘‘= 7’0,' ‘ 

56 3SUEF=. 01*SLBEDO (Jrl) 

^ 

SN0TE5=SQET ((GW (J, I) -.5) /I 000. ) 

■ ■ SSO'PF'^,'5^NOTEH~ 

IF (SNOTFM.LT. 1.) PSU3F=RSU5 F+. 0 1 *SLB3D0 (J,I) * ( 1 SNOTFU) 

”F;sTn " amtn"t (s^'uEF,.8) 

66 WET^AHINI (1. , (AMOB (GW (J^I) ^1.) *2) ) 

IF VALBPG IS T MUST HAVE FOLLOWING CABB 11/2/77 
SNO W = LA ND. AND. (GW(J,I) . G.T. 1 . 0 .OF . SLBEBC (J, II .65.40) 

LAND=LANB. and".. NOT, SNOW 

MIXWI=GI (J,I) ,LE. 1, 1.AND. (LAND.OR.SNOK) 


OLf103120 
OLN03130 
OLH03140 
OLH03150 
OLM03160 
'“OLM03170 
OLM03180 
OLW03190 
OLS03200 
OLH03210 
OLH03220 
T (JNPOLM03230 
0L503240 
AT(1) OLN032'50 
OLN03260 
OLM03270 
OLH03280 
OLN03290 
OLM03300 
OLH03310 
OLN03320 
OLM03330 
OLH03340 
OLN03350 
OLM03360 
OLM.03370 
OLN03380 
OLM03390 
01^103400 
OLM03410 
OLM03420 
OLM03430 
0L«03440 
OLM03450 
OLM03460 
OLN03470 
OLM03480 
OLM03490 
OLH03500 



w 

1 

IS3 


fsost=land!Tnd. (GT (J#D .lfTtics) . Tnoi. «ixwi 

, '^L EV?/tI ON ’ ^ 

c 

IF(OGBAN.AND. (PHISrJ,I)"sQ.O.) )‘ Z=0".'~ " 

__ iF-f.NOi«ocFAH) z=ph ; s (J.I^ /GRAV •_ ... . 

C " 

__ C_ ....... 3B SORFACB T F MPESATUKS 

c" 

IF ( OCEAH) TG=TOPOG (J,I ) 

IF“(“InOT.OCFAN) 'J!G=GT(J,I) 

IF( Hi rrfl ) TG = T ICF 

IF(: not.'land) 'wet=i."o ' ■ 

c 

"c ■ ■ 3C SURFACE 'OFaG COFFFICIFNT 

C _ _ _ _ _ 

TMINC=I-1 ' “ 

W?l AG = .5*SQRT{F g(J,I) +FD {a,IHrNC) t FP ( 1 ,1) •«• FD ( J+1 ,IMXNC)) 

TF(OCEAN) ~C'0= AHINl ( {1~. +'. 07 =f'W~aAG)''i. 001 , .0"025) 

IF (. NOT. OCEAN) CD= . 0 02 +. 006*Z/5 JOO. 

"“c “ 

IF (IGVET.EQ.15) GO TO 71 _ _ 

‘■‘iF{GDH. ANprOCSAN) Cb=CDC5' 

IFCGDH.AN D. . NOT. OCEAN) CD=CD#2.0 

‘GO‘t6~72" ■■■ ■“ 

I 71 CD=.OO32*(1.+3.*Z/5D0O.) 

•“IF (OCEAN) CD= . 1 3/ (ALOG ( 1 . 0F067T5. 56/lMAXf {1 . ^'w «aT)”“. 98^+ 

1 .00687*WMAG*W.1AG) ) ) **2 

■72'CONflNbS 

C 

c ' ' ■" 

C STEP U: VSF.T. ARRAYS FOB SOME VARIABLES 

C' 

I C 4A PBFSSUHES 

' C 

SP=PC(J,I) 

C’ ~ PC (J,I)"=ps~pTFbP=i'I=P (J,I) =sp 

C P ( J, I) =PTBOP+SIG (NLAY) *PT 

.‘“C PS=PTRGP+SIG (Tt)l * (?S-PTPOP)=PS-~ ''' 


• OLM03510 
OL.M03520 
OLM03530 
OLM0354O 
‘OLM03550 
OLM03560 
OLH03570 
OLM03 580 
OLM03590 
OLM03600 
OLM03610 
OLH03620 
■ OLM03630 
OLK03640 
OLM03650 
OLH03660 
OLM03670 
OLM03680 
OLM03690 
OLM03700 
OLM03710 
OLM03720 
OL MO 37 30 
OLM037U0 
OLM03753 
OLM03760 
OLM03770 
OLM03780 
OLM03790 
OLH03800 
OL MO 38 id 
OLM03820 
OLM03830 
OLM03840 
OLM03850 
OLM03860 
OLM03870 
OLM03880 
■■■■OLM03890 



H 

I 

Oi 


: c 


PS=SP+PTHOP=PLE (NLAY + 1) 

OLM03900 

, ^ 


DO 30 L=1 ,NLAY 

' OLN03910 

, \ 


PL (L) =SIG (L) *SP+PTROP 

OLM03920 


30 

PLS (L) =SIGE (L) «'SP+PTROP 
PLK(L)= PL(L)**KAPA 

OLM'OE'^SO 

OLH03940 

c 


PLj?Xm Y+ 1) =SP4;p.TROP ~ 

OLH03960 

c 


ITB temp'^RAtuRbs and cobvbcTIve ad^ITsTHent 

• OLH03970 
OLM03980 



DO 35 L=1,KLAY 
NC=L-6 

“ ■ ■ ""OrrMOBTSU' 
OLM04000 


35 

J. f ( L • 11& xL{l) =tc (xc {kJ f 1. f z) d 4 J. f 1 ) ^ ( NIj AY~2'*LJ'1 

IF(L.GE.7) TL {L) =TC {J,I,HC) 

tjrH’owoTir 

OLM04020 



IF(GCH) GTT' TO 4-1 ' 

DO 40 N = 1 ,3 

OXBWTTJCr 

OLM04040 



DO Tnr'i;=2,NLA^ 

C5"I7K1040"5^0 


L«1=L-1 OLH04060 


j rr rn.ir- 1 ) /l->lk (x- i| " .« E ':Tnx) ‘ 7 PiK (T. ' ) ) c;(r T r tr " 4 " U ‘OXfi 04 “DT 0 

I TRSTA=(TL (LM 1 ) *DSIG (LH 1 ) +TL (L) #DSIG(L) ) / (PLK (LM 1 ) ♦DSIG (LMl ) + OLM 04080 

i «PLKTO'#DS1G (L) ) ' ^OrHOtra'90 

; TL (L- 1 ) =TIIBTA*PLK (L- 1 ) OLM 04100 

; T L ( L ) =i‘ HiiTTl* P L KIT) OIWITI TO" 

I 40 TH(L-I) =TL(L- 1 ) /PLK (L- 1 ) 'OLH 04120 

! 41 TH (NLAY)=TL (NLAY) /PIK (NLAT) OLW 04130 

‘ C OLM 04.140 

C 4 C HOISTOH 2 VABIABLSS ' ” OLMOflSO 

j DO 50 L -1 .NLAY : OLM 0416 Q 

NC=L -6 OLM 04170 

IFg.LE. 6 ) SHL fP =SHC (J, 1 . 1 ) - fSHC -SHC fJ.I. 11 » rNLAY- 2 -L ) ) OLM 04 1 80 

IP(L.GS. 7 ) 5 HL(L) =SHC (J,I,NC) OLM 04190 

• I? fSHLfLI ,LB.O.) SHIfl 1 = 0 . 0 L MM 200 

50 CONTINUE OLM 04210 

IF fGcm GO TO 175 ; OLM 0 42 2 Q 

C OLM 04230 

_ C STEP 5 QLM 04240 

'C VEHTICAL DIFFUSION OF HEAT AND KO ISTURE (AIS- EARTH ■^"bLH 04250 ’ 

INTERACTION OLM 0426 Q 

C ‘ _ - OLH 04270 

65 LMIN=NLAYM 1 OLM 04280 



92-a 




DO 80 L=LMIN,NLAy 

OLM04290 



LM1=L-1 

OLM04300 



DTBTA={IH (LH1)~TH (L) ) ♦ (PLK(LCll) +PLK (L) ) *.5 

oTM0"43ib 



DZOP=SP*DSIG (LM1) *RGAS*!TL (LMl) / (PL (L«1) *GRAV) 

OLM04320 



D2DN=SP*DSIG (L) *RGAS*TL (L) / (PL (L) »GRAV) 

OLM04330 



EDLE=-2#ED 

OLH04340 



TEt5P=DTC3*( (DSIG (LH1) +DSIG (L) ) / { ( DZUP+ DZDIT) ♦ (DZOP+D ZDH) ) ) 

OLM04350 



FLP=-2.»EDlE*D1?ETa*TEHP 

OLM04360 

- 


TL(LMI) =T1 (LMl) +FLE/DSIG (LM 1) 

6lV04370" 



TL(L)=TL (L)-FLE/DSIG (L) 

OLM0438O 



TH (LBI) =!TL (LMl) /PLK (LMl) 

OLH04 390 



TH(L) =TL(L)/PLK(L) 

OLM04400 



DSH=SHL (LMl) -SHL(L) 

OLM04410 



ELB=-2, *EDLE*DSH*TEMP 

OLM04420 



SHL (LM1)=SHL (LMl) +ELS/DSIG (LMl) 

6lh04430 


80 

SHI (L) =SHL (L) -ELE/DSIG (1) 

OLM0444O 


175 

CONTINUE 

OLMb44’50 

c 



OLM04460 

c 


STEP 6: RELATIVE HUMIDITY 

0LH044'70 

c 



OLM04480 



DO 90 L=7,NLAY 

OLM04490 



TSHP=QSfiT (TL (L) ,PL (L) ) 

OLH04500 


90 

RH (L) =SHL (L) /TEMP 

6LM045lb 

c 



OLM04520 



IPE=1 

‘ "0LHb4530" 



I?(NCSLL.EQ, 1.AND.IHH.EQ. 1) 

OLM04540 



♦ CALL RITE (IWRS,I,J,IHH,HOUR,IPS, 

"‘OLM04550 



♦UCLA ,UCNI,UCNJJ A, UCNII , I2C, J2C, IMA X ,aM A X,IMINUS,J MINUS , MONTH) 

OLM04560 

c 



OLMOasTo" 

c 


STEP 7; DETEP.MIHATIOH OF SURF TEMP 

0LM04580 

c 



dLM04“590 

c 


A INITIALIZATION 

01H04600 

' "c 



' OLM056'1Q" 



TLE(NLAYPI) =TS(J,I) 

0LM04620 

, 


PS=PLE(NLAY+1) 

OLMO'46'3d 

J 


PSK=PS**KAPA 

OLM04640 

1 


THG = TS(t7,I)/PSK 

OLH04650 

c 



OLM04660 




B RH SCALS 





3 


EHS=0. 

IF (RET+EH (NL?.Y) .NS.O) EHS=2. *WET*SH (KiaY) / (WET+RH (NLAr) f 
C 

_ C bulk aehodynahic cobff. 

c 


c 

c 


c 

c 


OLfl04680 

olmO'46'90 

_OLM04700 

orffd4Tin 

OLM04720 


■t)'XS=T(3-'rreTmTPT) 

IF(DTS.GB.O.) DHAW=CD* (WMAG+SQST (DTS) ) 

TF(I>TS.UTy:)' DRA«=CD*WMAG*(ttMA'G'#W«XG)'/TTraS=i;wMAGy"'- 

D SUBFACB DIFFUSION COBFF 

ZLNCO = .5*DSIG (NLAY) *5Gas/GRAV 

SP^fPLB nO) -PLf9) )/f1.-17./18.) 

E ORIGINAL FORaOLATION FROM TSANG 6 K ABN 


WMAGN=1 


— m'M(nr7w 

DLN04740 
orHW75U 
OLM04760 

dl¥OTT7'0' 

OLM04780 
OLH04790 

j[y,H04800 

' OLM048fd 

OLM04 82.Q, 

OLMd483d 

OLN04840 

0LMd4 8*50 
O1M04860 
OLM04870 

O1NQ4 880 

O1M04890 
: OIM04.9Q0 

C OlMOif’dTd" 

SDNS=10.-f10Q.» ( 1,-SXP (1200.* ( (TH (NLA Y) -IHG) / ZLN) ) ) OLM04920 

“C ■ ' TH(9) GT TUG THE FOLLOWING EDNS PBEDICTION IS BAbLM04930 

IF(TH(9) .GT.THG) SPNS=EOBIG OLM04940 

pc ■ “ ■" ' OLM 04950 

BDV=EI)NS/2LN OLH04960 

■■ ■ » t n ut , i ii i 1 1 ii i wi i m i i— i ^ <p<um 9* ^^1 

C OLH04970 


IF{DTS,LT.O.) WMAGN = . 2* (WMAG+ . 01 ) 
EDNS=AHIH1 (100. ,ED»EXP (. 32*DTS/WMaGN»HaAGN) ) 


FOHIG^EDNS 
ZLN = ZLNCO*TLE (NIAYP1) *SP/PS 


NEW EDNS FOFHULATION (1/30/78) 
NCAB W/NO COrjNTBRGBADIENT 


.^•l C WEITE(6,7123) OLM04980 

C7123 F05MAT(10X,’BOBIGS6X,'EI)NS’,8X, 'TH(9) »',8X, 'THG'»)"“ ' OLRd4990 

_C W BITE(6,7119) E0RIG,3 DN S, TH (9) ,THG OLM05000 

C7119 F0BHAT(5X,4 (B12.4) ) ' ' " OlN050‘l0 

C OLM05020 

•'■'"c ■" ' ' ‘ . . ■ '-oLgosdio 

EVE=1. OLH05040 

"" IF(LAND.AND-,NOT.SNOW) GO TO 7034 ~ ‘ OLMOSOSO 

GO TO 7235 OLM05060 



8Z-a 




-• 7034 

] C 

CONTINUE 

olm'osoyo 

OLH050S0 

c 

F EVAPORATION SCALING FACTOR 

OLM05090 

C 


OLM05100 


IF(OCNI) EVE=AMIN1 {1. , (HBT-PWP*C2) / (FC^CI -P»P*C1) ) 

OLM05110 


TEMP=2. *GWC(J,I) 

OLM05120 

t 

IF(UCLA) EVB=AMIN1 (1. ,2.*TSHP) 

OLH05130 

1 

PLANT=AHIN1( 1. , COSZ/. 2588) 

OLM05140 

i 

IF(COSZ.LT.O) PLANT=0. 

OLHOSlSO 



OLM05160 

' 7233^ 

CONTINUE 

OLl3’'0'5170 

! 

GO TO 7236 

OLH05180 

TJTT 

EVE=1.0 

■ QLM0519'0 

7236 

BVACO=EVE*DRAW 

OLH05200 


FVAC0=AMAX1 (EVACO, 1.E-40) 

OLHO'SSfO 


SHSATS=QSAT (TLE(NLAYPI) ,PS) 

OLM05220 


SHG=QSAT(TG,PS)' 

OLM05230 


TEMP=SHSATS/TLE (NLAYPl) 

OLM05240 


GAHC=ZLN^EHS*9.8* (1.- (Tr+CLH*TSMP/KAPA) /TT. +GAMFAC*TEMP/ 

OLM0S250 


* TLE (NLAYPl) )) +KAPA/SGAS 

OIM05260 


• TLB (HLA^r f ) = (Dgi ' ro s yFpy^TrH (MLA ' yy » ' P5 T C-gracTr ^^^ ^ 

■ ! SHLE (»L AYP1) ={EVACO*SHG+EDV*SHL (NLfiY) ) /(EVACO+EDV) 

^ (TLB (NLAYP1 )' 

IF(SH1E (NLAYP1) .LE.3HSATS) GO TO 7121 

' S=Gk HP AC*'3HSATS/T1T (NIA YP1 ) *'*2 • 

TEMP= (SHLE{NLAYP1) -SHSATS)/{1+GAMS) 

T L E { IfL'A YF'Tr-TTE '(iniirrP’Tr'-F ' 'T EHP^UIH 

SHLE<ULAYP1) = SHLE (SLAYP1) -TEMP 
" •”7T2T"'ErHUS=P3 / {KGA3^'i'LE'{KL'ATP'’n"y 
C 

-"C G FOTENTrAI E V APUH'OTaNTfimTHEl 


t3'LT'{T5'2ro “ 

OLM05280 

•DLaD52'90' 

OLM05300 
0005310 
OLH05320 
'• 'GLHOSaBO' 
OLM05340 

0‘Lm)5350* 

0LM05360 

D0033TO‘ 

OLH05380 


PLE ' V AP^ ' DHAWgEHUb* VOO T ^ «j3AT (TG , P3 ' ) ~'S Hi rB~ (i r LAY y ~'TOr'fi 05300 

IF{PHSVAP.LS. 0. .AND.EVE.lt. .P99) GO TO 7237 OLH05400 


OLH05U10 

C H POT EVAP USING LOGS INPUT (MK/HE) OLK05420 

— c xyotmECT 

I T=GTC(J,I) OLM05440 


DGCM=GTC (J,I) -TSC (J,I) . . „ OLM05450 



62-^9 


r 

HBTGCM=GWC (J-.Il #2. 

OLM05460 


BVCM=AHi"Nl (1. ,2.*WETGCM) 
SGCM=QSAT (T,PS1 

"6lH05470 

OLM05480 


IF(DGCH.GE. 0.) DRCM-CD* (MMAG+SQRT (DGCM) ) 

IFfDGCM.LT. O.J DECM=CD*WMAG* ( WMAG*W!1AG) / ( (aMAG*»MAG) -7 . 

OLM05490 
*DGCM) OLM05500 


EVACOG=EVC«*DRCM 

OLH055i6 


SH10 =rFVACOG*SGCM+EDV*SHL (9) ) / (E VACOG+ SDV) 

OLM05520 


DGCMSH=SGCi1-SHlO 

OLm65530 


EGCM=DGCHSH*EVCM*DTC3»DRCM«RHOS^10. 

OLM05540 

C 

C 

STEP 8A: 

6lm05550' 

OLM05560 

c 

SENSIBLE HEAT(LY/DAY) & E VAPOR ATXON (GE/C M**27SEC) 

OLHOBSYd 

OLM05580 



DTS=-5<‘ (DTS+TG-TLS (NLAYP1) ) 
DSHS=SHG-SH1E (NLAYPI) 

OLM05590 

OLM05600 


IF (DTS. GE. 0.) DRAW=CD* (HHAG + SQRT (DTS) ) 

IFfDTS, LT.O.) DRAH=CD* (WMAG*=«'3) /( WMAG*=«'2-7 . *DTS) 

~OLM056lb 

OLM05620 


FSURF=COSFS*EHOS*DEAW/SP 
TL{NLAY)=TL (NLAY) +FSUEF*DTS 

OLH05636 

OLM05640 


TH (NLAY) =TL (NLAY) /PLK (NLAY) 

'■ OLMb565'0 


C 

c 


SaVSHL^SHL (NL&Y) 


OLM05660 


SHL (NL-ay) -SHL (NLAY) + FSUPF*EVE*DSHS 
SAVSHL=SfiVSHl/SHL (NLRY) 


IF(SHL(NLAY) .LT.O.) SHL (NLAY) =0 . 


IP? = 2 


CALL RITE (IHRS,I,J,IHH,HOUE,XPP, 


*DCLA ,DCNI ,UCNJJA, OCNII ,I2C,J2C, IMAX,JMAX,IMINUS,JMINUS, MONTH) 


HTSTO=WEI 
1405 IF (OCEAN) 


IF(SNOW) GO TO 1410 
IF (ICB) GO TO 1420 


IF(FBOST) 

IF(HIXWI) 


GO 

GO 


FVAL=600. 
IF (UCLA) 


OLM05670 
011305680 
"'■■OLM556 9b' 
OLH05700 
OLM05710 
OLH05720 
■■■ OLHb573b 

_OLM05740 

STEP Viz GEOOND 0Pt>i~TE' OF MOISTORS~ A ND“T5MPOrM05750 

_ OLM05760 

OLM05770 
OLM05780 

~~~ ■■■ " ■ ■■ OLH'65790 

OLH05800 

OLHOSSlO 

OLM05820 
b'LM05830’ 

CZH= ( (.386 + .15*WET) » (1 . + HET) *2. B-3 »SD AY/TWOPI ) **.5 _OLB505840 


GO TO 1490 


TO 

TO 


11130 

1440 


B-30 



r — 

1410 

GO TO 1450 
EVAL= 60 O. 

0 X 805850 

0 X 805860 

r 

CZH = 2.3 
GO TO 1450 

0 X 805870 

OLM 05880 

1420 

E 7 AL= 880 . 

CZH= 5.1 

oxtiOs’eTO 

0 X 805900 

1430 

GO TO 1450 
EVAL= 680 , 

6 X 805910 

0 LM 05920 


CZH=SQRT( (,331 + . 075 *WET) * ( 2 . + ^. S^^TET) ^1 . E- 3 *SDrvTVK)PI) ' 
GO TO 1450 

'“■ 6 L¥ 0 S^ 3 ir 

OLM 05940 

144 0 

: C 

EVAX= 680 , 

' olmOsoTo 

OL 805960 

C 

C 

A HI CALCULATION 

o'lhP’spTo 

OL 805980 

1445 

irr=GT(J,I) 

CZH=SQET( (. 276 +(. 11 +. 15 *WET) * { WET-. 5 <=»I) ) * ( 1 . +WBT-HI + 1 . 25 *WI) * 

OLM 05 ^ 9 ' 90 ' 

0 X 806000 


* 2.E~3*^SDAy/TWOPI) - - - -- QLf^QgorO 

1450 TSM=0, OLM06020 


¥1T=WTST0 ' " ' ■" OLM0603"0 


(T 




C 

C 

C 

C 




WT=GT (J,I) 


B GROOND TEMPEHATOEE UPDATE 


IF(ICS.AND.Z.LT. . 1) TEM-CTID/HICE 
TGSQ=TG*TG 

DEAD=4. ♦STBO^TGSQ’i'TG 

DSQG=5418.*SHG/TGSQ 

TBMP=1Q.*D5AH*RHOS*SDAY 


TSEH1 = SG*NHOGAN 

TEBM2 = -CP*TEMP»DTS 

TERM 3 = -EVAX*SVB=*'TS!iP*D3HS 
TERH4= TE"1» (TICE-TG) 


Eai 5=.3 

BE (NLAYP1) = (i. 17E“7*TG**4) »EMIS 

TBEH5 == -SE (KLAYP1) *HHOGA N . 

DENOM = SDAY*CZH/DTC3“+ DBA'd + TEM^ + TEHP<=EDV *1cP/d’0.W+¥dvT 


OLM06040 

OLM06060 
'OLM06070' 
0X806080 
OLH'0‘6’090 
0X806100 
OLM061 10 
OLM06120 
OXM06130 
. „.QL806 140 
dxKoeiso 

^0X806160 
6X806170 
OLH06180 
OXM06190 
. —QJi^Joeppo 
0X806210 
^ _OLM06220 . 
+ 'ox 806230 ’ 



n 


L 


I 

1 

] 

I 

!■ 

i 


1 



*EVAL*l!:V3*DSQG/f^VBC0+EDV) 

3HS (J,I) =SHSC (J,I) ♦NHOGAN “ 

C_^ SHS(J^T ) = SG-R3 (MLAYP1) 

RADTR»=SHS (J,I) 

GT(J,I) =TG+ (gAD!^RW+TE5H2+TBSM3^•TB£M4) /DEKOa 

1490 'continue 

14 75 CONTINUE 

C ~ 

IPE=4 

C 


CALL RITE (IWRS,I,J,IHH,HOOS,rpa, 
*0CLA,UCNI,0CNJJA,UCNII,X2C,J2C,IHAX,JIiAX,IHINUS,JHINUS, MONTH) 


OLM06240 

'’o’lm 06250 
OLH06260 
■'OLM06270 
OLH06280 
"OLM06290 
•01H06300 
■ 0LM063‘l0 
01M06320 
1)LM06336 
Om06340 
'bLH06“350 


C C ACTUAL EVAPORATION (CM/HH) 


E6DT=DS HS*EVS*DTC3*DHAH*FHOS*10.0 


C D GROUND MOISTURE UPDATE 

C 

HETE=WET 

SNR=0. 


_OLM_06^0 
0'LM06370 
O LM06380 
blM06390 
_ OLM06UOO 
6 ih064’10 
_ __OLM06420 
‘ OLM06430' 
OLM06440 


«IR=0. 


OLM'06450' 


TGS=TG 
^'651 EAl'N=(i, 
SNFAL=0. 

PEEBZ = Q ; 

OVER=0. 

■■ SMBLT=0. 
FLDBK=0. 
DENG^. 


I 



BUNOVV=0. 

c 


TF (OCEAN) GO TO 18 8^' 
TGE=GT(J,I) 


"TbT]rLT-^nnrcT37r,"3j' 


THEN RAIN INPUT HOULD BE IN 


CM/S A MPLS TIME 
MM/SAM PLS TIME 


IF(TOTALP,LT..O) TOTALP=0. 


OLM0^460 

'~0LMb647b' 

OLM06480 

OLM06490 

OLH06500 

OLM‘06510 

OLM06520 

■"OLH06S30' 

OLM06540 

~bLtiO?"5"57)" 

OLM06560 

bLH06570 

OLM06580 

OLM06590 

OLM06600 

OLMD’66TO 

OLH05620^ 




* - £ 


K3, , 


SNR^FIOAT (ISNH) /1000. 0 

r F ’CS N K . LT . , 0 U 1 ) Sim^O . 

SNR1=SNS+SNFAL 

SNR=SNH 1-E6 l)T ^ 

IF{SNH1 .LB.O.) GO TO 1850 

iF(i>NH. LE.O":') Ga'U'O TBSO" 

IF(XGR.LE.TICE) GO TO 1800 

-5TTE1‘T=C'ZHJ«‘ (TGR-TICEr/aO. 

SMS1T=AMIN1 (SMELT, SNR) 

TGJi=TGH-3 0. ♦SHEL'i’/L’ ZH 

SNR=SNS-SMELT 
IF (ICE) GO TO 1870 


E SOIL MOISTUBE AND FEEEZING 


aO;4-''Kij:v!‘ ’=.>.vss^ . 


1 SNOW COVEEED LAND 


^%SW 


1800 CONTINUE 

WETB^WET^- (RAIN-fSMELTI /WCAS 


EUN0FF=AMAX1 (0.0, WETS- 1.) 
HETR=WETR-RDNOFF 



GO TO 1880 


2 SNOW FREE. LAND 


1850 SNR=0. 

IF (ICE) GO TO 1870 


EAIN=FLUXC(J,I,3) *100. 

1892 E6DT=DSHS*EVE=»DTC3*DRAW*RflOS»10.0 


o o 

-3L3a- 


T7 Q 





'OtMCf6'630 
OLM06640 
T5tlSr0-6'650 
OLH06660 
"UTHOTeTO ■ 
OLM06680 
“OLTrOr669(T' 
01M06700 

"ia‘Mcr6Tir ' 

OLM06720 
■“OrH0'67T0" 
OLM06740 
■••■OL"MD'6750" 
OLM06760 
OLM06770 
_ OLa.0^780 
OLM06790 
OT.MQ8 80 Q 
OLH06810 
OLM06820 


OLM06830 

OLM06840 


OLM06850 

0LH06860 


WETS=WST 

IF(RAIN.GT. 0.) OVF,B=RAIN»(.1 + .9*WET*»1.5) 


IF(UCLAVW) WgTB=HET 


IF(WETR.GT, 1.) WETS=1,' 
IF(M£TR.LT. Q.) WETR=0. 


MOISTURE CONTINUITY EQUATION 
+ (BAIN“-E6DT -• OVER) /WCAS 


OLB06870 
OLM06880 
OLm66890' 
^LM0690p 
OLM069“i6‘ 
OLM6920 
OLM06930' 
OLM06940__ 
Ol¥ 06‘950 
OLM06960 
OLiib6970 
OL106980 
OLM06 990' 
OLM07000 


ADDW= (EAIN-E6DT-0VER) /WCAS 


OLM07010 


















5?%S» ** *^-V'’^¥s^‘ 
•r .-?-•: . 

kiV *‘’r^'''^0^ »>T’*^ '^.'■^^S?' 










k*^- 'V,*- < 






a^-Tz * .? . '^‘'*^*i*S 

•“^ • " ' i> O' 

.•-..•= •.. ‘..' ''M 

♦, * *'' A'^s * , 

■ -T ^ X**v » 




IF (FROST) TGR=TGH+80.<'ADDW/CZH 
IP(.NOT.MIXWI) HI=0. " 

IF (FROST) WI=«ET5 

IF(.NOI. (MIXWI.OR. (FR0ST.AND.TGR.GT.TIC3) .OR. (TG, GT. TICE. AND. 

J TGR.LT.TICE) ) ) GO TO 1880 

FREEZ=CZH* (TICE-TGB) /80./WCAS 
WIR=HI+FRESZ 
IF(HIR.GT.HETR) WIE=WEIR 
IF(WIR. LT.O.) WIE-O. 

FREEZ=WIR-WI ^To 

TGR=TGa+8Q.»FREEZ»WCAS/ CZH 2 

GO TO 1880 !F® 


= TGE 


<o 2 
& ® 

■ " ■ ^ '-■ asz ‘ ■ 

*< m 


1870 TGS=ABIN1 (TGR,TIC2) 

WETR=1, 

1880 IF(IGVST.GS. 10) GT(J,I) = TGE 

IF(WIE.GT.. O.AND.HIB.LT.HBTR) 6T(J,I) = HIE 
SNR = AMIS1 (SNR,307(T) 

ISNH = IFIX (AMAX1 (0.,SNE)»1000.) 

IF(. NOT. LAND) G0“T0 300 
IPR=5 


CALL RITE (IWRS,I>J, IfiH, HOUR. IPR, 

*DCLA ,aCNI, OCNJJ A, DCNII,I2C,J2C, TJlAX,JHAX,iai NU si J MINUS, MONTH) 

SOMAV=SUHAV+l ^ 

AV3TC=A VGTC+GTC (J ,1) 

AVTSC=AVTSC+TSC (J,I) 

AVGWC=AVGWC+GWC (J,I) 

AVGT = AVGT + GT(J,I) 

AVTS=AVTS + TLE (10) 

AVGW=AVGW+G3 (J,I) 

PLOT DATA 


IH2=IHH 

CONTINUE 

PLTGTC(IH2) =GTC (J,I) 


ONLX 


0_LM07020 

■ olmoToso' 

O LM 0 704Q 
OLM07050 
OLM07Q6Q 
OLMOfoTo" 
OLMQ 7Q80 
OLM07090 
OLH07100 
"OLHO? 110“ 
OLMQ71 2Q 
OLM07130' 
OLM07140 
OLM07150 
OLH07160 
OLM07170 
OLM07180 
■“OLMOffPO" 
0LH Q720Q 
OLMO'7'210' 
OLH07220 
OLM0723"0 
OLM07240 

"OL'Hffrair'' 
qLM0726_0 
"OLH07270 
OLK 07280 
'0LM07290" 
OLM07300 

OLM07320 

‘mTOT3T)“ 

OLP107340 
■'OLM 07350' 
OLM07360 
■olmD73'70“ 
OLM07380 
OLH07T9O 
OLK07400 











iiM 




PLITS (I1T2) =T^S (J,!) 

PLTGT (IH 2 ) =TG 

PLTGWC (IH 2 ) =SHLE (NLAYPI) 

PLTGW gH 2 ) =SHL( 9 ) 

PLTSH 1 (IH 2 ) =BORIG 

PIS HE fIH 2 ) =EDHS 

PLT 9 TG (IH 2 ) =TH ( 9 ) /THG 

PLTX fIH 2 ) =IH 2 


IF(DCLA) PSSVAP= (S 6 DT/EVS) * 20 . 

PLTPE (IH 2 > =PBE 7 AP/ 2 Q. 

PLTAE (IH 2 ) =E 6 DT 

PL 3 ?EVE(IH 2 ) =EVE 

PLEGCM (IH 2 ) =EGCM 

PLIBGS{IH 2 ) =!TIB ( 10 ) . 

P AEGO (I H 2 ) =PLEGCfl (IH 2 ) /PITAS (IH 2 ) 

IP(IH 2 . GT. 95 , AND.IH 2 .lt. 101 ) IH 2 =IH 2+1 

IF(IH 2 . GT. 95 . AND.IH 2 .lt. 101 ) GO TO 2553 

CONTINUE 

CONTINUE 

TEMP=SUHAV 


AVGHC^AVGWC/TSHP 

AVTSC=AVTSC/TEMP 

AVGTC^A VGTC/TEMP ^ 

AVGW=AV GW /TEMP 
AVTS=AVTS/TE«P 
AVGT=AVGT/TEHP 


WRITE (6,2571) 

2571 FORMAT(6X, 'TIME STEP* , IX, *T6C ,6X, »TG* ,7X , ' iTC^r'SX , *TS (N + 1) 

»»SHLE(1Q) » ,4X, *SHL(9) 

*4X,‘E0RIG» ,3X,*EDNS‘, 6X, «PE > , 8X, • AE » ,?X , ’ EVE' ,6X, » AB-GCM • , 

*2X,»EG/EO») 

PQ “2574 i=i^ijhH 

WRITE (6 ,2572) I,PLTGTC (I) ,PLTGT (I) , PITTS (1) ,PLIRGS (I) , PLTGWC (I) , 

*PLTGW (I) ,PLTSHL (I) ,PLSHE{I) ,PLTPE (I) , PLTAE (f) ,PLTEVS (I) ,PIEGCM"(IV 
*,PAEGO(T) 

2574 CONTINUE 

2572 FORMAT{8X,I3,2X,4 (F7.2,2X) ,2 (P9.7,2X) ,2 (F5. 1 ,2X) , 4 (P8.5,2X) ,F5.2) 

T 






■ 6 iM 07410 
OLM 07420 
OLM 07430 

OIH Q 7440 

OLf 367450 " 

OLM 0746 Q 

0 LK 07470 * 

OLM 07480 


OLM07490 

OIW Q75QQ 

OLM075~10" 

OIMQ752Q 

OIH07530 

OLM0 7540 

“OLH07550 

ii-MPiseo 

OLiy!0757“0“ 

0LM07580 

'OIH07590“ 

OLK07600 

OIM07610"' 

q^LH0762J) 

'bl¥07630“' 
_ OIM07640 

“6“LM0765b" 

OLM07660 
OLM07670 
__ _piH07680 

,'2Xr OLM0769"d ’ 

PLM07700 

r ■ OIH077~10”' 

_ OLM07720 
■' ■■‘oLM07'73b" 

C) , OLM07740 

r(f)" dlH07750 
OIM07760 
■"OIK 07770 
5.2) OLM07780 
OIK 0779 O' 








w 




M 



















cd 

1 

u> 

Ln 


PNHH=NHH/1. 

TB'0T^1XW1?+ FV GlTC + A VT S)’7T: -2 ST 
7BOT=265. 


“TETFZf 


'VTOP=VBOI+50. 
MFITE (6,2564) 

F HAT (“A Y3‘ X , • Fv“ 


OlK07e00 

OLf-OTSIO 

0LM07820 

diMdVe'ao 

OLM.07840 


AVGT,RVTS,AVGK, !IVGTC,AVTSC,AVGWC 

GR TE'flP" = *';"BT(T;4 , 1 4X , » AV SURFACE TEMP =''“**‘,ET0‘V4 VOLM07'850 
115X,’AV GS WETNESS = • ,B1 0, 4 ,/, 3X ,• ASSIMILATION GB TEMP = ' ,B1 0.4 , OLM07860 

7 fT074“5X , • A'SSiM IL'A Tl (!) N 'G B W ET NE S OL M 07 8 7 0 




35 = ' ,E10.4) 


CALL G5DPLT (NHH,0. ,RNHH, VBOT , 7TOP ,PLTX , PLTGTC, PERIOD, 1 , 

GR T n?xT ' (yH Hrg' r7™F r7 VBor ;'g Yg F 7 'pr!rjrrn?i'ST';5TAR",T;Tr 


0) 


•cKnr 

WRITE (6,2580) 

72 5BTr^7)RMAT('20T7*PEHITnr^'TnfPUr"eT^5^^ = CTL GTf~~/7777) 

CALL GHDPLT(NHH,0.,SNHH,7EOT,VTOP , FLTX,PLTTS , POUND , 1 ,0) 

cArE-GRUP LT ( NMHVO .■,TrFrHTr7W)Tr,irTm?Tn?jr^^^ FiTrs7Tr ,"i ) 

WRITE (6,2595) 

TFyB—FDTT H AT C 2 JX7 ' pu iTFD - ' =-x ' NP ' U ' T ' ' T ir»' ;3X ' ; '» pi TT?^'~irrrci)Tnrrx!r^ 

CALL GRDPLT (NHH,0. ,RNHH, .000,250. ,PLTX, PLTSHL,STAR, 1 ,0) 

7?ry\) ,U10,PLTX,P LT9 TG ',P OU W70",Tr7 


2585 


CALL GRWLT(NH'ff77).»HNHb, 

CALL GRDPLT (NHH,0. ,RNHH,, 

■ WRITE (672’5B'5) 

FORMAT( 18X, • PLUS = T S K SDNS’ 

THTT 


000,250. ,PLTX,PLSHB, PLUS, 0,1) 


5X, ’STARS = NCAR SDNS WHSN 


*TH(9)/T - HG LS LESS THA ' r ~ T:0 ," 'TrNE TFBI E " CAS E ' 777777) 

CALL GRDPLT (NHH , 0 . , RNHH 05 ,0 . 20 , PLXX , PLSGCH, PER 10 D, 1 , 0) 

-”CAXL GRDPLT ( NHR7’’0ITR‘Nini7~U5VDTirO';PLTX , PXTPirrPLTrS7D‘,"D)‘ ~ ' 

CALL GHDPLT(NHH,0. ,RNHH,-. 05, 0.20, PLTX,PLTAE, STAR ,0,1) 

WICrTirTST2S9Dl 

259 0 FORMAT (20X, ’STARS = ACTUAL B7AP 

*(CM/HH) *,5X,»PLUS = POTENTIAL E7AP (CM/HH) ',/, 20X, ' PERIOD = ACTUAOLH08090 


OLH07880 
OL¥07890 
OLM07900 
'OL'Ma79TO 
OLH07920 
■ Ol'M07930 
OIM07940 
OLK07950 
OLK07960 
"OLMO7970 
OLM07980 
‘01M07990- 
OLM08000 
■■D'LM'OSDIO ■ 
OLM08020 
‘arNOSElD" 
OLM08040 
OLM080‘50 
OLM08060 
OLM0E070 
OLH08080 


*L EVAPOBATION FROM- INPUT TAPE (CM/HH) ». /////> 


OL«08J,00 

oiniosTio 

OL MO 8.120 
6lM08130 


CALL GRDPLT (NHH,0.,FNHH,. 000,1. 00, PLTX,PLTE7E, POUND, 1 ,0) 

CALL GRDPLT f NHH . 0 . . EN HH. 0 . . 2 . 5 . P LTX.PASGO .PERIOD. 0 . 1) 

WRITE (6 ,2609} 

2609 F0?«AT(20X. ’POUND = S7AP SCALING FACTOR *. 5X .» PERIOD = A SGC M/A BOLM ' OL M08 1 40., 
*,/////) OLMOSISO 

CALL GRDPLT f NHH, 0 . , RNHH ,0. , . 05 . PLTX .PITGWC .STAE.1 ,0) _0 L M 08.1 6 0 

OLMOSilO 
OLM08180 


CALL GRDPLT (NHH,0. ,RNHH,0. 
WRITE f6.26Q5) 


05,PLTX,PLIGW,PLUS,0, 1) 



•ge-e 


2605 FORMAT (20X, ’STAR = SHL1?(10) : TK' ,_5X, ’PLUS = SEL (9) ' //7) 

STOP ^ 

BND 

SUBROUTINE SDET fPI, DPI ,COSD , SIHD , EOT,DLON , ASSI H.IHH ,LAT , J D AY, 

INCELL, TOFDAY) 

DIMENSION LAT (46) ,DAYSPM (12) ,JMONTH (12) ,AMONTH (2,12) 


0LM08190 

OLM Q8200 

olmospTo 

OLM0 8220 

oLMoeYid 

OLM08240 


COMMON 

COMMON 


$ 

$ 


$ 


COMMON 


PHIS (08,14) ,FD (08,14) ,TOPOG (08 , 1 4 ) ,SINL(46) ,COSL (46) OLM08250 

OLM08260 

““‘“OLM 08^70 
^OM08280 
bLM08290 
OLM08300 


PC (8, 14) ,TSC(8,14) ,SHSC(8,14) ,GTC (8,14) ,GWC(8,14) , 
_U^(8,14,3) , VC (8,14,3) ,TC (8 , 14, 3) , SHC (8 , 1 4 , 3) ,PLK(10) 
SGC(8,14) ,US(08,14) ,VS(08,14) 


* 

♦ 


PLS(IO), TLE(IO), BH (10) ,TH(10) ,DSIG(10) ,SH1(10) , 
SHLE(IO) ,RF, (10) ,SH(8,14,3) ,PL(10) , 


OLH08310 

OLM08320 


C 

■■c* 


* 


COMMON 


GT (08,14) ,SHS (08,14) ,TL (10) , 

3LB5D 0 (08,14) , TS(08,14), GS (08, 14) ,ROK (08 ,14) , 
IWE (10) ,STGE (10) 7sIG(10) 


»RHQS,CZ 


3HS,DTS,CD,BDV,ZLN,DEAH,PSSVAP,PLANT,EVE,B6DT, 
H, DT5BBT,TEBM1 ,TBB M2 ,TES M3 , TEEM4 , TERMS, DSNOH,RAD TRM 


SNPAL, BAIN, OVER, WETS, WET, TAU, WMAG, EL E 

LOGICAL OCEAN,SNOW,ICF,LAHD,MIXWI, FROST 

LOGICAL GISSJK,TJCLA ,DCNI ,UCNII , UCN JOA , ASSIM, TINPT 




INTEGER ABONTH,DAYSPM,DAYSPY 
DATA AHONTH/’ JANUARY FEBRUARYMAPCH 


APEIL 


MAY 


JUNE 


JULY 


AUGUST SEPTEMB3R0CT032S NOVEMBERDECEMBER ’/ 

JJA PUT COOMENTS ON AMONTH 10/22/77 
DATA AMONTH/ 1 , 2 , 3 , 4 , 5', 6 , 7 , 8, 9, 1 0, 11 , 1 2 , 1 , 2 ,3 , 4 ,5 ,“6V7 r8T'9 ,T0, 

111,12/ 

DATA DAYSPM/31,28,31,30,31,30,31,31,30,31,30,31/ “ 

DATA IDAYO/302/,DAYS?Y/365/,LASTM/0/ 

INTEGES*2 SLBEDO " ’ “ 

PEAL*4 LAI 


0LH08330 
OLM08340 
OLM08350 
OLB08360 
OLHd837"d 
OLM08380 
0YMb’839b 
_pLM0840_0 
OLM08410 
0_LH08420 
■bLMO'8430 
OLM08440 
OLM08450 
OLH08460 
OLM084Y0 
_pLH08480 
OL Mbs 4 9 O' 
OLJ08500 
dLMO‘8510 
0LH08520 
"bLM'08 530 


INTFX (XTAU) =INT(XTAU*XINT+.5) 
j " ■ XINT-=600, 

C 

r” c 


OLM08540 

"OLM'bgSSO 

OLM08560 

OL'M'oes'^b 
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JSP=1 

OLH03580 

i 


JSB=i 

JKP=46 

OLM08590 

OLH08600 



JNB=46 

0L«0861‘0 



GT.i) GO TO 211 

OLM08620 

c 



.. — ^-jjj,og-g3o 



FJEQ=.5*(JSP+JNP) 

OLM08640 

i 


DLAT=180/ (JNP-JSP ) 

0“LM08e50 



DLAT=DL AT+PI/180. 

OLM08660 

i ■ 


DO 2U0 J=JSB,JNB 

“ - 0LM0'8F70 


200 

LAT ( J) =DLAT* (J-FJEQ) 

OLM08680 

S 


LAl: (JSP) =-. fj^DPj 

OLMOTeS'O' 

1 


LAT (JNP) =-LAT (JSP) 

OLM08700 



DO 210 J=JSB,JNB 

0”LM0871 0 

5 


TEMP=LAT(J) 

OLM08720 



SINL (J) =SIN (TEMP) 

OLKO8730 



COSL (J) =COS (TEMP) 

OLM08740 

C 


WHITE(6,30) J,LAT (J) ,SINL (J) ,COSL (J) 

OLM08750 

C 

30 

FORMAT (2X,‘ POE J =' , I4 ,5X , • LAT ( J) = *.E12.4,5X, 

OLM08760 

c 


1 ‘SINKJ) = » ,F12.4,5X, ‘COSL (J) = ',312.4) 

OLM 08770 


210 

CONTINUE 

OLM08780 

c 



OLm’68796 

c 



^QL^PJSOO 


211 

DT=3 00. 

. oXkbs'aio 



IF(ASSIM) DT=60.«60.*3. 

.. 0LH08820 


IDTHE=INTFX (DTHR) 


I24=INTFX (24.) 
NSTEP=.5-t-TAq/DTHE 
ITAU=NSTEP*IDTHB 
T&q=FLOAT figau) /XINT 


IDAY=ITaa/I24 
I OFDAY= FLOAT (ITAa--IDAY*I24) /XIN T 
ROT=2*PI*TOFDAY/24. “ 

IH=72 

DION =2.*PI/IK 


WRITE (5,212) DTHE,IDTHR,NSTEP,ITA.U,TAU,IDAY,TOFDAY,nOT,DLON 
212 FORHAT { 1 0X , » DTHR_= * ,E1 2. 4 ,5 X , < IDTHR = * ,_1 6 , 5X , « NSTEP = 


■I 


OLM08830 

OLH08840 

OLM08850 

OIM08860 

OLM0B870 

0LM08880 

OLH08890 

OLM08900 

OLHO8910 

OLH08920 

OLM08930 

OLM089UO 

OLM08950 

OLK08960 



1I6','5X, '■*TtAO ‘»'rI10,'5X/'"TAU = %F12.3,5X, 'XDAY= ',16^7 OLM08970 

210X, 'TOFnAY (Z HRS) = • ,21 2. 4 ,5X , • EOT = • ,F1 2. 4 ,5X , • DLON = • , S 12. 4 ,OIM 08 980 

3 32X, 'STEP = » ,l53 OLH08990 

C OLM09000 

C " ■ 'OLH090i0 

SOLS = 173. ^LM09020 

ArH12L=183. ■ ‘ OLM09030 

DECflAX=PI*(23.V180) OLM09040 

ECCN=.0178 " ' 6LM09050 

C OLM09060 

C CALCULATE CALBNDAP ' ‘ “ OLM0'90'70 

C OLK09080 

JYFAB=1 + {IDAY+IDAYO-1)/0AYSPY '"OLMOPOgO 

JDAY=IDAY+IDAYO-DAYSpy* (JYEAP-1) 01M09100 

J=0 ■ " . OLM09110 

DO 10 MONTH=1,12 OLM09120 

J=J + DAYSPM (MONTH) ‘ " ' " " OLH09130 

IF (JDAY.LE. J) GO TO 20 0LM09140 

TCF^NTINUF ■ ' ■ ' ■ OLM09'150 

20 JDATF=JDAY- (J-DAYSPB (MONTH) ) OLM09160 

JMONTH (1) =AMONTH (1 ,MONTH) OLM‘09'170 

JHONTH(2) =AMONTH(2,MONTH) OLM09180 

' ~ " OLM09190 

WRITE(5,213) JDAY,JDATE OLH09200 

... .. 

C OLM09220 


'0LM09230 


C 

”'TDTY=7n>T! 

SEASON= (FDAY-SOLS)/DAYSPY 

DT5T^(Fi)AX-APHEL) /DAYSPX 

DEC=DECHAX*COS (2. *PI»SEASON) 

P-3X Ti 5T 5:'T T' :T E' ccN ' ^r o s ' ' ( 2" ; ~ yp 1 y urs T) - ) - y»" 2 * 

SIND=SIN (DEC) 

-X'0'EnHCO'S~(WC) ^ - " 

IF(IHH,BQ. 1) WP.IIS(6,25) SIND,COSD 

2'5-FOH H AT"( /, 5 X, » si N D“==“"^ 727'2 . 4 VSX ^ ■~r ,ifT2T4T 

C 

■C ; E5 7 :D ~ l ' H ' "l TF P0GTJ7 I ) ''"' l 'S U ' TAX ' Ol inJir^'PTrrSTJ~^^^^ ' 


OLM09240 
OLH09250 
OLH09260 
■'OLMff9270 
OLH09280 
'“0LB09290 
OLH09300 
■‘OLM09310 
OLM09320 
"■ 0‘LM09‘330 
OLM09340 
OLH09350 
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■ ^ 



OLD509360 
■ OXN09370 
012109380 
"OX M 09 390 
OLM09400 
OLM094.10 
OLM09420 


' SDBBOUTINS SOEFWD (NLA Y, NL& YH 1 , JMP ,IN0 ,I M, J1 P, J BPP 1 , TINPT, ft SSI M) 6 lh0'9430 

• DIMENSION ZM(Q8.14) OLH094 40 

COMMON PHIS (08,14) ,FD (08, 14} ,TOPOG (08,14) ,SINL (46) ,COSL (46) OLM094ror 

COMMON OLMQ946 0 

$ PC (8,14) ,TSC (8,14) ,SHSC (8, 14) ,GIC (8, 14) ,GWC{8, 14) , 0LH094'70 

$ PC (8, 14, 3) , VC (8,14,3) ,g?C (8,14,3) ,SHC (8,14,3) ,PLKf 1 0) ,OL H09 480 

$ SGC (8,14) , OS (08,14) , VS (08, 14) OLM09490" 

COMMON OLM09 5PO 

* PLE(10|, TLS(10), P.H (10) ,TH (10) ,DSIG (10) ,SHL (10) , OLH095io” 

* S HIE (10) ,RH (10) ,SH (8, 14,3) ,P1 (10) , OLM09520 

* GT (08,14) ,SHS (08,14) ,TL (10) , " ‘ OLM095T6 

’ * SLBEDO (08,14) , TS(08,14), GW (08,1 4) ,R0K(08,14) , DLM09540 

* IWE (10) ,SIGE(10) ,SIG(10) ' ■ OLMOgs'SO' 

■ COMMON 0LM09 560 

* EHS, DTS, CD, EDV,Z1N, DRAW, PHEVftP, PLANT, EVE, B6DTi 'OLM09570 

*RHOS ,CZH,DTRBF.T,TEBH1 , TBR.M2 , TERM3 , TBBM4 ,TSB M5 , PHNOM , 3ADTRM , _ OJLHO^SBO 

* SNFAL, BAIN, OVER, WETS, NET, TAP, RMAG,ELE 'oLH69590 

LOGICAL OCB AN, SNOB, ICE, LAND, MIXg I,FBOST OLM09600 

LOGICAL GISSJK,UCLA,UCNI,UCNTI,uCNJJA,ASSiM,TINPT ' ' ' 6 lH09610 

I NTBG3R*2 SLBEDO OLM09620 

■ "c" ‘ . OLM09630' 

_C SURFACE WIND MAGNITODB OLM09640 

“““C ■ ■ OLM09 6'50" 

, C CMS; E SOB WIND OLH09660 

C ■ - - — - -Qj_f,Qgg=7Q- 

_ C__ NOTE FD (J,I) MUST B E CALCO LATED AT AUj_CELLS FOB 0LM09680 

C ■ WMAG CALCOLATION ’ OLM09 690 

IMIN=INC OLM09700 

• “ “ IMAX=IM ' “ " ■ “ ' 'OLK09710 

■ JHIN=J1P-1 OLM09720 

JMAX=JMPP1 ■ ■ OLK09730 

OLM09740 


C 

. —t:” 
. c 
!“c- 


SEF TSAN6 AND KAHN CODE LISTING PAGE A-23 

LINES HA 1 27 1 00“TO"‘HAT272T8no ' 

DID NOT INPUT INTO MODEL AT THIS TIME 10/29/77 


BETUBN 

■END 


C 



CQSF1=(1,-SIG (NLAYHI) )/{SIG (NLAY)-SIG'(NLAYMir) 

C0FF 2= (SIG (NLAY) - 1 ) / (SZG (NLA Y) -SIG (NLAYMl) ) 

DO 10 I=I«IN,I!1AX ■ ' 

DO 10 JBAX 

IF(TINPT) GO TO 10 

US(J,I) =C'03F1*0C (J,I,3) +C0EP2*0C(J,I,2) 

VS(J,I) =C0EF1*VC(J,I,3) +COEF2^?C (J,I,2) - 

10 FD (J,I) =DS(J,I) **2 +VS(J,I)**2 
RYTOBN • 

END 

SOBSOUTINE EITB (X8SS,I, J,IHH,HOUR,IPE, 
*UCLA,DCHI,UCNJJA, UCNII,I2C»J2C, IHAX,JMAX,i aiNUS >J MI NOS « MONTH) 
DIMENSION ~ 

* ZM (08,14) , CELL (08,14) 

COM MON PHIS (08,14) ,FD(08,1'4) ,TOPOG (88 , 1 4) ,SnHL {hq-jmsTlUJ 

COMMON 


COMMON 


PC (8, 14) , ISC (8/14) ,SHSC(6,14) ,GTC (S;i4) ,GHC (8',~1 4)', 
UC(8,14,3) , VC (8,14,3) ,TC (8,14,3) ,SHC (3,14,3) ,PLK(10) 
SGC(8,14) ,DS (08,14) , VS (08, 14) " 


■,SHL (TU) , 


PLE ( 1 U ) , TLE (10), HH(10), TH (TO) , DS I G ( 10) ,SHL”rTU^ 
SHLE(10) ,HE(10) ,SH(8,14,3) ,PL(10) , 

■gt"( 08,14) ,SHS (08,14 ) ', T1 ( TTiy/ 

SLBSDO (08, 14) , TS(08,14), GW (08 , 1 4) ,E0 K (08 , 14) , 
TWH (T0) ,SIGB(T0') ,SIG (TDT“ 


COMMON 


^ /'dTS, CD ,EdV , 2tN’ , dS AT,W svl'F,Ti:iNr,W1?/l;WT , 

♦EHpS,CZH,DTRBFT,TEEHl ,TERM2 ,TB5 M3 , TEE M4 , TERMS , DBNOM,E ADTRM, 

* SN F ALTR'SrN , 0 V EE , W E TP. , WE T , T A U , W H A G ,"EL E 

LOGICAL 0 CEAN, SNOW, ICE, LAND, MIX WI,FEOST 
LOGICAL UCLA,OCNI,UCNII,DCNJJA,ASSIM,TINPI 

LOGICAL PRIN,PRISKV,PESDEV,PBGST,PEMCSV 

INTEGER PARM10,PAFH1 1,PARH12,PAPM13,PABM50,PASM51 ,PARM52,PARM53 
INTEGER *2 SLBEDO 


CMS; E SUB PEINT 


LOGICAL PRINT CONTROL VARIABLES 
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'I'd.- 



PRIN=.TROS. 



OLM10 140 


paxsHV-.TBas. 

PSSORV=.a?BljR. 


O.LM10150 

OLM10160 


PSSR"t=. TBUB. ' 


OLM10170 


PEMGEV=.TRaB. 


OLM10180 

. c 



OLM10190 

c 



OLH10200 

c 

PRINT CONTROL 

VARIABLES 

OLH10210 

c 

' * 


OLM10220 


PPRH10=f2C 


OLM10230 

• 

PARMi i;=J2C 


OLM10240 


PABH12=IHH 




PARM13=1 


OLM10260 


PAR MS 0=1 


OLM10270 


PABH51=J 


OLH10280 


PABa52=IHH 


OLH10290 


PAP.M53=,Tpp 


OLH10300 

c' 



OLM10310 

c 

INITIALIZED VARIABLES 

OLM10320 

! c 



0LM10330 

1 

X S. 

IFflPR.NE.I.OE.'.NOT.PRIN) GO TO 305 


OLM10340 

.! C 



blH 10350 

• c 

INDEX TOR CELL 


OLM10360 


Tl=jMINUS-1; . ' 


■ OLM10370 


DO '1217 J=1,JHAX 

, 

OLM10380 


T1=T1+1 

' 

OLM10390 


T2= (IMINpS-1) /100. 


OLM10400 


D"onr2TTT=T^ I MA“X • 


“■6lM10410 


T2=T2 + . 01 


OLM10420 


CELL (J,I) =T1+T2 ■ 


'■■olm‘IO'430 

1217 

CONTINUE- 


OLB10440 

. C 



bi tiro’s 5*0 


WHITE (6,234) 


OLH10460 

■*'”2"3Tr 

FORMAT ( 7/7H3x7'GISS INDEX FOB CELL 

(J.I) •) 

btH10470 

. 

WEITE(6,90) (I,I=1,IMAX) 


OLM10480 

!" 

1 

J=JMAX+ 1 


OLnOS’^O’ 

\ 

1 

DO 1216 JJ=1,JMAX 


OLM10500 




' 'OL'M'IOSrO' 

1 

\ 

WRITE (6,233) J, (CELL (J,I) ,I=1,IMAX) 


OLM10520 





' ■ r2'iT 
c 

"CONITNUT ■ ~ 


201 

■TSITFTb";20TJ 

FOB«AT(5X,//,33X, ’SURFACE ALBEDO FOR JULY, DSNABS = 

I527ALBD’) 

* WHITE (57^)1) 

J=JMAX+1 

DO 'I203"UJ=1 
J=J-1 

1203 

""WE ITS 1572751 J, (SLBEDO(J,i) ,i=1,IMAX) 

CONTINUE 


2tTT 

C 

“FOKBAl’l 08 (2X, 13, 14 [Wl^rZ'rZTn TT 


;■ 207 

WRITE (6;207) ‘ 

FORMAT ( //, 4 OX, ’SURFACE ELEVATION (GSOPOTENTIAL HT) 

’) 


WK'.TE (6 ,9 1) (i ,1= 1 , iMAXJ 
J=JMAX+1 


i 

UO 1204 Jj=1,JMAx 
J=J-1 

Ss 

^ 5_ .. 


(6^2 1 8) J ( PHIS'XT,' i)"VT='l #' 


1204 CONTINUE 


3i 

Sl 


■■OLNI'0530 

OLM10540 

0L'&T0550 

OLM10560 

■'OtHl0570 

OLH10580 

•OT>n05?'0‘ 

OLM10600 

"oxHTcr'ei'o 

OLR10620 

■O'LFr’TO'6'30 

OLM10640 

o'L’TraeBo 

OLM10660 

'■‘or«ro'67o 

OLK10680 

OrM1069O 

OLH10700 

'orwiuTio 

OLR10720 

OLM10730 



write f6. 2091 MONTH 

lO— 3n_- 

OLM10740 

209 

FORMAT ( //,44X,’TOPOG ARRAY FOR MONTH = ’,12) 

ifeS? 

jC 

OLH10750 


WRITS(6.911 (I.I=1.IMAX) 

OLM1Q760 


J=JMAX+1 

DO 1206 JJ=1 . UMAX 

■^Sj 

biM 10770 
OLM10780 


J=J“1 


C)LHlb790 


WRITE (6,210) J, (TOPOG(U,I) ,I=l',IMAX) 


OLM10800 

TzoT 

CONTINUE 


OLM '10810 

210 

PORHAT( 08(2X.I3, 1X.14(F7.0,1X) )) 


OLH10820 

c 

WBITS(6,212) 


OLM 10830 
. „ , OLM10840 

__ 

FORMAT( //,43X,’ INITIAL SURFACE TEMPERATURE (K) ’) 


OLM 10850 


WRITE(6,90) (I,I=1,IMAX) 


OLK10860 


J=JMAX+1 


OLH10B70 


DO 1207 UJ=1,JMAX 


OLM10880 


J=J-1 


bL"Kl689b 


WRITE (6 ,233) J, (TSC (J,I) .1=1 ,IHAX) 


OLM10900 

Y27T 

continue 


OLH109lb 







WRITE (6,214) 

214 FOBMAT( //,43X, « IHITI&L GROUND TEMPERATURE (K) « ) 
WHITE(6,90) (I,I=1,IMAX) 

J=JBAX+ 1 ; 

DO 1208 JJ=1,JHAX 

J=J- 1 

WRITS (6,233). J, (GTC(J,I) ,I=1,IMAX) 

1208 CONTINUE ^ 


WHITS (6,216) 

216 FOEMAT( //,43X, • INITIAL SURFACE FEESSURE (MB) » ) ' 

WRITE(6,91) (1,1=1 ,7.HAX) ; 

J=JMAX+1 

DO 1209 JJ=1,JMAX 

J=J~1 

W PITE(6,218) J, (PC (J,I) ,I=1,IBAX) 

1209 CONTINUE 

218 FOBaAT( 8 (2X,I3, IX, 14 (F7. 1, IX) )) 

C 

WRITE(6,205) ^ 

205 FORHAT(/,2X,38X, ’INITIAL GROUND WETNESS MATRIX (GW C) » ) 

HRITE(6,91) (I,r=1,IMAX) ; 

90 F0PMAT(/,3X, »J I’ , 14(13, 5X)) 

■?1 FORMAT (/,3X,» J/T» ,14 (13, 5X) ) 

""J=JHAX+1 

. DO 1213 JJ=1,JHAX 

J=J-1 

WRITS (6,203) J, (GWC (J,I) ,I=1,IMAX) 

■' 1210 CONTINUE 

200 FOHMAT( 08 (2X ,13 , 14 (F7. 4 , 1X) )) 

c " ' ■ 

WRITE (6,237) 

237 FORMAT (//,45X, ’SUEFACS SPECIFIC ATION’ ,/,1 OX, • 1 OCR ANT“’,3r, 
1'2: LOAMY SOIL, ’ , 3X, ' 3 iSANDY SOIL, ’ ,4X , *4 : ’) 

ITsITET^T^ (I,I=1,IMAX) 

J=JHAX+ 1 

DO"TJ2'2 JJ=1,JMAX ■ 

J=J-1 


OLM10920 
" OLM10930 
■ 0LM10940 
OL Ml 0950" 
OLM10960 
OLM10970 
OLM10980 
''ol¥10 996'' 
_ O LM1 1000 

"olmiToTo 

O LM1 1020 
‘"'‘OLMfl030’" 
OJiMlI 040 
OLBIIOSO" 
OLM11060 
- OLMiY670 
_OLMl 1J)80 
OLMIIODO” 
OJbMlIYOO 
'OLMlTliO' 
OL H1112Q 
OLMl'li'30 
OLM11 140 
OLMIIISO'" 
OLMl'l 160 
"OLMI1 170’ 
OLM11 180 
OL'HI 1T90" 
OLM112O0 
0LM112"16' 
0LM11220 
‘OLMl 1T3’0 
OLM11240 
OLM11250 
OLMI 1260 
■■ OL'Mll270 
OLHH280 
“ OLM11290 
OLH1 1300 
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U 


1222 


WHIT3 (6,253) J, (T^OK (J,I)7^#IK1'.XT 


olHTsio 


. 

WRITE(6,236) 

_ “"-’0L«'11330 

OIM11340 

■~2'3ir 

“Fok1iat(//,' 34 X,»tranfOrmbd ground wetness 

MATRIX (GWT^^-.T SCALE) OOL"H1'1350 


W?ITE(6,91) (I,I=1,IMAX) 

0LM11360 


J=JMAXTT“ 

OLM11370 


DO 1220 JJ=1,JHAX 

OLM11380 


J=J-1 

OLM11390 


WRITE(6,200) J, (GW(J,I) ,I=1,IMAX) 

OLM11400 

1220 

CONTINUE 

6LMl14io' 

c 


OLM11420 


WRITE (6,235) 

"OLM11430 

235 

FORMAT (//,45X, ‘WMAG COMPONENT : FD(a,I)') 

0LM11440 


WRITE(6,90) (I,I=1,IMAX) 
a=JHAX+ 1 


DO 1221 JJ=1,JHAX 
J=3-1 


WRITE (6,233) J, (PD (J,I) ,1=1 ,IMAX) 
1221 CONTINUE 


IFdHH. EQ.1) WRITE (6,240) 


IF(IHH.GT.I) WRITE(6.,221) 

2U'0 FORMAT {//,45X, » S OFFACS SPECIFIC HOMIDITX*) 
'221 POB«AT( //,47X',-*NET RADI ATION "(LY/DAY) ') 
«PITE(6,91) (I,I=1,IMAX) 

‘j=JMAX+1 ” 

DO 1211 JJ=1,JMAX 


J=J-1 
I F(IHH.FQ.I ) 
IF(IHH.GT.I)" 
1211 CONTINUE 


WRITE (6,22 3) 
WRITE (6; 21 8) 


J, (SHSC (J,I) ,1= 1 , IMAjy _ 
J, (SHSC (J, I) ,i=T,Tmax) 


LL=0 

DO 229 L=1,3 


LL=L+6 


WRITE (6,222) LL 

222 FQHaAT{ //,45X,*X COMP O F VEL AT LEVEL*, 12) 
WRI TE (6,96) (Id=1,IMAX) 


■ OLM11450 
OEMliaeO 
OLMliiiiO 

. PLM11480 

biiiiiugo 

, OLM11500 
'“ol'niiBio 
OLM11520 
OLfil 1530 
OL»11540 
OLM11550 
OLN11560 
Ol«1157d 
OLM11580 
OLH11590 
OLM11 600 
OLM116 10 
OLH1 1620 
OLHli630 
O1N11640 
OLM11650 
OLH11660 
OL Mil 670“ 
OLM11680 
0LA!11M6 



3l* 


J=JHAX+1 


w 

-p- 

Oi 


DO 1212 JJ=1,JMAX 
J=J-1 


OLM11700 

OLK11710 

0LH11720 


1212 

W?ITE{6,233) 

CONTINUE 

J, (UC (J,I,L> ,1=1 ,IMAXj 

■ OLH11730 
OLH11 740 

2JT 

C 

FORMAT ( 

08(1X‘,I3,1X, 14(.F7.2,1X) )) 

OLMfiySO 

OIH11760 


WRITE(672'24) LL 


OLM11770 


224 FORMAT( //,45X,'Y COKP OF VEL AT LEVEL’ , 12) 
MRITETFTW (r,I=1,IMAX) 

J=JHAX+1 

DO 1213 JJ=1,J!1AX 
J=J-1 


1213 CONTINUE 


j, (VC (0 ;T7lir7l = iTXlA'xr 


J, (TC (J,I,L) ,1=1 ,IHAX) 


WHITE (6,225) LL 
WRITE(6,90) (I,I=1,IHAX) 

j=Tfnrx"i 

DO 1214 JJ=1,JMAX 

J“=U=1 

WRITE (6,233) 

'•^TTU'XUNTXNirr'' 

C 

226 FORKAT(^ //,45X,’SPEC HUMIDITY AT LEVEL’, 12) 

5rRTTET6T9l^' ~(x','l=r,‘XMTT' 

J=JMAX+1 

“■’■■D'D“T215'-J'J=T,'J Ml X 

J=J-1 

-FEITE (6,223) 

1215 CONTINUE 
AT ( 

229 CONTINUE 
C 

305 CONTINUE 


"J”rs7ic7u~,“r,‘: 


ubttt;t 371'X7"1 .■'-T( F7 747TX) — n 


c 

£ 


S URFAC E _V A,E.T AB LXS_ 


OLM11780 
0Lrfl1‘790 
OLM11800 
'O'LH'irSlO 
0LM11820 
OLM1T830 
OLM11R40 
OLM11850 
OLM11860 
01M11870 
OLM11880 
OLM11890 
OLH11900 
'OLK'11910 
OLM11920 
'orriio3o 
OLH11940 
0LK11950 
OLM11960 
01'M11970 
OIMlI 980 
OLH11990 
OLM12000 
OLM12010 
OLM12020 
OLK12030 
OLM12040 
OLB12050 
OLM12060 
OLM12070 
OLH12080. 



IFfIP3.gI^.2.0P..N0T.PEISMV) GO TO 395 ; 

1 FARM 13=2 

IF fPARfnO.EO.PARa50.aND.PARMl1.5Q.PAHM51.ftND.PAR«12.B( 
AND.PARH 13.BQ.PARH53) 

*GO TO 32Q ; 

GO TO 395 

) CONTINUE 

WPITE(6,490) 

) FORMAT (/,Q4X, 07X, »DTS',07X, «ELS«,08X, 'EPy*, 

*05X, ’SHLFCIO) ',03X, •SHL(9)*, 04X, *DRAW« ,05X, 'TLE (10) 
WRITE (6,500) BHS.DTS, E LE ,ED V, ZLM, SHLE (10) ,SHL (9) ,DRA H, 
FORMAT (2X,09 (Ffi.4 ,2X) ) 

CONTINUE 


PARMf 


0JX^.._ 

*) 

TLE 


GROUND TERMPBEAIUPS CALCULATION VARIABLES 


IF (IPE. HF.4 .OS. .NOT.PEGRT) GO TO 605 

. PAR Ml 3=4 

IF(PARH10,SQ.PARM50.ANP,PARH11 ♦EQ.PAEM51.AND.PARH12. 
♦ . AND.PAPM13.F.Q.PAEM53) 

*GO TO 510 

GO TO 605 

I CONTINUE ^ 

WRITE (6',595) 

i F03aAT(6X , *EHOS*,07X, *DTS *,Q6X,«CZH», 

"♦»WMAG *,06X, 'TEEM2', 05X,*TFRK3*, 06X, 'TESM4»,06X 

« , 06 X , » DENOM* ,06X, »RADTRa* ,06X, »GT «) • 

"WPITE (6 ,600) BHOS,DTS ,CZH, ‘ WMAG",TSRM2 , TERH3', 
♦TL(9) ,DENOM,RADTRM,GT (0,1) 

> FOSaAT( IX, 1 1 (F101 3, IX) ) 

J CONTINUE 


EQ.PAEK! 


7’TL(9)' 
TERM4 , 


MOISTURE CONTINUITY EQ. VARIABLES 


IF (IPR.NE.5.0R..N0T. P5HCE7) GO TO 695 _ 

i PAR Ml 3= 5 

IF(PAEMl.0,3Q.PARM50.AND.PAEMl1.BQ.PAEM51.AND.PARm2.EQ.PARH' 
♦ . AND.PARM13.EQ.PARM53) 


ORJGffilAL PAGE IS 
OF POOR QUALITY 





w 

I 

*p> 

•vj 


♦GO TO 610 01M12480 

GO TO 695 """ OLH1'24'9"0" 

610 COHTINIJE OLB12500 

WRITE (6, 690) ' " 01Ml25"iO 

690 FORMAT (5X, » SNESL' , 05X, ‘SAIN * ,0 6X, * OVER • , 06X, • PL ANT * , OLM12520 

♦ Oax, *P0T-S' ,04X,*ACTUAL-E» ,04X, »HSTS*,06X, »»ET“", OLM1253'0 

♦7 X, *RH (9) « ,6X^*PL(9) 6X,»SHL(9) ») 0LM12540 

WEiTE(6,700) SNPAL, RAIN, OVER, PLANT, PSEVAP,E6Df,WETH,W2T7 " ’ 6 l«12550 

♦RH{9) ,PL (9) ,SHL (9) OLil12560 

'““700 F0PaAT(2X,11 (r9.4,lX) ,//) ' " D1K12570 

695 CONTINUE OLM12580 

“5~‘ y/BXEC LDl ■' OLift259’0 ' 

C /DD 008 DSNAME=DJJA ‘ OLM12600 

C~" /EXEC ♦ " ' ■ 0'LKT2610' 

C OLK12620 

'I=PARM5T5 ■ > -- - - • oLfn‘2630 

J=PABM51 OLM12640 

^ OiH"126'5'0 

RETURN OLM12660 

— OLM12670 



